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Preface 
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the  choosing  of  problem  areas  which  could  be  profitable  investigated. 

I  also  want  to  express  ny  thanks  to  I-Iiss  Gayle  Mortensen  for  her 
help  in  compiling  the  data  presented  in  Appendix  B,  and  to  lirs, 
Henrietta  Thomas  for  accomplishing  the  tedious  task  of  typing  the 
decay  schemes. 

The  support  given  to  me  by  the  Air  Force  Weapons  laboratory, 
the  Research  Division,  and  Biophysics  Branch  in  jnrticular  was 
excellent  and  materially  aided  in  the  completion  of  the  research 
effort. 
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Abstract 

This  thesis  develops  a  method  by  which  the  time  dependent  beta 
spectrum  from  four  different  types  of  nuclear  and  thermonuclear 
fission  can  be  predicted.  From  the  mass  chain  yield  and  the  Equal 
Charge  Displacement  theory,  the  initial  isotopic  concentrations  are 
predicted.  The  time  dependent  relationships  vdthin  each  mass  chain 
are  solved  exactly.  The  shape  of  the  beta  decciy  spectra  for  each 
isotope  is  predicted  by  use  of  the  Femd.  theory.  An  estimation 
tec})nique  is  used  to  predict  the  spectra  from  isotopes  vdth  unknown 
decay  schemes.  Attempts  are  made  to  predict  the  percent  of  the  total 
available  decay  energy  which  is  accounted  for  by  beta  decay,  for 
those  isotopes  with  unknown  decay  schemes,  Tnese  procedures  are 
combined  in  a  set  of  computer  progr.ams  written  for  the  IBM  704A, 

The  spectra  predicted  using  the  programs  are  compared  with  West's 
experimental  results  and  agree  within  a  factor  of  1,4, 
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I,  Introduction 


Purpose  and  Scope 

The  purpose  of  the  research  reported  in  this  paper  is  to  con¬ 
tinue  development  of  a  method  by  which  the  time  dependence  of  the 
beta  spectrum  resulting  from  a  nuclear  or  thermonuclear  detonation 
can  be  predicted.  This  effort  was  started  by  Lt,  Robert  L,  Stovall, 
AFIT  GSP-64,  and  reported  on  in  his  thesis  (Ref  17).  Since  Lt. 
Stovall's  thesis  is  classified  SiSCHIJr,  in  accordcuice  with  AFR  205-1, 
and  its  availability  is  therefore  limited,  enough  of  his  work  is 
included  in  this  report  to  assure  continuity  of  thought.  Because  of 
the  complexity  and  number  of  the  processes  involved  in  these  detona¬ 
tions,  this  prediction  method  necessitates  the  use  of  a  digital  com¬ 
puter,  The  problem  encompasses  many  areas  where  experimental 
evidence  is  not  available;  in  these  areas  the  best  available  theo¬ 
retically  predicted  information  is  used,  and  where  no  theoretical 
predictions  are  possible,  the  author  has  made  approximations  which  he 
believes  are  logical  and  appropriate.  In  order  to  understand  better 
the  calculations  necessary  to  make  the  predictions  for  any  of  the 
reactions  of  interest  (U-235  thermal  neutron  fission;  U-235  fission 
spectrum  neutron  fission;  Pu-239  thermal  neutron  fission;  thermo¬ 
nuclear  neutron  fission),  a  brief  review  of  the  typical  reactions  and 
processes  which  are  taking  place  will  be  helpful, 

'.^'hen  a  compound  nucleus  fissions,  it  splits  into  two  fission 
fragments  and  emits  a  number  of  neutrons,  depending  on  the  nucleus 
undergoing  fission  and  the  energy  of  the  neutron  causing  the  fission. 
After  a  large  number  of  such  fission  events,  a  plot  showing  the 
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statistical  assembly  of  the  fission  products  according  to  their  masses 
exliibits  the  familiar  "double  humped"  shape.  This  curve  indicates  the 
percent  of  the  fission  product  assembly  characterized  by  each  mass 
number.  In  many  cases  these  curves  have  been  determined  experi¬ 
mentally,  In  tills  study  the  mass  chains  with  atomic  weights  of  77 
through  159  will  be  considered.  In  each  mass  chain  there  are  several 
different  charge  values  (atomic  number)  which  can  be  formed.  They 
are  formed  in  vai^’^ing  amounts  v/hich  can  be  predicted  by  a  theoretical 
approach;  each  of  the  members  of  each  of  the  chains  then  decays 
according  to  the  radioactive  decay  lav/.  Since  the  fission  products 
are  in  general  neutron-heavy,  these  decays  are  accomplished  by  the 
emission  of  beta  particles.  Each  of  the  members  of  the  mass  chain 
•fiay  emit  one  or  several  beta  particles  which  leave  the  resulting 
nucleus  in  a  variety  of  excited  states  which  then  decay  to  the  ground 
state  by  the  subsequent  emission  of  one  or  more  gamma  rays.  In  beta 
decay  the  beta  particles  are  observed  to  exhibit  a  spectrum  of 
energies  characterized  by  a  well  defined  upper  limit,  referred  to  as 
the  end  point  energy.  This  beta  spectrum  for  each  individual  beta 
decay  is  predictable  by  use  of  the  Fermi  theory  of  beta  decay. 

Therefore,  to  predict  the  gross  beta  spectrum  of  a  mixture  of 
fission  products  at  any  time  of  interest  follov/ing  the  fission  event, 
it  is  obvious  that  a  great  deal  of  information  must  be  available 
experimentally  or  through  theoretical  means.  This  information  con¬ 
sists  of  first,  the  percent  of  the  gross  r.±xture  of  fission  fragments 
’.;hich  are  contained  in  each  of  the  mass  chains;  second,  the  percent 
of  the  total  chain  yield  ’.diich  is  accounted  for  ir.mediately  after 
fission  by  each  member  of  the  chain;  third,  the  structure  of  each 
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decay  chain  including  all  branching  ratios;  fourth,  the  decay  con¬ 
stants  of  each  member  of  the  chain;  and  fifth,  the  number  and  energy 
of  the  different  betas  by  which  each  isotope  decays  and  the  end  point 
energies  of  these  beta  decays. 

Previous  Work  in  This  Area 

The  first  effort  to  attempt  to  solve  this  problem  took  place  in 
195B  (Ref  14),  The  approach  used  consisted  of  taking  the  then  avail¬ 
able  data  from  BoUes  and  Ballou  (Ref  l)  concerning  the  time  dependent 
activity  of  the  fission  products  from  thermal  neutron  fission  of  U-235 
and  then  using  the  current  information  concerning  the  decay  schemes 
and  half  lives  for  the  fission  products  and  solving  the  decay  equa¬ 
tions  to  give  the  spectra.  The  time  period  over  which  the  calcula¬ 
tions  were  made  was  from  31»2  minutes  to  119  years.  No  apparent 
effort  was  nade  to  predict  the  spectra  at  shorter  times.  There  also 
was  no  attempt  nade  to  theoretically  predict  or  estimate  any  of  the 
unknown  data. 

The  next  attempt  to  solve  this  problem  occurred  in  I96I  (Ref  12), 
This  effort  contained  a  more  thorough  attack  on  the  experimentally 
unknown  data.  In  treating  the  isotopes  with  unknown  decay  schemes, 
however,  the  authors  considered  that  all  of  the  available  energy  was 
taken  by  the  beta  particle  and  none  by  the  gamma  emission  from  excited 
states  of  the  resultant  nucleus.  This  approach  was  taken  in  order  to 
arrive  at  an  upper  limit  for  the  high  energy  portion  of  the  spectrum. 
The  results  of  this  effort  were  presented  as  graphs  of  total  beta 
activity  per  fission  versus  time,  the  percent  of  total  beta  activity 
versus  beta  energy,  and  the  number  of  beta  rays  per  fission  per 
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energy  interval. 

The  most  recent  effort  in  this  area  is  presently  being  conducted 
by  the  b.  3,  Naval  Radiological  Defense  Laboratory  (Ref  o).  The  irain 
purpose  of  the  work  is  to  develop  a  means  to  predict  the  radiological 
properties  from  fractionated  nuclear  or  thermonuclear  v/eapon  debris. 

In  the  course  of  developing  this  prediction  model,  it  v^as  necessary 
to  develop  a  prediction  method  by  wiiich  the  abundances  and  activities 
of  all  of  the  fission  product  nuclides  could  be  calculated,  A  great 
many  of  the  assumptions  made  by  Stovall  and  in  this  thesis  are  identi¬ 
cal  with  those  made  by  the  researchers  at  the  U.3.N.R.D.L.,  especially 
in  the  areas  of  chain  and  isomeric  yield  determination,  partition  of 
yield  among  isomeric  pairs,  and  solution  of  the  decay  equations.  It 
will  be  interesting  to  compare  the  solutions  arrived  at  by  both  groups 
of  researchers. 

Results  Obtained 

Briefly  this  paper  has  taken  the  prediction  method  developed 
earlier  (Ref  17)  and  improved  its  accuracy,  efficiency,  and  versa¬ 
tility,  Its  versatility  has  been  enhanced  by  the  inclusion  of  a 
fourth  type  of  fission;  U-235  fission  spectrum  neutron  fission.  Its 
efficiency  has  been  improved  by  the  inclusion  of  the  computer  programs 
of  the  exact  solutions  to  the  differential  equations  describing  the 
isotopic  decays.  The  efficiency  has  also  been  improved  by  the 
development  of  a  method  by  which  the  thennonuclear  fission  case  can 
be  treated  in  a  completely  unclassified  manner,  thereby  allowing  the 
widest  dissemination  of  the  prediction  method,  its  results  and  con¬ 
clusions,  As  a  result  of  this  wide  dissemination  it  should  be 
possible  to  more  easily  obtain  independent  critical  analyses  of  the 
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ijietliod  and  its  results.  The  accui’acy  of  the  method  has  been  Liiproved 
by  the  incoi’poration  of  new  data  concerning  the  half  lives,  the  decay 
paths,  and  the  beta  decay  schemes  of  the  individual  isotopes.  In 
addition  to  these  ijnprovements,  this  paper  also  reports  on  a  method 
which,  although  it  has  not  yet  yielded  any  useful  results,  should 
allow  the  determination  of  the  percent  of  the  total  energy  available 
in  an  isotopic  decay  which  is  accounted  for  by  beta  decays. 

Report  Presentation 

The  reminder  of  this  report  will  discuss  the  topic  areas  in  the 
logical  sequence  indicated  in  the  preceding  discussion.  Chapter  II 
will  consider  the  experimental  evidence  available,  and  the  theore¬ 
tical  techniques  used,  to  predict  the  distribution  of  the  fission 
products  according  to  the  nass  chains  and  also  the  distribution  of 
the  fission  products  formed  within  each  mass  chain.  Chapter  III  will 
discuss  the  characteristics  of  the  decay  cliains  and  how  their  time 
dependence  is  solved,  Ciiapter  IV  will  present  a  discussion  of  the 
theory  by  'which  the  individual  spectrum  associated  with  each  beta 
decay  is  computed,  for  the  isotopes  ’.vith  known  decay  schemes,  and 
v/ill  present  the  estimation  technique  '*ath  which  those  isotopes  with 
unknown  decay  schem.es  will  be  dealt.  Chapter  V  will  then  describe  the 
procedures  by  which  all  of  these  calculations  are  brought  together  by 
the  computer  codes  which  have  been  developed  and  hoi^  the  end  pixiduct, 
the  time  dependence  of  the  beta  spectrum,  is  obtained.  Chapter  VI 
will  present  the  results  obtained  in  this  study  and  will  discuss  the 
limitations  on  these  results.  This  chapter  will  also  present  a 
discussion  of  the  utility  of  tiiis  research,  and  its  prediction 
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technique,  to  the  Air  Force  and  to  the  rest  of  the  scientific 
community  and  also  suggest  areas  in  which  further  effort  on  this 
problan  may  be  fruitfully  directed.  The  mass  chain  yields  for  the 
four  types  of  fission  considered  in  this  report  are  contained  in 
Appendix  A,  Appendix  B  contains  all  of  the  data  pertinent  to  the 
decay  schemes  for  mass  chains  numbers  77  through  159.  In  appendix  C 
are  foimd  the  PORTRAN  listings,  the  definitions  of  all  variables  used, 
and  the  flow  charts  for  the  computer  codes  developed  for  this  report. 
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II.  Chain  and  Isotopic  Yields 


Ikss  Chain  Yields 

One  of  the  basic  pieces  of  inforiTiation  necessary  for  the  solution 
of  this  problem  is  the  percent  of  the  total  number  of  fission  frag¬ 
ments  which  are  formed  with  each  atomic  weight.  This  data  is  con¬ 
tained  in  what  are  commonly  known  as  mss  yield  curves.  The  shape  of 
the  mass  yield  curve  varies  as  a  function  of  the  fissioning  nucleus 
and  the  energy  of  the  neutron  causing  the  fission.  In  Fig.  1  are 
presented  the  mss  yield  ciorves  for  (1)  thermal  neutron  fission  of 
U-235;  (2)  fission  spectrum  neutron  fission  of  U-235j  (3)  therml 
neutron  fission  of  Pu-239j  and  (4)  therraonuclear  neutron  fission, 

(Ref  18:17-37)  (Ref  6:7-15)  The  values  for  the  mss  chain  yields  per 
10,000  fissions  for  each  of  the  cases  are  presented  in  Appendix  A, 

From  the  shape  of  the  mass  yield  curves  for  the  four  different  cases, 
it  can  be  seen  that  the  trough  area  rises  and  the  twin  peaks  v/iden  as 
the  energy  of  the  fission  causing  neutron  increases.  It  can  also  be 
observed  that  the  curves  are  approximtely  symmetrical  about  a  mss 
number  equal  to  approximtely  one  half  of  the  mass  number  of  the  com¬ 
pound  nucleus  (the  sum  of  the  mss  number  of  the  nucleus  undergoing 
fission  plus  the  mass  of  the  fission  causing  pi-ojectile)  less  the  mass 
of  the  prompt  neutrons  emitted  in  the  fission  process. 

Isotopic  Yields 

Once  the  mass  chain  yields  are  knovai,  uhe  next  problem  area  is 
the  apportionment  of  tiiis  yield  among  the  various  isotopes  which  make 
up  this  mass  chain.  The  theory  which  vn.ll  be  used  to  predict  the 
isotopic  yields  in  this  report  is  the  one  proposed  by  Glendenin  which 
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is  be  sed  on  the  postulate  of  equal  charge  displacejnent  (J,  C.  D) 

(Ref  95439-515),  (Ref  4:646-63).  Its  initial  assumption  is  that  the 
most  probable  charge  formed  in  fission  (Zp)  for  any  atomic  i.iass  is 
equally  displaced  from  the  charge  of  maximum  stability  (Z;^)  for  both 
the  light  and  the  heavy  fragments. 


Heavy  fragment  (l) 


Light  fragment 


The  particulars  of  this  theoiy  are  described  in  completeness  in 
an  earlier  report  (Ref  17).  The  result  of  this  theory,  as  modified 
by  VJahl  (Ref  4:649),  is  the  prediction  of  the  fraction  of  the  iiass 
chain  yield  which  is  initially  forraed  as  each  member  of  the  mass 
chain.  This  isotopic  yield  distribution  curve  as  a  function  of  the 
distance  of  the  isotopic  charge  from  the  most  probable  charge  is  pre¬ 
sented  in  Fig.  2.  The  curve  is  seen  to  be  Gaussian  in  natiu'e  out  to 
an  absolute  value  of  (Z  -  Zp)  of  2,2,  and  exponential  in  form  from 
there  out  to  an  absolute  value  of  (Z  -  Zp)  of  4,0, 

Therefore,  the  problem  then  becomes  one  of  calculating  the  value 
of  Zp  for  each  of  the  four  types  of  fission  under  investigation. 
Ilnpirical  values  for  Zp  have  been  calculated  for  U-235  thermal  neutron 
fission  ’xsing  Wahl’s  approach.  Using  these  U-235  thermal  neutron 
fission  values  of  Zp  as  reference  values,  a  means  by  which  the  Zp 
values  for  any  other  type  of  fission  could  be  calculated  was  desired 
A  term  A  Zp  was  defined: 

A  2p  =  2p  (fission  type  of  interest) 


■Zp  (U-235  thermal  neutron  fission) 


(2) 


which  could  be  calculated  by  two  equations: 
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Isotopic  Yield  Distribution  Curve 
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AZp  =  l/2(Zc-92)  -  0.19(Ac-236) 

+  0.19  (^-2. 5)  (3) 

and,AZp  =  1/2(Z^-92)  -  0.19(A^-236)  +  0.19(0. 12 

+C^-2.5)!  where  (4) 

Zg  5  the  nuclear  charge  of  the  compotuid  nucleus 

Ag  s  the  mss  of  the  compound  nucleus  before  the  emission 
of  prompt  neutrons 

=  the  number  of  neutrons  emitted  when  the  compound  nucleus 
breaks  into  fission  fragments 
=  0.12S  + 

S*'  s  excitation  energy  of  the  fissioning  nucleus  after  bombard¬ 
ment 

E  =  the  energy  of  the  pirojectile  causing  fission 

and  C2  =  constants  being  defined  for  the  various  fission 
processes  as: 


Nucleus 

Cx 

Cz 

Pu-239 

2.9 

6.38 

U-235 

2.5 

6.5 

U-238 

2.5 

6.5 

By  use  of  Eqs  (3)  and  (4),  the  values  forAZp  for  the  fission  types 
of  interest  were  determined,  except  for  the  case  of  thermonuclear 
fission.  Previously,  (Ref  17),  this  case  was  handled  by  considering 
the  contribution  to  the  yield  from  each  of  the  two  isotopic  forms  of 
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uraniiJEi  v/hich  make  up  the  thermonuclear  device.  The  inclusion  of 
this  data  caused  the  report  to  be  classified.  In  order  to  allow  the 
widest  dissemination  of  this  report,  it  was  decided  to  calculate  the 
value  of  AZp  for  the  case  of  thermonuclear  fission  considering  an 
average  fissioning  nucleus. 

Besides  its  disadvantage  of  being  classified,  Stovall's  method 
is  not  rigorously  correct  in  using  the  experimental  data  (Ref:6), 

It  assumes  in  calculating  the  U-238  therraonuclear/i  Zp  that  the  value 
of  ^  is  3»0  ^nd  likewise  for  the  U-235  thermonuclear A  Zp,  In  both 
of  these  calculations  the  choice  of  =  3.0  is  taken  from  the 
symmetry  in  the  mass  yield  curve  of  Crocker,  Fig,  1,  ffowever,  this 
value  of  is  only  expressly  valid  for,  and  was  calculated  from,  the 
gross  assembly  of  fission  products  frcan  the  thermonuclear  fission  of 
a  mixture  of  U-238  and  U-235.  Therefore,  xmless  the  value  of  is 

independently  equal  to  3.0  for  both  the  cases  of  U-235  and  U-238 
thermonuclear  fission,  then  both  of  the  calculated  values  of  A  Zp  are 
incorrect. 

The  method  which  was  developed  to  determine  a  single  value  for 
AZp  for  the  thermonuclear  case  rested  on  the  fact  that  the  sum  of 
the  mss  numbers  of  the  light  and  heavy  fission  fragtients  plus  the 
mass  of  the  prompt  neutrons  einitted  in  fission  equals  the  mss  of 
the  compound  nucleus.  Also,  mowiedge  that  the  thermonuclear  fissile 
material  is  a  combination  of  only  uranium  isotopes  allowed  the 
deteriuination  of  the  value  of  Z^  as  being  92.  A  small  computer  pro¬ 
gram  irvas  developed  vrhich  picked  the  mass  ji-eld  corresponding  to  an 
integer  mass  nxomber  on  one  side  of  the  mass  yield  curve  and  then 
searched  the  other  side  of  the  mass  yield  curve  for  the  complementary 
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mass  number  vdth  the  identical  yield*  In  this  way  all  of  the  equal 
yield  mass  combinations  were  found,  summed,  and  an  average  was  com¬ 
puted*  This  average  was  therefore  equal  to  the  mass  of  the  compound 
nucleus  after  emission  of  the  three  prompt  neutrons  as  reported  by 
Crocker*  This  yielded  a  value  for  A  of  239*  Using  these  values  of 
Zg,  A^,  and  7/^  in  i3q  (3)  an  unclassified  value  for  thermonuclearAi^p 
was  determined.  The  values  for  A  Zp  for  the  four  cases  of  interest 
are  presented  in  Table  I, 


TABIE  I 
Values  of  Zp 


Fissioning  Nuclide 

Neutron  ihergy  (E) 

4  Z 

P 

U-235 

Thermal  (0,025  ev) 

2.5 

0 

U-235 

1  Mev 

2.7 

0.023 

Pu-239 

Thermal  (0,025  ev) 

2.9 

0.316 

( Thermonuc lear 
Mixture) 

Thermonuclear 

3.0 

-0.475 

Isomeric  Pair  Formation 

In  some  instances  it 

is  possible  for  two 

isomeric 

states  of 

same  isotope  to  be  formed.  At  present  there  is  no  method  available 
which  v/ill  predict  the  relative  amounts  of  the  tvjo  isomers  that  are 
formed  initially  in  fission.  In  this  study,  the  arbitrary,  simpli¬ 
fying  assumption  of  a  50“50  split  is  made.  This  will  be  discussed 
further  in  the  discussion  on  possible  futiure  work  in  the  further 
refinemoit  of  the  whole  prediction  method. 

Since  the  chain  and  isotopic  yields  must  be  calculated  only 
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once,  for  each  type  of  fission,  this  data  is  computed  and  put  on 
punched  computer  cards  by  a  siuall  computer  program  and  used  as  input 
data  in  the  rfiain  pix)gram.  Simple  linear  scaling  can  be  used  to  allow 
the  use  of  any  fission  yield  weapon.  The  present  method  assumes  a 
nominal  1  MT  fission  yield. 
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III ,  Decay  Chain  Characteristics 

Decay  Chain  Solutions 

The  solution  of  the  beta  spectrum  as  a  function  of  time  after 
fission  requires  the  knowledge  of  the  amovint  of  each  isotope  of 
each  mass  chain  which  is  present  at  any  time  after  fission  and  also 
its  activity.  The  decay  chains  for  the  mass  chains  of  77  through 
159  are  presented  in  Appendix  I  (Ref  10&15)*  Once  the  amount  of 
each  isotope  initially  present  is  known  then  the  problem  becomes  one 
of  solving  the  radioactive  decay  equations.  Stovall  (Ref  17s33~36) 
used  a  numerical  method  to  solve  the  differential  equations  arising 
from  the  radioactive  decay  laws  to  determine  the  concentration  of 
each  isotope  at  the  times  of  interest.  Tliis  method  had  the  severe 
practical  drawbacks  of:  (l)  taking  a  great  deal  of  computer  time; 

(2)  taking  progressively  longer  times  to  solve  for  times  of 
interest  which  were  very  large;  and  (3)  the  answer  arrived  at  was 
only  approximate  and  the  error  accianulated  the  further  the  time  of 
interest  was  from  the  time  of  fission.  It  was  therefore  decided 
that  for  increased  flexibility  and  accuracy,  it  would  be  profitable 
to  spend  the  necessary  time  to  solve  and  program  the  exact  solutions 
to  the  differential  equations. 

Once  the  decision  was  made  to  solve  the  equations  exactly,  it 
became  necessary  to  decide  which  of  two  methods  should  be  used.  The 
first  method  considered  involved  the  construction  of  a  conglomerate 
decay  chain  which  contained  within  it  all  of  the  different  decay  chain 
combinations  which  occur  in  the  mass  chains  77  through  159*  It  then 
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would  be  possible  by  setting  certain  branching  ratios  equal  to  zero, 
to  synthesize  any  of  the  actual  decay  chains.  This  conglomerate  decay 
chain  when  once  constructed  then  would  have  to  be  solved  exactly  for 
each  of  its  members.  It  was  decided  that  this  method  was  both  too 
complicated  to  solve  in  the  time  period  alloted  and  furthermore  was 
too  inflexible  for  use  in  case  some  new  decay  chain  was  encountered 
which  was  not  contained  in  the  conglomerate  decay  chain. 

The  second  method  considered,  and  the  one  which  was  actually 
utilized,  consisted  of  reducing  all  decay  chains  to  an  arithmetic 
summation  of  straight  chain  decays.  For  example,  the  decay  path  of 
mass  chain  77  is  schematically  presented  in  Fig.  3. 


at  any  time  is  a  relatively  simple  application  of  the  radioactive 
decay  laws.  The  concentration  of  D  is  computed  by  considering  a  four 
member  straight  decay  chain  ;a.th  percent  of  the  contribution  of 
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isotopes  A,  3,  and  C  and  100  percent  contribution  from  D,  The  con¬ 
centration  of  3,  v.'hich  is  a  different  isomeric  state  of  D,  is  the 
sura  of  the  contributions  of  the  five  mei;.ber  straight  chain  of  A,  B, 

C,  D,  and  3  and  the  four  laeiaber  straight  chain  A,  3,  C,  d.  In  the 
five  nember  chain  the  contribution  from  A,  B,  and  C  are  multiplied 

by  the  oroduct  of  P-  and  P, :  and  in  tho  four  raember  chain,  the  con- 

i  4 

tribution  from  A,  3  and  C  are  multiplied  by  ^  sumraing  these  two, 

care  must  be  taken  to  see  that  the  decay  of  the  isotope  3  is  only 
considered  once.  The  concentration  of  F  is  likerd.se  a  combination  of 
three  straight  chains:  The  five  nember  chains  A,  3,  C,  D,  F  and  A, 

3,  C,  3,  F;  and  the  six  member  chain  A,  B,  C,  D,  3,  F,  3ach  member's 
contribution  is  of  course  r'/eighted  by  the  appropriate  branching  per¬ 
centages,  As  can  be  seen,  this  method  involved  the  initially  large 
task  of  detemiining  all  of  the  different  paths  which  can  lead  to  the 
for.Tation  of  each  of  tlie  isotopes  contained  in  each  of  the  33  decay 
chains,  and  the  correspondingly  tedious  task  of  puttiiig  this  infor¬ 
mation  on  punched  cards  which  can  be  accepted  by  a  digital  computer. 
The  maximum  length  straight  c.aain  which  occurs  in  the  decay 
chains  as  presented  in  Appendix  3  is  ten.  Fortunately,  the  solution 
for  a  straight  chain  decay  of  any  length  folxows  a  predictable,  though 
involved  for:.a,  w’lich  increases  greatly  in  length  for  each  additional 
member  in  the  chain.  The  equations  to  be  solved  are  of  the  form: 


(5) 


whore:  X.  =  the  decay  constant  of  the  isotope  =  ln2 


T  1/2 


1 
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Tj  =  half  life  of  the  isotopes 
The  familiar  solution  for  a  single  member  decaying  is: 


o  ->t 
N  =  N,  e 
1  1 


(6) 


where:  =  the  concentration  of  the  isotope  at  time  =  t 

=  the  concentration  of  the  isotope  at  time  =  o 
X  =  the  decay  constant  of  the  isotope  as  defined  above 
t  =  the  time  of  interest 


The  solution  for  the  second  member  of  a  two  member  chain  is: 


N  =  A  C 


o  “  XT^t  “  ^  2^  .  o  ot 

°  'e—± _  +  e _  )  +  n"  e  2 


(7) 


(VV 


The  solution  for  the  third  member  of  a  three  member  chain  is: 


o  ”^1*' 

N  =  X  AN  (  ® 


+  _ §_ 


-■^2^ 


+  _ e_ 


— 


3  3=3  (X.-VtyV  (V3,){X3-y 


o  ,  -X2t 

+  X.N.  (-S 


— Xot 

a 


=  =  (X3-X3)  (X3 .3 


—Xot 


o  -Aol 
- )  ^  N  e  ^ 

-aJ  3 


(8) 
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As  can  be  seen,  the  length  of  the  solution  increases  very  rapidly. 
The  solution  for  the  tenth  meraber  of  a  ten  member  chain  contains  ten 
major  terms,  the  first  term  of  which  consists  of  ten  members  each 
with  a  denominator  nade  up  of  the  product  of  nine  differences,  the 
second  term  consists  of  nine  members  each  with  a  denominator  made  up 
of  the  product  of  eight  differences,  and  so  forth.  All  of  these 
solutions  were  written  out,  programmed  for  the  computer,  and  the 
computer  program  cards  prepared.  Since  the  types  of  sources  to  be 
considered  by  this  code  are  nuclear  and  thermonuclear  eaqjlosions, 
the  equations  were  solved  assuming  an  instantaneous,  non-time 
dependent  source  term.  This  implied  that  the  values  of  N?  in  the 
equations  were  those  concentrations  of  the  isotopes  calculated  as 
being  initially  formed  in  fission,  as  described  in  the  previous 
chapter. 

Estimation  of  lAil-aiown  Ifalf -Lives 

The  method  used  to  predict  half-lives  for  those  isotopes  for 
which  no  measured  value  exists  is  the  same  one  used  previously, 

(Ref  17:23-25), 
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IV,  Theory  of  Beta  Decay 

Description  of  Beta  Decay 

In  the  process  of  beta  decay,  the  nucleus  undergoing-  decay  eriits 
a  beta  particle  and  an  antineutrino.  This  causes  the  resultant 
nucleus  to  increase  its  charge  by  one  unit.  In  tliis  decay  process, 
the  mss  number  of  the  original  and  final  isotope  is  considered  the 
same.  This  process  is  presented  in  equation  form  below: 


gA 

Z+1  +  +  V  ° 

-1  o 


(9) 


Since  the  decay  process  involves  a  three  body  problem  v^ith  the  final 
nucleus,  beta  particle,  and  antineutrinoj  it  is  not  possible  to  pre¬ 
dict  a  nonoenergetic  value  for  the  energy  of  the  beta  pcirticle. 

Rather,  since  the  atomic  mass  difference  beWeen  the  initial  nucleus 

(  A  )  and  the  mass  of  the  final  nucleus  )  represents  the  total 

Z  ^ ^ 

kinetic  energy  available  for  partition  between  the  beta  particle 
(  and  the  antineutrino  (Jo  °),  the  process  of  beta  decay  can  be 

characterized  by  the  instance  in  -which  practically  all  of  this  kinetic 
energy  is  possessed  by  the  beta  particle  and  essentially  none  by  the 
antineutrino.  This  maximum  kinetic  energy  is  commonly  referred  to  as 
the  beta  end  point  energy  (jj],P,E),  The  observed  beta  particle 
energies  therefore  form  a  spectrum  starting  at  zero  energy  and  ex¬ 
tending  to  a  maximum  value,  the  beta  end  point  ener,gy.  The  shape  of 
this  beta  spectrum  has  been  studied  and  a  theoi-y  by  which  it  can  be 
calculated  and  predicted  has  been  developed  by  Fenrd  (Ref  7i548-62), 
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The  mathecatical  derivation  of  the  Fermi  theory  vjas  discussed  in 
Stovall’s  report  and  vdll  not  be  reproduced  here.  The  spectral  shapes 
predicted  by  this  theory  correspond  to  vdiat  is  known  as  "allowed" 
transitions.  There  are  other  spectral  shapes  observed  experimentally 
which  are  known  as  "forbidden"  transitions.  The  terns  "allowed"  or 
"forbidden"  transitions  do  not  indicate  vrtiether  or  not  these  modes  of 
decay  are  any  more  or  less  likely  to  occur  in  naturej  they  simply  refer 
to  those  transitions  whose  spectral  shapes  are  directly  predicted  by 
the  Fermi  theory  and  those  whose  spectral  shapes  are  modified  ver¬ 
sions  of  the  spectral  shapes  predicted  by  the  Fermi  theory. 

Shape  Correction  Factors 

In  order  to  produce  the  actual  spectral  shapes  for  the  forbidden 
transitions,  a  set  of  correction  factors  are  used  (Ref  ll:557-S)j 
these  correction  factors  are  commonly  known  as  spectral  shape  factors. 
The  determination  of  the  degree  of  forbiddeness  of  any  specific  beta 
decay  and  therefore  the  deterndnation  of  the  appropriate  correction 
factor  to  be  applied  to  the  results  of  the  Fermi  theory  is  accom¬ 
plished  by  the  calculation  of  a  "log  ft"  value.  The  mathematical 
expression  for  the  calculation  of  this  "log  ft"  value  (Ref  17:64)  is 
a  function  of  the  atomic  number  (Z),  of  the  nucleus  undergoing  decay, 
the  maximum  beta  energy  of  the  decay  (E,P,E. ),  and  the  lialf-life 

(Ti)  of  the  decay, 

2 

Table  II  lists  the  log  ft  values  and  the  appropriate  shape 

factors  (C^^),  In  the  shape  factors,  the  following  definitions  apply: 

£  =  E/m^c*^  where  E  =  total  energy  of  the  beta  (kinetic  +  rest  mass), 

2 

m^  =  rest  mass  of  the  beta,  c  =  speed  of  light;  and  ^  ^ taax/^o^ 

where  E  max  =  T  +  T  ,  T _ _  =  maximum  kinetic  enerar  of  the  beta 

max  o^  max  ^ 
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''FORBIDDMSSS'' 

LOG  ft  1/2 

10 

1 

1 

SHAPE  FACTOR 

(c  ) 

ft 

SUPxJfiA  LL0.1/ III) 

3.5jp.5 

C  ~  Ifriity 
ft 

1 

ALL0;»'2D 

5.0j:i.o 

1 

1 

C  ~  Unity 
ft 

1st  FORBIDDEN 

7. 0+1,0 

- 1 

O  2  2 

Cf-  P  +1 

UI.'IQUE 

1st  FORBIDDEN 

8.5jp.5 

C  'N-'  l/l2(p^+q^) 
ft 

2nd  FORBIDDi:; 

13.0 

pVlO/3p^q^+q^ 

xL 

1 

3rd  FOIiBIDDBN 

~  13.0  ; 

p° 

+7p^q^(p^-«5^)  T 

Trarisi“ornation  Relations  j 

q2.c  (f.-e)2 


Table  II 

iinpirical  Log^^ft  Ranges  and  Shape  Factors 
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(j;.P,E. ),  and  =  rest  mass  energy  of  the  beta. 

Therefore,  in  the  calculation  of  the  beta  spectrum  due  to  any 
single  isotope  in  any  of  the  mass  chains,  the  follovring  steps  are 
taken.  Consider  the  hypothetical  decay  scheme  in  Fig,  4, 


Consider  first  the  beta  decay  characterized  by  an  E.P.E,  of  1,0  I-fev. 
The  first  thing  which  is  done  is  to  compute  the  spectral  shape  by  the 
Fermi  theory.  Then  the  log  ft,  value  is  computed,  the  appropriate 
shape  factor  determined,  and  the  spectral  shape  is  then  modified  by 
this  shape  factor.  The  same  procedure  is  followed  for  the  beta  decays 
with  E,P,E,  's  of  2,0  and  3,0  llev,  IVhen  this  is  done,  then  the 
spectral  shapes  for  these  three  betas  are  weighted  according  to  the 
percent  of  the  decays  going  by  each  route  and  summed.  This  sum  is 
then  normalized  to  give  an  area  of  one  under  the  conglomerate  spectral 
curve.  This  result  is  then  referred  to  as  the  normalized  spectral 
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shape  for  the  deca^r  of  isotope  A  to  isotope  B, 

Treatment  of  Isotopes  v-dth  Unknovm  Decay  Schcr.ies 

In  the  case  in  which  on.lj'’  the  identity  of  the  initial  and  final 
isotopes  are  knovm,  the  respective  half-life  values  are  either  known 
or  estimated,  and  the  scheme  of  beta  decays  by  vdiich  isotope  A  decays 
to  isotope  B  is  not  known,  an  approxrLmation  technique  is  eiiiployod,  A 
quantity  o( , 


(X  = 

Srp 


Aliere  E  =  The  total  beta-antineutrino  energy  (E.P.E. ) 

P— 

=  The  total  disintegration  energy 

as  defined  above,  is  used  in  ;diich  the  value  of  E^  is  taken  to  be  the 
difference  in  the  atomic  masses  of  isotopes  A  and  B,  Stovall 
(Ref  17:54)  used  two  different  empirical  nass  formulas,  the  Ouaeron 
(Ref  2:1021-32)  and  the  Coryell  (Ref  3:305-334)#  No  attempt  was  made 
to  show  that  one  or  the  other  of  these  nass  formulas  was  superior  to 
the  other,  A  method  to  do  this  is  presented  in  the  section  on  the 
discussion  of  possible  future  work  in  the  further  refinement  of  the 
whole  prediction  technique.  The  Coryell  mass  formula  was  employed  in 
the  actual  calculations  since  it  yielded  solutions  which  agreed  well 
with  the  available  experimental  data. 

For  unimown  decays,  therefore,  the  value  of  the  mass  difference 
betv/een  the  parent  and  daughter  isotopes  in  the  decay  was  multiplied 
by  the  value  of  to  determine  the  hypothetical  beta  E,P.E,  By  this 
method,  the  parent  isotope  is  assumed  to  decay  by  one  beta  to  an 
excited  state  of  the  daughter  isotope.  By  assuming  one  value  of 
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for  all  isotopes  vdth  unimovm  decay  scheraes,  it  v/as  illustrated  that 
a  value  of  =  2/3  aatched  the  observed  experiiaeiital  data  fairly  well. 
Tfiis  was  the  value  of  used  in  the  coraputer  programs. 

It  was  decided  that  one  of  the  jrain  objectives  of  this  research 
would  be  the  discovery  of  an  alternate  technique  by  v/hich  the  value 
of  <X  could  be  detci’iidned,  Tlie  approaches  investigated  did  not  lead 
to  any  usable  teclinique,  although  one  showed  sufficient  promise  of 
being  of  value  to  warrant  further  work.  The  approaches  investigated 
are  described  below. 

Alternate  Approaches  to  the  Determination  of 

The  first  approach  was  to  try  to  determine  the  variation  of 
directly  as  a  function  of  half  life.  This  was  accomplished  by  con¬ 
sidering  all  of  the  known  decays  contained  in  the  d3  decay  chains  of 
interest.  For  each  of  the  known  decays  a  value  of  o(  was  determined. 

For  example,  consider  the  decay  scheme  illustrated  in  Fig,  4»  In 
this  case  cX  would  be  determined  as  follows: 

o<  =  0.0).  (3^.). 

(4.0)  4 


=  0,4375  (  24) 

A  total  of  sixty-two  points  were  calculated  and  plotted  on  a  graph  of 
log  Ti  vs,  ,  There  was  no  observable  correlation.  It  was  therefore 
concluded  that  there  was  no  simple  relationship  between  the  half-life 
and  the  value  of  c<  , 

The  second  approach  was  to  attempt  to  find  a  correlation  between 
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the  calculated  values  of  and  the  quantities  and  |Z-Zpj  for 

the  known  decay  schemes  and  assume  that  they  would  be  valid  for  the 
unknown  decays.  Plots  were  made  in  which  all  of  the  data  was  pre¬ 
sented  on  a  single  sheet  and  also  in  which  the  data  were  presented  on 
four  sheets  according  to  whether  the  decaying  nuclide  had  values  of 
Z  (number  of  protons)  and  N  (number  of  neutrons)  which  were  odd-odd, 
odd-even,  even-odd,  or  even-even,  »Vhen  the  points  were  plotted,  the 
distribution  was  random. 

The  third  method  used  to  try  to  determine  more  rigorously  the 
value  of  ck  was  suggested  by  a  report  prepared  by  James  Griffin,  of 
the  Los  Alamos  Scientific  laboratory,  of  the  decayed  gammas  from 
fission  fragments  (Ref  11),  Although  his  method  for  theoretically 
synthesizing  the  fission  process  was  not  considered  to  be  as  rigorous 
as  the  method  described  in  this  paper  and  in  Stovall’s,  the  report 
contained  plots  of  the  total  decay  energy  per  fission  per  second  and 
the  gamma  energy  per  fission  per  second  from  a  time  of  approximately 
1  X  10“^  seconds  to  45  seconds  after  fission  for  the  cases  of  several 
different  fissioning  nuclides.  These  theoretical  plots  compared  well 
with  the  experimental  points  he  used. 

By  computing  the  difference  betv/een  the  total  decay  energy  and 
the  gamma  energy  at  the  sane  times  it  was  possible  to  determine  the 
beta  energy  per  fission  per  second  at  any  time  of  interest  in  the  time 
region  covered  by  the  original  graphs.  The  desired  result  of  this 
approach  was  to  be  the  determination  of  the  value  of  cX  as  a  function 
not  of  half-life  but  rather  to  determine  values  of  o<  for  a  finite 
number  of  half-life  groups.  To  do  this  it  was  necessary  to  divide  the 
isotopes  into  six  groups  according  to  their  half-life,  as  shown  in 
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Table  III, 


TA3IE  III 


Half-Life 

Grouping 

Gp 

1 

1 

oo 

>  T 

1/2 

> 

67 

sec 

Gp 

2 

67 

>  T 

1/2 

2 

21 

sec 

Gp 

3 

^3 

21 

>  T 

1/2 

2 

8 

sec 

Gp 

4 

^4 

8 

>  T 

1/2 

4 

sec 

Gp 

5 

^5 

4 

>  T 

1/2 

> 

2 

sec 

Gp 

6 

“^6 

2 

>  T 

1/2 

> 

0 

By  sa-xiing  the  products  of  the  activity  (A^),  the  end  point 
energy  (E.P.ij. ),  and  the  appropriate  value  of  for  all  isotopes 
v/ith  unknown  decay  schemes  present  at  several  times  of  interest,  and 
subtracting  the  contribution  from  those  isotopes  whose  decay  schemes 
are  known,  it  should  be  possible  to  create  an  overdetermined  set  of 
equations  vd.th  the  values  of  the  « 's  as  the  unloiowns.  From  the  graphs 
contained  in  the  Griffin  reports  (Ref  U)  the  values  of  the  beta 
energy,  for  various  times,  was  computed  and  a  small  computer  program 
was  vn*itten  to  solve  the  overdetermined  set  of  equations, 

A  set  of  thirteen  eqiiations,  (Eqs,  12  through  24),  was  generated 
coinciding  with  thirteen  values  of  the  total  beta  energy  calculated 
from  Griffin's  plots,  j^ch  equation  represented  the  sitviation  at  one 
particular  tiiae.  The  times  chosen  were  from  1  to  45  seconds  after  the 
instantaneous  fission  event.  The  set  of  equations  is  presented  below. 
The  values  were  calculated  for  one  l-xH  of  fission  and  a  common  factor 
of  1  X  10“^^  was  applied  to  each  term  in  order  to  reduce  the  size  of 


the  numbers. 
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(sec 

) 

0^ 

o( 

cX 

cx 

c 

1 

2 

4. 

6 

2 

!• : 

+  7.04 

3.28 

4- 

5.46 

4- 

26.1 

4- 

72.5 

4- 

54.4 

= 

136. 

(12) 

1.5: 

3.45 

+  6,83 

3.43 

+ 

5.66 

4“ 

25.7 

+ 

64.6 

4- 

42.5 

= 

103. 

(13) 

2,: 

3.46 

+  6,60 

■f 

3.56 

4- 

5.31 

4- 

25.1 

4- 

57.4 

4- 

33.4 

38.2 

(14) 

3.: 

3.48 

+  6,11 

-f 

3.73 

4- 

6.00 

4- 

23.6 

+ 

44.8 

4- 

20,8 

= 

67.3 

(15) 

4.: 

3.47 

+  5.61 

•f 

3.35 

4- 

6.08 

4- 

21.8 

4- 

34.7 

4- 

13.1 

= 

47.1 

(16) 

6, : 

3.40 

+  4.65 

+ 

3.94 

4- 

5.99 

-f 

18.1 

4- 

20.6 

+ 

5.29 

30.7 

(17) 

8,: 

3.31 

+  3.79 

+ 

3.92 

4- 

5.71 

4- 

14.6 

4- 

12.1 

4- 

2.16 

24.4 

(18) 

10,: 

3.21 

+  3.07 

+ 

3.85 

4- 

5.34 

4- 

U.7 

4- 

7.12 

4* 

.891 

= 

18.7 

(19) 

15.: 

2,96 

+  1.81 

+ 

3.57 

4- 

4.33 

4- 

6.44 

4- 

1.95 

4- 

.100 

= 

9.86 

(20) 

20,: 

2,76 

+  1,08 

+ 

3.26 

4- 

3.43 

4- 

3.51 

+ 

.564 

4- 

.012 

= 

5.57 

(21) 

30,; 

2,46 

+  .456 

-f- 

2.68 

4- 

2.11 

+ 

1.06 

+ 

.061 

+^^.00 

= 

2.80 

(22) 

40,: 

2,25 

+  .257 

+ 

2.21 

4- 

1.30 

4“ 

.330 

4- 

.011 

+  A/,00 

1.90 

(23) 

45.: 

2,17 

+  .213 

2.01 

4- 

1.03 

4- 

.187 

4- 

.007 

4- 

'\j  *00 

1.52 

(24) 

The  term  refers  to  the  contribution  resulting  from  the  decays 
of  the  isotopes  with  known  decay  schemes.  This  term  is  computed  as: 


C 


1 


N 

i=l 


Con.  X  X. 
i  1 


M 

j=l 


Frac  X  E.P.E. 
ij  ij 


(25) 


Where  Con.  = 
1 


1 

Frac. 


ij 


E. P .  E» ,  , 

IJ 

N 


Concentration  of  isotope  i  at  the  time  of  interest 
Decay  constant  of  isotope  i 

=  The  fraction  of  isotope  i  which  decays  by  a 
specific  end  point  energy, 

=  A  beta  end  point  energy  by  which  isotope  i  decays 
=  The  number  of  isotopes  with  knov/n  decay  schemes 
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H  =  The  number  of  different  betas  by  v/hich  the  isotope 
decays. 

The  term  Cp  is  the  constant  derived  from  Griffin’s  Jieport 
(Refill),  and  it  is  calculated  as; 


Cp  —  O.m,  —  G,ii>,  (26) 

Vihere 

D,E,  ==  Total  decay  energy  from  fission 
G,E,  =  Total  gamma  energy  from  fission 
The  terms  in  columns  through  are  computed  as; 

P 

^  Con  X  M,D.  .  X  X.  .  (27) 

J  i  =  1 

Vihere  o{^  =  The  alpha  value  in  half  life  group  j 

v) 

Con. .  =  Concentration  of  isotope  i  in  half  life  group  j  which 
decays  in  an  unlmown  manner, 

=  Decay  constant  of  isotope  i 

M,D,^^  =  The  mass  difference  energy  between  isotope  i  and  its 
decay  product 


In  order  to  have  a  meaningful  solution  to  the  equations  viiich 
result  from  condensing  the  two  constants,  in  other  words  one  in  which 
the  values  of  alpha  range  betv/een  zero  and  one,  it  is  necessary  that 
be  smaller  than  C^,  This  is  the  case  in  Eqs  (12)  through  (22), 
however,  in  the  last  two  equations  the  opposite  is  true.  This 
occurrence  forces  the  solutions  to  contain  negative  values  for  some 
of  the  alphas. 
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There  are  several  possibilities  waicii  could  account  Tor  this 
apparent  discrepancy.  The  first  is  the  data  used  to  calculate  Cp, 

The  energy  data  is  reinforced  by  the  experiiaental  points  ii>dth 

which  it  is  compared;  however,  the  total  decay  ener^  curve  is  not 
reinforced  by  experimental  data.  Upon  trying  to  check  on  the  reference 
given  by  Griffin  for  this  curve  it  was  discovered  that  the  original 
work  was  presented  in  Gernian  and  no  translation  was  available  in  the 
time  remaining  for  the  research. 

The  second  possible  explanation  for  the  discrepancy  is  based  on 
the  half-life  prediction  method.  If  this  prediction  method,  vdiich 
was  originally  done  by  Holies  and  Ballou  (Ref  l),  predicts  half-lives 
which  are  too  short,  it  would  cause  the  concentration  of  those 
isotopes  which  are  included  in  to  be  higher  than  it  should  be.  By 
increasing  this  concentration,  the  value  of  as  calculated  by  i^(25) 
would  increase  and  overwhelm  the  value  of  Cp, 

The  lack  of  computer  time  available  and  the  lack  of  time  re^iain- 
ing  in  the  research  period  precluded  doing  any  further  work  in  tiiis 
area.  In  view  of  the  difficulties  encountered  in  this  effort,  it  i.ay 
be  that  the  solution  to  this  problem  may  have  to  await  further  experi¬ 
mental  evidence. 
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V,  Spectrmn  Calculations 


Xethod  of  Solution 

As  explained  in  Chapter  II,  the  initial  concentration  of  the 
fission  products  from  any  of  the  types  of  fission  considered  is  com¬ 
puted  and  put  on  punched  cards.  Then  the  normalized  spectral  shapes 
for  all  of  the  isotopic  decays  contained  in  all  of  the  eighty-three 
decay  chains  are  computed  and  stored  on  mgnetic  tape.  These  two 
data  records  are  then  combined  by  a  program  which  puts  this  data  on 
another  magnetic  tape  which  compiles  the  explicit  data  needed  to 
solve  problems  involving  any  of  the  four  types  of  fission  separately. 
Therefore,  in  order  to  solve  a  problem  it  is  only  necessary  to  use 
this  data  tape  as  input  into  the  program,  described  in  Chapter  III, 
which  contains  the  exact  solutions  to  the  differential  equations 
resulting  from  the  radioactive  decay  law.  By  these  means  it  is  no 
longer  necessary  to  compute  the  isotopic  yields  and  the  spectral 
shapes  each  and  every  time  a  problem  is  to  be  solved  as  was  necessary 
using  Stovall's  programming  arrangement  (Ref  17),  This  reduces  the 
running  time  necessary  on  the  computer  from  approximately  an  hour  for 
a  typical  problem  to  twenty  minutes  to  prepare  the  composite  data 
tape  one  time  only  and  thereafter  only  six  minutes  per  problem.  This 
time  advantage  becomes  much  greater  if  the  time  after  burst  at  which 
the  spectrum  is  desired  is  much  longer  than  five  or  ten  minutes. 

This  is  true  because  in  the  previous  method  using  the  NDA  subroutine 
(Ref  17:34-36)  the  iterative  solution  had  to  start  at  zero  time  and 
proceed  to  the  tine  of  interest  by  siaall  incremental  steps  while  the 
present  program  uses  the  exact  solutions  to  the  differential 
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equations  vd.th  the  tine  of  interest  used  as  a  sinple  variable. 


A  schematic  presentation  of  the  data  used  and  the  means  by  vdiich 
the  spectrum  is  calculated  is  presented  in  Pd-g.  5» 


lass  Yield  Curve  +  Glendenin  S, D, C»  Theory  »  Isotopic  Yields 

Decay  Chain  Data  +  Ferini  Theory  +  Shape  Noriialized  Spec- 

Factors  +  Empirical  lass  Formulas  •  tral  Shapes 

Isotopic  Yields  +  Normalized  Spectral  Composite  Data 

Shapes  for  all  I^pes  of  Fission  - -  Tape 

Composite  Data  Tape  +  Decay  Equations  _ _  Total  Fission 

+  Parameters  of  Specific  Problems  Spectrum  of 

Specific  Problems 

Fig.  5 

A  Sumraaiy  of  the  Beta  Spectrum  Calculations 
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VI .  Results  and  Conclusions 

Results 

Predictions  of  the  time  dependant  beta  spectrum  at  various  times 
for  the  four  cases  of  fission  are  presented  in  Figs.  6,  7»  and  9. 

At  the  present  time,  predictions  of  the  spectrum  resulting  from  high 
altitude  detonations  are  based  on  the  steady— state  fission  spectrum 
of  Uranium  235  (Ref  20).  The  shape  of  this  spectrum  is  shown  on  Fig 
6.  At  the  early  times,  it  represents  fairly  well  the  high  energy 
portion  of  the  spectrum,  but  grossly  overestimates  the  low  energy 
portion.  At  late  times  it  does  a  fairly  accurate  job  of  representing 
the  low  and  medium  energy  portions  of  the  spectrum  but  leads  to  an 
overestimation  of  the  high  energy  portion.  Therefore,  depending  on 
the  time  after  fission  when  the  beta  particles  become  trapped  by  the 
magnetic  field,  the  predicted  spectrum  will  be  in  error  at  either  the 
low  or  high  end. 

In  order  for  any  prediction  to  be  useful,  it  is  most  desirable 
that  its  result  should  compare  favorably  with  an  independent  source. 

In  this  case  the  independent  source  is  the  experimentally  determined 
results  of  West  (Ref  19).  The  grapliic  comparisons  are  presented  in 
Figs.  10  through  l6,  and  are  discussed  below. 

Conclusions 

In  comparing  the  predictions  of  the  model  with  the  experimental 
results,  it  is  very  important  to  consider  exactly  how  the  two  results 
were  obtained.  The  method  by  which  the  experimental  data  was  obtained 
is  very  poorly  described.  In  a  telephone  conversation  with  West  made 
in  the  Fall  of  1964  by  Stovall,  it  was  discovered  that  the  data 
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Fig.  10 

Comparison  with  West's  data  at  t  =  1.04  sec 
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Fig.  12 

Comparison  with  West's  data  at  t  =  18,9  sec 


kO 


No,  Beta/Siergy 


G;^VP‘Vs/65-3 


Fig,  13 

Cnmparison  with  VJest’s  data  at  t  =  39,0  sec 
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reported  for  1.04  seconds  was  obtained  by  irradiating  the  sample  for 
1/2  second,  waiting  for  l/2  second  and  then  counting  over  a  l/2  second 
period.  Likewise  the  second  set  of  data  was  obtained  with  a  one 
second  irradiation,  one  second  wait,  and  a  one  second  counting  period. 
Tlie  exact  method  by  wliich  the  remaining  sets  of  data  were  obtained 
was  not  made  clear  during  the  conversation.  In  order  to  synthesize 
this  experimental  approach  with  th.e  prediction  method,  it  was 
necessary  to  make  some  assumptions.  Since  the  prediction  method 
assumes  instantaneous  formation,  it  was  assumed  that  this  occurred  at 
the  middle  of  the  exposure  time,  A  series  of  instantaneous  spectra 
were  calculated  at  times  ranging  from  3/4  to  1  I/4  seconds  after  for¬ 
mation,  These  instantaneous  spectra  were  integrated  using  the  trape¬ 
zoidal  method.  The  results  are  presented  in  Figure  10,  The  areas 
under  the  predicted  points  and  under  West's  points  were  computed  using 
the  trapezoidal  rule  and  West's  points  were  adjusted  to  yield  the 
same  area.  The  factor  by  wtiich  West's  points  had  to  be  divided  was 
designated  as  C,  In  the  first  instance  the  value  of  C  was  2,99.  The 
same  procedure  was  followed  for  the  second  data  set  yielding  a  value 

of  C  of  1.42. 

When  the  exposure  or  measuring  period  is  small  compared  to  the 
waiting  period,  the  data  sets  should  agree  well  with  the  instantaneous 
spectrum.  This  hypothesis  was  tested  by  plotting  the  instantaneous 
spectra  at  times  used  by  West  against  the  experimental  points.  Once 
again  the  value  of  C  was  computed.  For  the  times  of  298>  147»  and  78 
seconds  the  values  of  C  ranged  between  1.30  and  1,51,  a  consistent 
difference  between  theoretical  and  experimental  results.  At  39 
seconds  the  value  of  C  =  1,66  which  tended  to  suggest  that  the 
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hypothesis  was  not  valid  at  this  time.  The  comparison  at  18,9  seconds 
supported  this  belief  since  the  value  of  C  =  2.76. 

It  is  encouraging  to  note  the  close  agreement  between  the  pre¬ 
dicted  spectra  and  the  adjusted  experimental  points.  In  order  to 
make  a  more  thorough  analysis  of  the  comparisons  it  is  imperative 
that  the  exact  method  by  which  the  data  points  were  obtained  be 
determined.  Unfort\inately  West  never  published  a  full  report  of  his 
work,  only  an  abstract.  Further  suggested  areas  of  work  which  could 
improve  the  making  of  comparisons  are  presented  in  the  section  on 
recommended  areas  for  further  effort. 

Applications  of  the  Results 

In  addition  to  the  academic  knowledge  of  the  beta  spectrum  as  a 
function  of  time  after  fission,  the  predictor  method  has  an  important 
practical  usefulness  as  a  research  tool. 

In  recent  years  much  has  been  learned  about  the  phenomenon  of 
trapped  radiation  in  the  area  of  outer  space  near  the  earth.  The  area 
in  which  this  trapped  radiation  plays  an  important  role  is  also  an 
area  in  which  many  scientific  and  military  manned  and  unmanned 
satellites  are  now  traversing  and  are  programmed  to  be  using  in  the 
future,  V/hile  the  beta  particles,  which  may  be  trapped  as  a  result 
of  a  hi^  altitude  or  near  earth  space  nuclear  or  thermonuclear  deto¬ 
nation,  will  not  themselves  constitute  an  appreciable  direct  hazard 
to  satellites  traversing  this  region,  their  interaction  with  the 
structural  materials  of  the  shell  of  the  satellite  can  result  in  an 
appreciable  X-ray  dose  as  a  result  of  bremsstrahlung  reactions. 

The  predictor  model  described  in  this  report  can  be  used,  with 
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minor  changes  in  the  output  section,  to  supply  data  vrtiich  can  be  used 
to  further  other  research  studies  in  the  field  of  nuclear  weapon 
effects.  One  problem  to  which  this  method  can  address  itself  concerns 
the  study  of  fission  fallout  following  a  nuclear  or  thermonuclear 
detonation.  In  the  study  of  this  problem  one  of  the  basic  items  of 
information  necessaiy  is  the  composition  of  the  nuclear  debris  at  the 
time  that  condensation  begins.  By  use  of  the  predictor  program,  the 
concentrations  of  each  of  the  isotopes  formed  directly  in  fission  and 
their  subsequent  decay  products  can  be  determined.  By  a  knowledge  of 
the  yield  of  the  device,  the  temperature  at  which  different  elements 
solidify,  and  the  time  history  of  the  temperature  of  the  fireball,  it 
should  be  possible  to  determine  the  concentration  of  each  of  the 
radioisotopes  at  the  time  of  solidification  and  therefore  to  determine 
the  composition  of  the  fallout  at  its  start. 

Recommended  Areas  for  Further  Effort 

Although  the  predictor  model  as  presented  in  this  report  gives 
the  beta  spectrum  as  a  function  of  time,  there  are  still  many  areas 
in  which  it  can  be  substantially  improved  in  terms  of  accuracy  and 
flexibility. 

Tne  first  area  is  the  continuing  effort  to  update  the  input  data. 
As  more  and  more  decay  schemes  and  half-lives  are  inserted  into  the 
program  to  replace  unknown  decay  schemes  and  predicted  half-lives, 
less  and  less  importance  will  have  to  be  placed  on  the  empirical  mass 
formulas  and  on  the  assumed  value  of  In  conjunction  with  this  is 
the  incorporation  of  new  data  which  may  show  that  certain  isotopes 
do  not  appear  in  the  decay  schemes  shown  in  Appendix  A,  It  would  also 
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be  of  :’reat  assistance  if  some  theoretical  method  could  be  developed 
to  predict  the  fraction  of  initial  yield  which  is  accounted  for  by 
each  member  of  the  isomeric  pairs.  It  was  shown  by  Porile  (Ref  16) 
that  in  the  fo miration  of  such  isomeric  pairs  the  member  with  the 
higher  spin  state  is  usually  formed  in  a  greater  amount.  The  deter¬ 
mination  of  the  ratio  of  formation,  hov/ever,  has  not  been  refined  to 
such  a  state,  and  complete  experimental  evidence  has  not  been 
sufficiently  obtained  to  justify  discarding  the  50-50  hypothesis. 

Some  recent  papers  (Ref  5  and  Ref  13)  also  indicate  that  the  50-50 
split  used  in  this  paper  is  incorrect. 

The  second  area  concerns  the  choosing  of  a  single  empirical  mass 
formula  for  use  with  those  very  short  half-life  isotopes  for  which 
no  decay  scheme  is  now  known,  and  for  which  no  decay  scheme  will 
probably  be  found.  A  comparison  could  be  obtained,  and  a  ch.oice  made 
between  the  Cameron  and  Coryell  formulas  in  the  following  manner. 

A  spectrum  could  be  obtained  for  a  time  of  interest  at  which  the 
activity  of  those  isotopes  with  unknown  decay  schemes  is  very  small. 
Since  practically  all  of  the  isotopes  which  fall  into  this  category 
have  short  half-lives,  this  would  mean  picking  a  time  of  interest  of 
perhaps  ten  minutes.  Next  a  similar  spectrum  could  be  obtained 
except  that  in  this  case  all  of  the  isotopes  will  be  considered  to 
have  unknown  decay  schemes  and  the  empirical  mass  difference  deter¬ 
mined  by  one  of  the  mass  formulas  should  be  multiplied  by  the  value  of 
in  order  to  determine  the  E.P.E.  and  therefore  to  determine  its 
spectrum.  This  procedure  should  be  followed  using  each  of  the  mass 
formulas  and  the  results  compared  with  the  first  spectrum  obtained 
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using  the  known  decay  schemes,  A  comparison  of  this  type  could  give 
a  sound  basis  by  wtiich  one  of  the  nic.ss  formulas  could  be  eliminated 
in  future  work  in  this  field. 

The  third  area  of  future  work  is  in  the  area  of  making  the  pre¬ 
dictor  more  flexible.  As  it  now  stands,  the  source  is  considered  to 
be  of  an  instantaneous  nature.  This  is  of  course  the  correct  approach 
to  take  when  consideiring  a  nuclear  or  thermonuclear  weapon  detonation. 
However,  since  the  validity  of  this  model  is  ultimately  based  upon  its 
comparison  with  the  available  experimental  data^  it  is  not  un¬ 
reasonable  to  attempt  to  adapt  the  model  so  that  it  duplicates,  as 
closely  as  possible,  the  conditions  under  which  the  experimental 
evidence  was  obtained.  In  the  experimental  cases  the  fissioning 
material  was  not  exposed  to  an  instantaneous  burst  of  neutrons  but 
rather  a  steady  source  of  neutrons.  This  means  that  the  fission 
occurred  at  an  approximately  constant  rate  over  a  finite  time  interval. 
In  the  comparison  which  was  made  with  West's  data  (Ref  19)  this 
exposure  time  was  one-half  of  a  second.  In  this  case  the  assumption 
of  instantaneous  formation  does  not  appear  to  be  gross;  however,  in 
some  of  the  later  measurements  this  exposure  time  became  long  enough 
so  that  an  appreciable  amount  of  decay  took  place,  thereby  making  the 
assumption  that  the  decay  started  with  the  isotopic  concentrations 
typical  of  instantaneous  formation  in  error.  To  correct  this  defi¬ 
ciency,  it  is  recommended  that  the  differential  equations  of  the 
radioactive  decay  law  be  solved  considering  a  constant  rate  production 
term.  This  will  have  the  added  advantage  of  allowing  the  predictor 
method  to  also  handle  the  problem  of  fission  product  inventory  in  a 
nuclear  reactor. 
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The  eauations  v/hich  must  be  solved  are  of  the  form: 


dN. 


i  =  -XN  +  R 


(28) 


dt  i  i  i 

The  solution  for  the  one  and  two  member  chains  are  shown  in  Eqs 
(29)  and  (30). 


(29) 


As  was  the  case  in  the  solutions  of  the  decay  equations  considering 
instantaneous  fonmation,  the  solutions  become  very  long  very  quickly. 

The  fourth  area  of  recommended  effort  concerns  the  problem  of 
computing  the  integrated  beta  spectrum.  Solution  of  this  problem 
would  also  add  to  the  flexibility  of  the  program.  As  in  the  example 
of  the  West  data,  the  experimental  evidence  available  as  to  the  beta 
spectrum  is  not  measured  instantaneously  but  ratlier  over  a  finite 
time  interval.  In  order  to  simulate  this  for  purposes  of  comparison, 
the  present  method  calculates  the  instantaneous  spectra  at  many  times 
durin,^  this  period  and  sums  them.  By  integrating  the  exact  solutions 
to  the  differential  equations,  an  admittedly  tedious  but  relatively 
simple  task,  it  would  therefore,  be  possible  to  duplicate  the  exact 
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conditions  under  which  the  experimental  evidence  was  obtained. 

After  the  recommended  improvements  discussed  above  are  incorpo¬ 
rated  into  the  predictor  model,  the  next  logical  step  would  be  to 
attempt  to  combine  this  program  with  existing  programs  which  predict 
the  extent  to  which  the  beta  particles  can  be  expected  to  be  trapped 
by  the  earth's  magnetic  field.  The  results  of  such  a  predictor  model 
could  then  be  used  in  other  programs  which  compute  the  radiation 
dosage,  dose  rate,  and  physiological  or  physical  effects  of  the 
bremsstrahlung  radiation  on  the  equipment  and/or  personnel  in  the 
satellite. 
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Appendix  A 
Chain  Yield 


Table  IV  Mass  Chain  Yields/lO^  fissions 


A 

U -23 5th 

U-235  f.sp. 

Pu-239th 

T.  N. 

77 

0.8249 

4.8235 

0.5223 

1.05 

78 

2.0098 

7.6204 

1.1745 

2.02 

79 

5.5622 

11.9882 

2.4580 

3.66 

80 

9.4156 

18.7126 

4.8I5I 

6.97 

81 

14.1037 

30.6543 

9.5527 

10.51 

82 

24. 1998 

42.7472 

16.5159 

22.42 

83 

54.2830 

77.2805 

29.0089 

50.12 

84 

100.5002 

91.5459 

46.7012 

78.21 

85 

131.0524 

121.9460 

54.0153 

110.3 

86 

201.0324 

172.7850 

75.8235 

145.3 

87 

249.0276 

250.7498 

90.8349 

186.5 

88 

356.8627 

343.5819 

140.3722 

226.7 

89 

471.7097 

421.6817 

171.1613 

268.0 

90 

572.0288 

507.3745 

223.6839 

308.8 

91 

582.5530 

552.4071 

289.3741 

366.8 

92 

597.6441 

581.9531 

308.6625 

411.0 

93 

643.5574 

609.3539 

395.6414 

470.9 

94 

638.5754 

619. 5415 

445.8155 

504.0 

95 

625.7403 

627.0754 

499.3532 

510.4 

96 

625.5302 

635.8294 

514.2104 

521.7 

97 

608.0432 

614.4070 

558.6327 

532.4 

98 

578.6314 

613.3496 

581. 2a9 

540.2 
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Table  IV  ( con' t) 

A 

U -23 5th 

U-235  f.sp. 

Pu-239th 

T.  N. 

99 

606.6437 

619.9075 

609.0357 

547.9 

lUO 

630.2021 

609.5O64 

706.0404 

542.5 

1 

101 

502.6711 

595.2214 

587.1927 

532.2 

102 

411.8747 

498.5613 

596.4574 

511.1 

103 

297.8291 

293.8295 

559.8567 

490.4 

104 

180. 0420 

170.8505 

588.1182 

433.2 

105 

90.4322 

101.6695 

518. 0615 

389.2 

106 

39.2247 

61.0513 

454.1944 

335.1 

107 

19.0766 

35.5966 

343.8860 

275.6 

108 

7.0038 

24.4014 

255.2732 

220.8 

109 

3.0077 

14.8411 

138.3521 

155.5 

110 

1.7777 

10.9786 

74.6542 

79.37 

111 

I.8O92 

7.2152 

22.9635 

42.17 

112 

1.0008 

4.1761 

11.9365 

29.99 

113 

1.6046 

4.0181 

12.0233 

26.03 

114 

1.3736 

3.9115 

5.4930 

26.02 

115 

1.0990 

3.8252 

4.0780 

25.02 

116 

1.8525 

3.7743 

3.7349 

24.90 

117 

1.1002 

3.7600 

3.5823 

24.55 

118 

1.4236 

3.7812 

3.5265 

24.09 

119 

1.4176 

3.8413 

3.5885 

24.53 

120 

1.4091 

3.9397 

3.7287 

24.83 

121 

I.406I 

4.0615 

4.3421 

25.03 

122 

1.5075 

6.2118 

4.7449 

25.23 
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Table  (con't) 

A 

U -23 5th 

U-235  f.sp. 

Pu-239th 

T.  N. 

123 

1.6146 

9.U94 

5.6251 

25.96 

124 

1.7807 

13.1158 

6.9888 

30.09 

125 

2.1619 

20.9510 

9.4606 

42.12 

126 

3.1930 

34.3669 

16.3245 

79.04 

127 

13.0052 

47.2269 

38.4008 

155.6 

128 

37.4850 

71.1708 

84.5268 

220.9 

129 

80.2351 

126.9981 

169.0715 

277.9 

130 

198.8415 

213.6048 

268.7451 

334.2 

131 

293.1623 

366.3871 

375.8085 

389.9 

132 

438.0352 

503.8890 

520.7085 

434.5 

133 

658.0632 

553.6068 

684.5557 

491.0 

134 

803.3213 

590.6258 

744.2233 

512.3 

135 

646.9588 

614.6103 

708.4287 

531.8 

136 

641.2265 

639.7539 

657.9431 

542.5 

137 

614.9160 

640.8418 

659.2099 

545.6 

138 

573.6294 

629.9629 

628.5799 

5a. 6 

139 

652. a80 

6I4.5IO7 

584.1377 

532.4 

140 

642.5670 

594.2453 

559.7373 

522.1 

141 

643.8075 

585.0948 

530.9185 

509.8 

142 

611.8407 

562.0456 

497.8704 

497.8 

143 

602. 7ai 

535.2041 

454.27a 

472.5 

144 

557.5730 

508.7896 

388.4964 

ao.7 

145 

395.3681 

406.6189 

309.1452 

363.6 

146 

308.2233 

322.3120 

259.5084 

308.8 
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Table  IV  (con' t) 

A 

U -23  5th 

U-235  f.sp. 

Pu-239th 

T.N. 

147 

237.1949 

234.3551 

205.9887 

267.4 

148 

171.0584 

181.0381 

169.1711 

227.3 

149 

112.9712 

111.8032 

131.3564 

186.9 

150 

63.6805 

78.8216 

99.6577 

145.6 

151 

44.2537 

53.3984 

79.7632 

110.6 

152 

27.8603 

35.6057 

61.3763 

78.06 

153 

17.1238 

21.3461 

36.8504 

49.67 

154 

7.6331 

13. 4208 

28.9538 

22.49 

155 

3.2602 

7.1121 

22.9694 

10.48 

156 

1.4096 

2.54O8 

10.9796 

6.98 

157 

0.7729 

1.7284 

7.5613 

3.68 

158 

0.2015 

0.8382 

4. 2266 

2.02 

159 

0.1047 

0.3461 

2.0957 

1.06 

From  Ref  16:17-37 
Ref  6:7-15 
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APPENDIX  B 

Nuclear  Decay  Schemes 

This  appendix  contains  the  information  necessary  to  describe 
the  paths  by  which  the  isotopes  formed  in  fission  decay  and  the 
beta  decay  schemes  for  all  of  the  known  decays.  The  half-lives  for 
the  isotopes  are  given  in  seconds  with  those  isotopes  whose  half- 
lives  are  only  estimates  identified  by  an  asterisk.  As  described 
in  Chapter  III,  the  mass  differences  as  computed  by  the  Cameron  and 
Coryell  formulas  are  Included  for  use  in  computing  the  spectral 
shapes  for  those  isotopes  with  undetermined  beta  decay  schemes.  In 
the  description  of  the  beta  decay  schemes  the  first  number  given  is 
the  E.P.E.,  in  Mev,  of  the  beta  particle  and  the  number  in  paren¬ 
thesis  is  the  percent  of  the  total  number  of  isotopes  which  decay 
by  emitting  beta  particles  with  this  E.P.E.  In  the  decay  diagrams, 
the  vertical  arrow  connecting  two  isotopes  with  the  same  atomic 
number  indicates  an  isomeric  transition  and  this  is  indicated  by 
the  letters  IT  in  the  descriptions  of  the  decay  schemes.  When 
branching  occurs  in  the  decay  diagrams,  the  percentage  going  by  each 
route  is  indicated  by  the  small  number  close  to  the  decay  path  to 
which  it  applies.  Where  no  percentage  is  given  it  is  implied  to  be 
100  percent. 
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77. 


32^®l 


29CU - to^n 


34 


Se 


Half- 

-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Cu 

1.5  * 

Cu - f  Zn 

10.5 

9.63 

Zn 

3.5  "" 

Zn - f  Ga 

8.36 

7.92 

Ga 

15 

Ga  - >  Ge 

5.49 

5.01 

GSj^ 

54 

Ge - *  As 

3.11 

3.30 

Geg 

3.96x10* 

As - ►Se 

.7a 

.593 

As 

1.4x10® 

Seg 

stable 

Decay  Schemes: 


Cu - ►  Zn 

Zn - #•  Ga 

Ga - >  Gej^ 

Ga - ►  Geg 

Gej^ - »As 

Geg - ►■As 

As  - ^Se 


Undetermined 
Undetermined 
Undetermined  (50%) 

Undetermined  (50%) 

2. 7(36%),  2.9(48%),  IT  (36%) 

.38(3%),  .76(15%),  1.2(11%),  1.3(4%),  1- 
2.12(29%),  2.27(14%) 

.16(2.8%),  .438(2.7%),  .684(94.4%) 


56(24%), 
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78. 


39^^ - - 3^*^® - 3^^® - 3^^® 


Half- 

-lives  ( sec) : 

Mass 

Diff . (Mev) : 

Cameron 

Coryell 

Cu 

1.8 

Cu 

- ■>  Zn 

13.8 

13.6 

Zn 

2.5 

Zn 

- ^Ga 

6.72 

5.61 

Ga 

8  ^ 

Ga 

- *Ge 

8.97 

9.04 

Ge 

7.56x10^ 

Ge 

- *-As 

1.28 

1.03 

As 

5.46x10"^ 

As 

- ►Se 

4.50 

4.46 

Se 

stable 

Decay 

Schemes : 

Cu' - 

-^Zn 

Undetermined 

Zn 

-►Ga 

Undetermined 

Ga^ - 

-Ge 

Undetermined 

Ge 

->As 

.9  (10056) 

As - 

-Se 

4.1  (50%),  1 
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Half- 

■lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Cu 

1  * 

Cu - > 

■Zn 

12.1 

11.3 

Zn 

1.5  * 

Zn - . 

•  Ga 

10.0 

9.61 

Ga 

4.5  * 

Ga - . 

■Ge 

7.32 

6.74 

Ge 

25  * 

Ge - . 

•As 

4.77 

5.07 

As 

540 

As - ► 

Se 

2.68 

2.40 

234 

Se - »- 

■  Br 

.358 

.327 

Seg 

2.205x1012 

Br 

stable 

Decay  Schemes: 

Cu 

—  Zn 

Undetermined 

Zn 

■^Ga 

Undetermined 

Ga 

— Ge 

Undetermined 

Ge 

—  As 

Undetermined 

As  — 

—  Se, 

2.2  (100^) 

Seg 

-^Br' 

.160  (100%) 
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80. 


32^®  33'^®  34^® 


Half- 

-lives  ( sec) : 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Cu 

1  «- 

Cu 

^Zn 

16.7 

15.2 

Zn 

1.5  ^ 

Zn 

-^Ga 

8.37 

7.22 

Ga 

3  ^ 

Ga 

-*Ge 

10.6 

10.7 

Ge 

25 

Ge  — 

-^As 

3.13 

2.74 

As 

15 

As 

—  Se 

6.17 

6.22 

Se 

stable 

Decay  Schemes: 

Cu - ^  Zn 

Undetermined 

Zn - Ga 

Undetermined 

Ga  ^  Ge 

Undetermined 

Ge - ^  As 

Undetermined 

As - ^  Se 

Undetermined 
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81. 


Half- 

-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Zn 

1.5  ^ 

Zn 

— >Ga 

13.0 

11.3 

Ga 

2 

Ga  - 

— >■  Ge 

8.99 

8.54 

Ge 

7 

Ge- 

— >  As 

6.45 

6.93 

As 

33 

As 

— >  Se 

4«  54 

4.32 

Se^ 

3.42x10^ 

Se 

Br 

2.03 

2.31 

Seg 

1.08x10^ 

Br 

stable 

Decay  Schemes: 

Zn - ►  Ga 

Undetermined 

Ga - »-Ge 

Undetermined 

Ge  ^  As 

Undetermined 

As - *•  Se^ 

Undetermined  (50$) 

As - Seg 

Undetermined  (5*^) 

Se^ - ^  Br 

1.38  (100$) 
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82. 


30^^  3^^^ 


Half- 

•lives  ( sec) : 

Mass 

Diff .  (Mev) : 

Cameron 

Coryell 

Zn 

1 

Zn 

- ^  Ga 

12.3 

9.12 

Ga 

1.5  * 

Ca 

- >  Ge 

13.6 

12.7 

Ge 

5  ^ 

Ge 

- >  As 

4.81 

4.76 

As 

18  ^ 

As 

- ^-Se 

7.86 

8.32 

Se 

stable 

Decay  Schemes: 


Zn 

— >  Ga 

Undetermined 

Ga 

Ge 

Undetermined 

Ga 

— >  As 

Undetermined 

As 

Se 

Undetermined 
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83. 


Half- 

•lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Zn 

1  * 

Zn - >  Ga 

14.5 

12.6 

Ga 

1.5 

Ga - >  Ge 

12.9 

9.85 

Ge 

2 

Ge - >As 

9.44 

8.29 

As 

7  ^ 

As - >  Se 

6.22 

5.73 

Se, 

1.5x10^ 

Se - ^  Br 

3.74 

3.77 

Seg 

69 

Br - >Kr 

1.77 

1.41 

Br 

8.64xlCP 

Kr, 

6.84x10^ 

Krg 

stable 

Decay  Schemes: 

Zn  — 

Ga 

Undetermined 

Ga  — 

Ge 

Undetermined 

Ge 

As 

Undetermined 

As 

Se^^ 

Undetermined  ( 50^) 

As  — 

— >  Se^ 

Undetermined  (50)6) 

-^Br 

1.5  (100%) 

Se^ 

Br 

1.5  (10%),  3.4  (90%) 

Br  — 

-^Kr^ 

.91  (20%),  .96  (80%) 
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84. 


Half-lives  ( sec) : 


Mass  Diff.(Mev):  Cameron  Corvell 


Ga 

1  ^ 

Ga 

-^Ge 

15.1 

14.0 

Ge 

1.5 

Ge 

—As 

8.72 

6.11 

As 

4  ^ 

As 

—  Se 

10.9 

9.73 

Se 

180 

Se 

—  Br 

2.10 

1.83 

Br, 

3.60x102 

Br 

— ^Kr 

5.10 

5.45 

Brs 

1.92x102 

Kr 

stable 

Decay  Schemes: 

Ga 

— Ge 

Undetermined 

Ge 

—►As 

Undetermined 

As  — 

-^.Se 

Undetermined 

Se 

->Br, 

Undetermined  (50%) 

Se 

— Br^ 

Undetermined  {50%) 

Br,— 

— Kr 

.80  (20%),  1.9  (72%),  3.2  (8%) 

-»Kr 

.5  (3%),  .77  (19%),  .98  (2%),  1.39  (14%), 
1.81  (1  5%),  2.80  (15%),  3.83  (14%),  4.6 

66 


GNJi/Phy  s/ 6  5-3 


85. 


Half- 

-lives  ( sec) : 

Mass 

Diff . (Mev) : 

Cameron 

Corvell 

Ga 

.8 

Ga- 

- ►Ge 

14.8 

11.4 

Ge 

1  * 

Ga 

- ►As 

11.0 

9.84 

As 

.43 

As 

- -Se 

10.2 

7.32 

Se 

39 

Se 

- ►Br 

6.74 

5.40 

Br 

1.8x102 

Br- 

- ►Kr 

3.47 

3.08 

Kr, 

1.584x10^ 

Kr 

- ►Rb 

1.03 

1.16 

Kra 

3.28x10^ 

Rb 

stable 

Decay  Schemes: 

Ga - ►  Ge 

Undetermined 

Ge - As 

Undetermined 

As - -  Se 

Undetermined 

Se - ►Br 

Undetermined 

Br - ►Kr^ 

2.5  (10056) 

Kr^ - ►  Rb 

.824  (81%),  IT  (19%) 

Kra - ..Rb 

.15  (.7%),  .672  (99. 
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86. 


3  3*^® - - 3^^® - 


Half- 

■lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Coryell 

Ge 

1.5  ^ 

Ge 

— >  As 

10.7 

7.57 

As 

2 

As 

— ^  Se 

12.4 

11.2 

Se 

16 

Se 

-.Br 

6.03 

3.39 

Br 

54 

Br 

-^Kr 

8.12 

7.06 

Kr 

stable 

Decay  Schemes: 

Ge  - •►As 

As - ^  Se 

Se - ►Br 

Br - .Kr 


Undetermined 

Undetermined 

Undetermined 


3.0  (33%),  5.0  (33%),  7.1  (33%) 
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87. 


32 


Ge- 


33 


As- 


34 


Se- 


35 


Br- 


36 


Kr- 


37 


Rb 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Coryell 

Ge 

1  * 

Ge - ►  As 

13.2 

11.4 

As 

1.5  * 

As - ►  Se 

12.1 

9.17 

Se 

16 

Se - 1-  Br 

8.28 

7.31 

Br 

55 

Br - ».Kr 

7.41 

5.05 

Kr 

4.68x10^ 

Kr - »Rb 

4.05 

3.19 

Rb 

stable 

Decay  Schemes: 

Ge  — 

As 

Undetermined 

As  — 

Se 

Undetermined 

Se 

-^Br 

Undetermined 

Br 

-*Kr 

2.6  (70%),  8.0  (30%) 

Kr 

—  Rb 

1.25  (25%),  3.3  (10%) 

,  3.8  (65%) 
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88. 


32*^®  34^®  35®^  ^6^^  37  3*8^^ 


Half- 

-lives  ( sec) : 

Mass  Diff , (Mev) : 

Cameron 

Ge 

1  » 

Ge - 9- As 

11.8 

As 

1.5  ^ 

As - >Se 

14.6 

Se 

2.5  ^ 

Se - i-Br 

7.99 

Br 

16 

Br - ^Kr 

9.67 

Kr 

1.01X10“* 

Kr - >Rb 

3.35 

Rb 

1.08x10'^ 

Rb - a-Sr 

5.17 

Sr 

stable 

Decay 

Schemes : 

Ge 

*  As 

Undetermined 

As 

-  Se 

Undetermined 

Se 

*  Br 

Undetermined 

Br 

♦  Kr 

Undetermined 

Kr 

-  Rb 

.52  (70%),  .9  (10%),  2.7  (20%) 

Rb 

-*Sr 

.34  (6%),  2.5  (14%),  3.4  (4%), 

Corvell 

9.30 

13.0 

5.24 

8.97 

1.18 

4.19 


2  (7656) 
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89. 


33-^®  34^®  35®^  38^^  3^^ 


Half- 

-lives  (sec); 

Mass  Diff.(Mev): 

Cameron 

Corvell 

As 

1.5  * 

As - 

-.►Se 

13.2 

11.0 

Se 

2 

Se 

—  Br 

10.5 

9.14 

Br 

4.5 

Br 

—  Kr 

9.38 

6.93 

Kr 

192 

Kr 

—  Rb 

5.60 

5.12 

Rb 

900 

Rb 

-^Sr 

4. 48 

3. a 

Sr 

4.36x10® 

Sr 

—  Y 

1.40 

1.00 

Y 

stable 

Decay 

Schemes : 

As - 

*  Se 

Undetermined 

Se 

-  Br 

Undetermined 

Br 

►  Kr 

Undetermined 

Kr 

-Rb 

2.0  (3556),  3.9  (6556) 

Rb - 

♦  Sr 

.4  (256),  .67  (2856),  1.17  (356),  1.33  (256), 
1.61  (5356),  2.87  (5%),  3.92  (756) 

Sr 

-Y 

1.463  (10056) 
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90. 


34^6 - TsBr — i^Kr— ^Rb - ^Sr 


39 


40  ^^3 


Half- 

■lives  ( sec)  : 

Mass  Diff . (Mev) : 

Cameron 

Coryell 

Se 

1.5  ^ 

Se - ^Br 

9.13 

6.87 

Br 

1.6 

Br - .Kr 

11.9 

10.6 

Kr 

33 

Kr - .Rb 

5.32 

2.89 

Rb 

174 

Rb - .Sr 

6.74 

6.66 

Sr 

8.86x10® 

Sr - -Y 

.714 

-]  .09 

Y 

2.31x10® 

Y  - .Zr 

2.25 

2.68 

Zr^ 

.8 

21-2 

stable 

Decay  Schemes: 

Se 

-^Br 

Undetermined 

Br 

-^Kr 

Undetermined 

Kr  - 

—  Rb 

3.2  (100^) 

Rb 

—  Sr 

1.2  (15^),  2.21  (15^) 

,  4.4  (15%) 

,  5.81  (: 

6.59  (39%) 

Sr 

— Y 

.544  (100%; 

Y 

—  Zr^ 

2.27  (50%) 

Y 

-^Zrg 

2.27  (50%) 
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91. 


Half- 

-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Se 

1.5  * 

Se 

-*Br 

12.3 

10.7 

Br 

2 

Br 

*Kr 

10.5 

8.57 

Kr 

10 

Kr 

^Rb 

7.85 

6.80 

Rb 

72 

Rb 

—.Sr 

6.45 

5.13 

Sr 

3.49x10* 

Sr 

—  Y 

2.98 

2.76 

3x10^ 

Y 

—  Zr 

1.56 

1.19 

5.10x10® 

Zr 

stable 

Decay 

Schemes: 

Se 

-^Br 

Undetermined 

Br - 

-Kr 

Undetermined 

Kr 

-Rb 

3.6  (100^) 

Rb 

-Sr 

4.6  (100%) 

Sr 

1.09  (29.5%),  1.36  (29.5%) 

Sr 

*^3 

.61  (7%),  1.09  (3%),  2.03  (4%),  2.67  (27%) 

^3  — 

-Zr 

.33  (.3%),  1.54  (99.7%) 
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92. 


^Kr - 3i-Rb - ^Sr - - ^Zr 


Half- 

-lives  ( sec) : 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Br 

1.5  ^ 

Br- 

—  Kr 

13.7 

12.3 

Kr 

3 

Kr 

-^Rb 

6.48 

4.57 

Rb 

5.3 

Rb 

—  Sr 

9.00 

8.37 

Sr 

9.72x10^ 

Sr 

-  Y 

2.70 

.648 

Y 

1.27x10“ 

Y 

—  Zr 

3.83 

4.45 

Zr 

stable 

Decay  Schemes: 

Br 

—  Kr 

Undetermined 

Kr 

—  Rb 

Undetermined 

Rb 

— ^Sr 

Undetermined 

Sr 

— Y 

.545  (9056),  1.5  (10%) 

Y  — 

Zr 

1.26  (9%),  1.75  (3%), 
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93. 


36 


Kr- 


37 


Rb- 


30 


Sr— r,i- 


2S 


% 


40 


Zr 


7  5^ 


4iNb2 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Coryell 

Kr 

2 

Kr - >Rb 

9.67 

8.25 

Rb 

6 

Rb - »-Sr 

7.62 

6.61 

Sr 

498 

Sr - ^Y 

5.24 

4.27 

Y 

3.64x10* 

Y  - ^Zr 

3.55 

2.73 

Zr 

1 . 2x10^ 

Zr - yNb 

.499 

-.120 

Nb^ 

1.17x10® 

Nbg 

stable 

Decay  Schemes: 

Kr 

—  Rb 

Undetermined 

Rb - -  Sr 

Undetermined 

Sr 

Y 

Undetermined 

Y  - 

Zr 

.45  (.155^),  .71  (1.8^),  1.47  ( 
2.62  (3.9$),  2.89  (90$) 

.9$),  1.95(3$), 

Zr 

-^Nb. 

.034  (25$) 

Zr - ►Nbg 

.063  (75$) 
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94. 


36 


Kr- 


37 


Rb- 


38 


Sr- 


39^ 


40 


Zr 


Half-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Coryell 

Kr  1.5 

Kr 

-*Rb 

7.95 

6.34 

Rb  3 

Rb 

—Sr 

10.8 

10.2 

Sr  78 

Sr - 

-^Y 

3.87 

2.50 

Y  1.2x10^ 

Y  — 

— Zr 

6.10 

6.34 

Zr  stable 


Decay  Schemes: 


Kr  — ^ 

Rb 

Undetermined 

Rb - - 

Sr 

Undetermined 

Sr - ^ 

Y 

Undetermined 

Y  - ► 

Zr 

5.0  (100)6) 
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95. 


36 


Kr- 


Rb- 


"3^ 


Sr- 


39-' 


Half- 

•lives  (sec); 

Mass  Diff.(Mev):  Cameron 

Corvell 

Kr 

1.3  * 

Kr - .-Rb 

11.1 

9.79 

Rb 

2 

Rb - ►Sr 

9.10 

8.19 

Sr 

48 

Sr - .Y 

7.07 

5.89 

Y 

660 

Y  - ►Zr 

4.73 

4.39 

Zr 

5.62x10® 

Zr - ►Nb 

2.78 

2.09 

Nbj 

3.24x10® 

Nb - -Mo 

1.09 

.590 

Nbg 

3.02x10® 

Mo 

stable 

Decay  Schemes: 

Kr 

-..Rb 

Undetermined 

Rb 

—  Sr 

Undetermined 

Sr 

—  Y 

Undetermined 

Y 

-^Zr 

Undetermined 

Zr 

—  Nb^ 

.885  {2%) 

Zr 

— Nbs 

.360  (4256),  .396  (55^) 

,  1.13 

(156) 

Nbg 

-^Mo 

.160  .93  (1^) 

77 


GNE/Phys/65-3 


96. 


3  sKr  Rb  Sr  Y - Zr— - jgMo 


Half- 

■lives  ( sec) ; 

Mass 

Diff . (Mev) : 

Cameron 

Corvell 

Kr 

.8 

Kr 

- -Rb 

9.21 

7.66 

Rb 

1 

Rb- 

- ►Sr 

12.3 

11.5 

Sr 

2.5  * 

Sr- 

- -  Y 

5.35 

3.90 

Y 

138 

Y  - 

- ►  Zr 

7.93 

7.76 

Zr 

stable 

Zr- 

- »Nb 

i.a 

.138 

Nb 

8.28x10* 

Nb. 

- — -Mo 

3.63 

4.00 

Mo 

stable 

Decay  Schemes: 

Kr - -  Rb 

Undetermined 

Rb - -Sr 

Undetermined 

Sr - -  Y 

Undetermined 

Y  - *  Zr 

3.5  (10056) 

Zr - -  Nb 

No  transitions  -  Shielded  Nuclide 

Nb - ►Mo 

.37  (856),  .7  (9256) 
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97. 


Half- 

-lives  (sec): 

Mass  Diff.(Mev); 

Cameron 

Coryell 

Rb 

1 

Rb 

— ^Sr 

10.4 

9.89 

Sr 

1.5  ^ 

Sr 

-^Y 

8.53 

7.63 

Y 

5  * 

Y 

— .Zr 

6.22 

6.17 

Zr 

6 . 12x10* 

Zr 

— .Nb 

4.62 

3.91 

Nb^ 

60 

Nb 

— yMo 

2.28 

2.45 

Nba 

4.32x103 

Mo 

stable 

Decay  Schemes: 

Rb  ■ 

- -Sr 

Undetermined 

Sr  - 

- -Y 

Undetermined 

Y  - 

- ►  Zr 

Undetermined 

Zr- 

_ _  Nb^ 

.45  (10$),  1.91  (90$) 

Nba 

— .Mo 

.93  (1$),  1.267  (99$) 
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98. 


37 


Rb- 


38 


Sr 


39" 


40 


Zr- 


41 


Nb- 


43 


Mo 


Half- 

-lives  (sec): 

Mass  Diff • (Mev) : 

Cameron 

Corvell 

Rb 

.7  ^ 

Rb 

— Sr 

13.9 

13.2 

Sr 

1  ^ 

Sr 

—  Y 

6.64 

5.66 

Y 

2.5  ^ 

Y  - 

—  Zr 

9.40 

9.51 

Zr 

60 

Zr 

Nb 

2.91 

1.96 

Nb 

3.09x10^ 

Nb 

Mo 

5.48 

5.81 

Mo 

stable 

Decay  Schemes: 

Rb - »Sr 

Undetermined 

Sr - -Y 

Undetermined 

Y - ►Zr 

Undetermined 

Zr - ►Nb 

Undetermined 

Nb - ►Mo 

3.1  (100%) 
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99. 


,Sr — 


39“^  40  41 


Half- 

•lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Sr 

1  * 

Sr 

—  Y 

10.1 

9.33 

Y 

1.5  "" 

Y  — 

—  Zr 

7.51 

7.90 

Zr 

30 

Zr 

—  Nb 

6.09 

5.67 

Nb 

144 

Nb 

—  Mo 

3.78 

4.24 

Mo 

2.38x10® 

Mo 

-^Tc 

1.97 

2.01 

Tc, 

2.16x10* 

Tc 

-►Ru 

-.227 

.583 

TCg 

6.63x1013 

Ru 

stable 

Decav  Schemes: 

Sr - ►  Y 

Undetermined 

Y 

Zr 

Undetermined 

Zr 

Nb 

Undetermined 

Nb 

—  Mo 

3.2  (10056) 

Mo 

—  Tc^ 

.45  (4J6) 

Mo 

— ^ 

.45  (1056),  .87  (156),  1.23  (8556) 

—  Ru 

.292  (100%) 
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100. 


36 


Sr— t&Y 


40 


Zr- 


41 


Nb- 


4S 


Mo- 


43 


Tc- 


44 


Ru 


Half-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Sr  1  ^ 

Sr 

Y 

8.11 

7.11 

Y  1.5 

Y  — 

Zr 

11.0 

11.0 

Zr  3.5  ^ 

Zr 

—  Nb 

4.20 

3.49 

Nb  180 

Nb 

—  Mo 

6.97 

7.36 

Mo  stable 

Mo 

-^Tc 

.265 

-.125 

Tc  16 

Ru  stable 

Tc 

—  Ru 

2.99 

3.74 

Decav  Schemes: 

Sr 

- -  Y 

Undetermined 

Y 

- -  Zr 

Undetermined 

Zr 

- »Nb 

Undetermined 

Nb 

- -  Mo 

3.1  (45%),  3.5  (45%),  4.2  (10%) 

Mo 

- -Tc 

No  transition  -  Shielded  Nuclide 

Tc 

- -Ru 

2.2  (33%),  2.89  (33%),  3.37  (33% 
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101. 


38Sr— r,! - 4o2r“^iNb— :j-aMo-^Tc— 


Half-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Sr  .7 

Sr- 

— 

11.5 

10.6 

Y  1.0  * 

Y  — 

-^Zr 

9.00 

9.17 

Zr  2.5  ^ 

Zr 

— ►Nb 

7.72 

6.98 

Nb  60 

Nb- 

—  Mo 

5.08 

5.59 

Mo  876 

Mo 

—  Tc 

3.46 

3.40 

Tc  840 

Ru  stable 

Tc- 

— -  Ru 

1.28 

2.01 

Decay 

Schemes: 

Sr 

►  Y 

Undetermined 

Y 

►  Zr 

Undetermined 

Zr - 

-Nb 

Undetermined 

Nb - 

♦  Mo 

Undetermined 

Mo 

-Tc 

.6  (3%),  .7  (3856),  .8  (13%),  1.2  (11%), 
1.6  (25%),  2.23  (10%) 

Tc - 

-Ru 

1.07  (8%),  1.32  (92%) 

GNE/t’hys/65-3 


102. 


Half- 

■lives  ( sec) ; 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Y 

1.2  * 

Y  - 

--^Zr 

12.4 

12.3 

Zr 

2  ^ 

Zr 

-^Nb 

5.70 

4.83 

Nb 

7  ^ 

Nb 

— .Mo 

8.60 

8.69 

Mo 

660 

Mo 

- .Tc 

1.57 

1.27 

Tc, 

270 

Tc 

- .Ru 

4.48 

5.13 

Tca 

5 

Ru 

stable 

Decay  Schemes: 

Y- - -  Zr 

Undetermined 

Zr - *  Nb 

Undetermined 

Nb - .Mo 

Undetermined 

Mo  -.Tcj 

Undetermined  (50%) 

Mo - .TCg 

Undetermined  (50%) 

T  Cj - .  Ru 

2.0  (100%) 

Tcg - .  Ru 

4.1  (100%) 
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103. 


Half-lives  (sec);  Mass  Diff .(Mev) ;  Cameron  Coryell 


Y 

1.0 

Y - -Zr 

10.3 

10.4 

Zr 

1.5  ^ 

Zr - -Nb 

9.12 

8.27 

Nb 

4 

Nb - -Mo 

6.58 

6.91 

Mo 

25  ^ 

Mo - -Tc 

5.10 

4.75 

To 

72 

Tc - -Ru 

2.60 

3.39 

Ru 

3 . 46x10® 

Ru - ^Rb 

.755 

1.23 

Rh^ 

3.42x10® 

Rha 

stable 

Decay  Schemes: 

Y 

-►Zr 

Undetermined 

Zr  — 

-*Nb 

Undetermined 

Nb 

^Mo 

Undetermined 

Mo 

-►Tc 

Undetermined 

Tc 

-^Ru 

2.5  (10056) 

Ru 

-►Rh 

.10  (7%),  .212  (89%), 

.71  (3%) 
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104. 


ggY  - 


Half-lives  (sec): 

Mass  Diff.(Mev): 

• 

Cameron 

Corvell 

Y  1  * 

Y  — 

—  Zr 

13.5 

13.5 

Zr  1.7  * 

Zr 

—  Nb 

7.04 

6.12 

Nb  3  ^ 

Nb  — 

—  Mo 

10.0 

9.97 

Mo  60 

Mo 

—  Tc 

3.08 

2.62 

Tc  1.08X10^ 

Ru  stable 

Tc  — 

— »Ru 

6.12 

6.47 

Decay  Schemes: 

Y  - -Zr 

Undetermined 

Zr — *-Nb 

Undetermined 

Nb — »-Mo 

Undetermined 

Mo — -Tc 

Undetermined 

Tc - -Ru 

2.4  (10056) 
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Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron  Corvell 

Zr 

1.2  * 

Zr - ►Nb 

10.2 

9.51 

Nb 

2  * 

Nb - ►Mo 

7.93 

8.18 

Mo 

40 

Mo - ►Tc 

6.51 

6.05 

To 

480 

Tc - ►Ru 

4.11 

4.72 

Ru 

1. 59x10* 

Ru - ►Rh 

2. a 

2.59 

Rhj 

30 

Rh - ►Pd 

.317 

1.06 

Rhj, 

1.296x10® 

Pd 

stable 

Decay  Schemes; 

a 

Zr — ►Nb 

Undetermined 

Nb — ►Mo 

Undetermined 

Mo  — ►Tc 

Undetermined 

Tc - *-Ru 

Undetermined 

Ru  — ►  Rhj^ 

1.08  (5%) 

Ru  — » Rho 

.525  (5%),  .915  (10^) 

,  1.080  (25%) 

,  1.145 

1.870  (10^) 

Rhg — ►Pd 

.25  (10%),  .565  (90%) 
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106. 


40^^  43'^®  44^^  46^^  46^^ 


Half-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Zr  1  * 

Zr  - 

-^Nb 

8.29 

Nb  1.5  * 

Nb- 

-^Mo 

11.1 

Mo  U  * 

Mo  - 

— »  Tc 

4*  43 

Tc  9  ^ 

Tc 

Ru 

7.54 

Ru  3 . 16x10’^ 

Ru- 

— »  Rh 

.393 

Rh  30 

Pd  stable 

Rh 

— »Pd 

3.85 

Decay  Schemes: 

Sr _ *  Nb 

Undetermined 

Nb  — >Mo 

Undetermined 

Mo  — »Tc 

Undetermined 

Tc — -Ru 

Undetermined 

Ru — ^Rh 

.039  (100%) 

Rh — ^Pd 

1.2  (.456),  1.3  (.256),  1.5 

{.5%) 

2.4  (1156),  3.0  (856),  3.54 

(78$) 

Corvell 

7.38 

11.2 

3.94 

7.76 

.498 

4.32 


2.0  (2%), 
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107. 


^qZt - i'lNb  43^® 


Half-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Zr 

9.6  ^ 

Zr - ►Nb 

11.5 

10.5 

Nb 

1.0 

Nb - ►Mo 

9.19 

9.24 

Mo 

2.5  ^ 

Mo - ►  Tc 

7.62 

7.14 

Tc 

270 

Tc - ►  Ru 

5.47 

5.84 

Ru 

252 

Ru - ▼  Rh 

3.83 

3.74 

Rh 

1.30x10® 

Rh - ►Pd 

1.84 

2.24 

Pd 

2.21x10* 

Pd - ►Ag 

-0.102 

0.340 

Ag 

stable 

Decay  Schemes: 

Zr 

— ..Nb 

Undetermined 

Nb 

— ^Mo 

Undetermined 

Mo 

— ►Tc 

Undetermined 

Tc 

- ►Ru 

Undetermined 

Ru 

— ►Rh 

3.0  (1256),  3.8  (9^) 

,  4.0  (7956) 

Rh 

Pd 

.9  (356),  1.1  (13%), 

1.2  (8456) 

Pd 

- ►Ag 

.035  (10056) 
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108. 


4  4 


- 45^^ - 


Half- 

-lives  (sec); 

Mass  Diff • (Mev) : 

Cameron 

Corvell 

Nb 

1.0  ^ 

Nb - »  Mo 

12.4 

12.3 

Mo 

2.0  * 

Mo - *Tc 

5.70 

5.16 

Tc 

3.5 

T  c - »  Ru 

8.65 

8.96 

Ru 

270 

Ru - ►Rh 

1.76 

1.80 

Rh 

17 

Rh  — --Pd 

5.27 

5.60 

Pd 

stable 

Decay  Schemes: 

Nb — vMo 

Undetermined 

Mo - >Tc 

Undetermined 

Tc  — ►Ru 

Undetermined 

Ru  — ►Rh 

1.15  (28%),  1.32  (78%) 

Rh  — ►Pd 

3.5  (22%),  4.1  (175&), 
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109. 


Half-lives  ( sec) : 

Mass  Diff.(Mev)s 

Cameron 

Corvell 

Nb 

1.0  * 

Nb - -Mo 

10.4 

10.4 

Mo 

1.5  * 

Mo - -Tc 

8.93 

8.35 

Tc 

2.5  ^ 

Tc - -Ru 

6.74 

7.08 

Ru 

16  * 

Ru - -Rh 

4.95 

5.01 

Rhj 

50 

Rh - -Pd 

3.21 

3.54 

Rha 

30 

Pd - -Ag 

1.33 

1.67 

Pdi 

288 

Pda 

4.86x10“ 

Agi 

40 

stable 

Decay 

Schemes: 

Nb — » 

Mo 

Undetermined 

Mo - ► 

■Tc 

Undetermined 

Tc - ' 

■Ru 

Undetermined 

Ru - ' 

•Rh^ 

Undetermined  (50%) 

Ru - ’ 

'Rh- 

Undetermined  (50%) 

Rh- — <■ 

<o 

•Pdi 

2.6  (50%) 

Rho — ' 

•Pdo 

2.6  (50%) 

Pda — 

1.025  (100%) 
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no. 


41^^ - Ts^o - - j^Ru - ^Rh - jgPd 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) ; 

Cameron 

Corvell 

Nb 

0.8 

Nb  — 

— ^Mo 

13.4 

13.2 

Mo 

1.0  * 

Mo 

—  Tc 

6.95 

6.08 

Tc 

2.0  ^ 

Tc  — 

—  Ru 

9.97 

9.86 

Ru 

10  ^ 

Ru- 

-^Rh 

3.04 

2.76 

Rh 

3.6 

Rh  — 

-^Pd 

6.40 

6.54 

Pd 

stable 

Decay  Schemes: 


Nb  —  Mo 

Undetermined 

Mo  — ^Tc 

Undetermined 

Tc  — »■  Ru 

Undetermined 

Ru  —  Rh 

Undetermined 

Rh  —Pd 

Undetermined 
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111. 


Half-lives  ( sec) ; 

Mo  1.0 

Tc  1.5  * 

Ru  4  * 

Rh  12 

Pd^  1.98x10* 

Pdg  1.32x10^ 

Agj  74 

Agg  6.48x10® 

Cd  stable 


Decay  Schemes: 

Mo — ‘Tc 
Tc  — -Ru 
Ru — -Rh 
Rh — -Pd^ 

Rh  — *Pda 
Pd,— Ag, 
Pdg— Ag, 

Agg-^Cd 


Mass  Diff . (Mev) ; 

Mo - -Tc 

Tc - -  Ru 

Ru - .  Rh 

Rh - .Pd 

Pd - .Ag 

Ag - »  Cd 


Undetermined 
Undetermined 
Undetermined 
Undetermined  ( 50/6) 
Undetermined  ( 50^) 
.61  (3256),  IT  (6856) 
2.13  (10056) 

.69  (6.256),  .79  (1. 


Cameron 

Corvell 

9.96 

9.24 

7.99 

7.99 

6.27 

5.94 

4.49 

4.49 

2.46 

2.64 

.886 

.990 

,  1.05  (92.7%) 
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112. 


42 


Mo- 


43 


Tc 


44 


Ru- 


45 


Rh- 


45 


Pd- 


47  45^*^ 


Half-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Coryell 

Mo 

1  ^ 

Mo - ^Tc 

8.08 

6.99 

Tc 

1.5  ^ 

Tc - <“Ru 

11.0 

10.7 

Ru 

3 

Ru - -Rh 

4.30 

3.71 

Rh 

7  * 

Rh - -Pd 

7.73 

7.45 

Pd 

7.56x10* 

Pd  - -Ag 

.551 

.426 

Ag 

1.15x10* 

Ag - r  Cd 

4.08 

4.17 

Cd 

stable 

Decay 

Schemes: 

Mo — - 

Tc 

Undetermined 

Tc — • 

Ru 

Undetermined 

Ru — ► 

Rh 

Undetermined 

Rh — . 

Pd 

Undetermined 

Pd — . 

•Ag 

.28  (100%) 

Ag  — » 

•  Cd 

.92  (.8%),  1.22  (2.6%),  1.34  (3.6%), 

1.50  (1.8%),  1.61  (2%),  1.78  (4.4%), 

2.63  (.95f),  2.73  (2.8%),  3.42  (18%),' 
4.04  (56%) 
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113. 


Half- 

-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

Tc 

1.2  ^ 

Tc 

Ru 

9.14 

5.30 

Ru 

2 

Ru— 

Rh 

7.31 

6.80 

Rh 

4  “ 

Rh  — 

—  Pd 

5.76 

5.38 

Pd 

90 

Pd 

Ag 

3.79 

3.56 

Agi 

72 

Ag  — 

Cd 

2.18 

1.94 

1.91x10* 

Cd  — 

—  In 

.253 

.324 

Cd^ 

4.41x10® 

Cdg 

stable 

In 

stable 

Decay  Schemes: 

Tc - Ru 

Undetermined 

Ru  Rh 

Undetermined 

Rh - Pd 

Undetermined 

Pd - Ag^ 

3.3  (100^) 

Aga  —  Cdg 

2.0  (10056) 

Cdj^ - In 

.575  (100%) 
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114. 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Tc 

1.0  * 

Tc  - 

— >Ru 

12.4 

11.5 

Ru 

1.5  ^ 

Ru  - 

— >  Rh 

5.45 

4.58 

Rh 

3  - 

Rh  - 

— 5.  P  d 

8.77 

8.29 

Pd 

144 

Pd  - 

— >  Ag 

1.82 

1.36 

Ag 

5 

Ag 

— »  Cd 

5.42 

5.07 

Cd 

stable 

Decay  Schemes: 


Tc 

- ►Ru 

Undetermined 

Ru 

- ^Rh 

Undetermined 

Rh 

- ►Pd 

Undetermined 

Pd 

- ►  Ag 

1.4  (10056) 

Ag 

- >Cd 

4.6  (10056) 

96 


GNE/Phys/ 65-3 


115. 


Half-lives  ( sec) : 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

Tc 

0.8 

Tc - »Ru 

10.5 

9.62 

Ru 

1  ^ 

Ru - ^Rh 

8.67 

7.43 

Rh 

2.5  ^ 

Rh - >Pd 

6.91 

6.04 

Pd 

45 

Pd - *kg 

4*  84 

4.45 

Agi 

20 

Ag - >.  Cd 

3.45 

2.86 

Aga 

1 . 26x10^ 

Cd - ►In 

1.60 

1.27 

Cd, 

3.72x10® 

In - *  Sn 

.158 

.182 

Cda 

1.987x10® 

In^ 

1 . 58x10^ 

stable 

Sn 

stable 

Decay  Schemes: 

Tc 

- »Ru 

Undetermined 

Ru- 

— >Rh 

Undetermined 

Rh- 

— »Pd 

Undetermined 

Pd- 

— ^•Ag, 

Undetermined  (28^) 

Pd  - 

1 

- ’■Agy 

Undetermined  (72^) 

Ag - 

— >Cd 

Undetermined  (28^), 

IT  {72%) 

Agg- 

-^Cd 

2.75  i9%) 

Ag  ■ 

— »  Cd^ 

2.9  (91^) 

Cd- 

In 

.2{.3%),  .335(1^),  . 

687(2%), 

1.63(96.7%) 

Cdo- 

— ^In, 

.59(24^),  .63(12.5^) 

,  .85(1%) 

Cdg- 

— *  li^s 

1. 11(62. 5/&) 

In,- 

— 5.  Sn 

.84(55^),  IT(95^) 
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116. 


44^^  — JfePd - J^Ag  — i?Cd 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Ru 

1.5  * 

Ru  — 

— >Rh 

6.77 

5.49 

Rh 

2  » 

Rh  — 

— •‘Pd 

10.1 

9.14 

Pd 

10 

Pd  — 

— >  Ag 

2.98 

2.31 

Ag 

150 

Ag 

Cd 

6.48 

5.96 

Cd 

stable 

Decay  Scheme; 


Ru- 

— ^Rh 

Undetermined 

Rh- 

— >Pd 

Undetermined 

Pd  - 

— ^Ag 

Undetermined 

Ag 

--►Cd 

5.0  (100%) 
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117. 


4*Ru 


Half-lives  (sec):  Mass  Diff.(Mev):  Cameron  Corvell 


Ru 

1.0  * 

Ru - ►Rh 

9.81 

8.53 

Rh 

1.5  ^ 

Rh - ►Pd 

8.24 

7.15 

Pd 

5  ^ 

Pd - ►Ag 

6.21 

5.37 

Ag 

66 

Ag - ►Cd 

4.62 

3.79 

Cd, 

1.15x10* 

Cd - >  In 

2.66 

2.21 

Cda 

1.01x10* 

In - ►Sn 

1.44 

1.13 

In, 

6.84x10^ 

2.7x10^ 

Sn 

stable 

Decay  Schemes: 

Ru  — 

-s-Rh 

Undetermined 

Rh  - 

-»-Pd 

Undetermined 

Pd  - 

— ►Ag 

Undetermined 

Ag 

-^Cd, 

Undetermined  (5056) 

Ag 

-^Cda 

Undetermined  (50)6) 

Cd  — 

1.0  (10056) 

Cdg 

—  In, 

1.8  (7056),  2.3  (3056) 

In, 

'  Sn 

.95  (456),  1.62  (2456) 

,  1.77  (5256) 

,  IT  (2056) 

1^3 

— ^  Sn 

.74  (10056) 
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118. 


Half- 

•lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Ru 

1  » 

Ru 

— -Rh 

7.80 

6.37 

Rh 

1.5 

Rh 

— -Pd 

11.3 

9.96 

Pd 

4.5  ^ 

Pd 

— -Ag 

4.32 

3.23 

Ag 

25  ^ 

Ag- 

Cd 

7.85 

6.82 

Cd 

3x10-^ 

Cd 

—  In 

.801 

.089 

Ini 

270 

In 

— Sn 

4*47 

3.68 

Ing 

5 

Sn 

stable 

Decay  Schemes? 


Ru — »Rh 

Undetermined 

Rh — -Pd 

Undetermined 

Pd  — -Ag 

Undetermined 

Ag  — Cd 

Undetermined 

Cd - -Ing 

.8  (10056) 

lUi — -Sn 

1.5  (100^) 

lug  — »Sn 

3.3  (20^),  U. 
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119. 


46 


Half-lives  (sec): 

Mass  Diff.(Mev)s 

Cameron 

Rh 

1  * 

Rh  — 

—Pd 

9.28 

Pd 

3  ^ 

Pd  — 

— Ag 

7.36 

Ag 

17  * 

Ag 

-^Cd 

5.96 

Cd, 

142 

Cd  — 

—  In 

4.04 

1 

Cdr, 

570 

In 

— r  Sn 

2.62 

Irij  1.08x10^ 
Ina  120 
Siij^  2 . 16x10’’’ 
Sng  stable 


Decav  Schemes J 

Rh - 'Pd 

Undetermined 

Pd - -Ag 

Undetermined 

Ag - ^Cd, 

Undetermined  (50%) 

o  1 

Ag - ►Cdp 

Undetermined  ( 50%) 

Cdj - ►Ing 

Cdo — ►In, 

Undetermined 

3.5  (80%),  4.0  (20%) 

Ihj^ - ►Sng 

1.8  (6%),  2.7  (90%),  IT  (4%) 

Ina — ►Snj^ 

1.6  (100%) 

Corvell 

8.00 

6.24 

4.68 

3.12 

2.06 
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120. 


Half-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Rh 

1.5  * 

Rh 

- ^  Pd 

12.4 

10.7 

Pd 

2 

Pd 

- »  Ag 

5.36 

4.07 

Ag 

6  ^ 

Ag 

- »  Cd 

9.00 

7.63 

Cd 

60 

Cd 

- -In 

2.15 

.992 

In, 

Sn 

3 

50 

stable 

In 

- ^  Sn 

5.85 

4.55 

Decay  Schemes: 

Rh — -Pd 

Undetermined 

Pd - -Ag 

Undetermined 

Ag_Cd 

Undetermined 

Cd  — -In, 

Undetermined  (50%) 

Cd - -Ing 

Undetermined  (50%) 

Ing — »Sn 

2.2  (100%) 
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121. 


46 


Rh 


tePd - J-,Ag 


Half-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Rh  1  * 

Rh - ►Pd 

10.8 

Pd  1.5  * 

Pd - ►Ag 

8.44 

Ag  4  * 

Ag - -Cd 

7.01 

Cd  210 

Cd - ►  In 

5.19 

Ini 

In - ►  Sn 

3.97 

Ing  30 

Sni  8.88x10® 

Sng  9.0x10* 

Sb  stable 

Sn - ►  Sb 

.427 

Decay  Schemes: 

Rh — ►Pd 

Undetermined 

Pd— Ag 

Undetermined 

Ag— Cd 

Undetermined 

Cd — ►Ini 

Undetermined  ( 50%) 

Cd  — ►Ing 

Undetermined  (5C^) 

Ini — ►Sni 

3.7  (100%) 

Ing  — Sng 

2.9  (100%) 

Sni — ►  Sb 

.42  (100%) 

Sng  — >•  Sb 

.383  (100%) 

Corvell 

8.77 

7.04 

5.51 

3.98 

2.95 

.918 
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122. 


46^^^  ~ 


Half- 

-lives  (sec): 

Mass 

Diff . (Mev) : 

Cameron 

Coryell 

Pd 

1.5  * 

Pd 

- ».Ag 

6.92 

4.94 

Ag 

3  ^ 

Ag 

- vCd 

10.1 

8.44 

Cd 

40  ^ 

Cd 

- ^In 

3.20 

1.90 

In 

7.5 

In 

- V  Sn 

7.02 

5.40 

Sn 

stable 

Decay  Schemes; 


Pd - *kg 

Undetermined 

Ag  — *Cd 

Undetermined 

Cd - *In 

Undetermined 

In  — »Sn 

Undetermined 
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123. 


Half-lives  ( sec) : 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Pd  1.5  * 

Pd 

-^Ag 

9.32 

7.85 

Ag  2.5  ^ 

Ag¬ 

— ►Cd 

8.57 

6.34 

Cd  9  * 

ed 

— ^In 

6.28 

4.83 

In,  36 

In 

-^Sn 

5.02 

3.82 

Ing  10 

Sn  — 

— >■  Sb 

1.59 

1.81 

Sn^  1.08x10’^ 
Sng  2.40x10^ 
Sb  stable 


Decay  Schemes; 


Pd— Ag 

Undetermined 

Ag  — Cd 

Undetermined 

Cd - ►In, 

Undetermined  (5056) 

Cd  - 

Undetermined  (50%) 

In,— Sng 

4.6  (10056) 

Ing — »Sng 

3.3  (10056) 

— ^Sb 

.34  (2%),  1.42  (9856) 

Sn^  ^  Sb 

1.26  (10056) 
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124. 


Pd  — ^Ag - - i^In - g^oSn 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Coryell 

Pd 

1 

* 

Pd 

— -Ag 

7.79 

5.83 

Ag 

2 

Ag 

— ..  Cd 

11.0 

9.26 

Cd 

10 

* 

Cd 

— ♦  In 

4.77 

2.81 

In 

20 

* 

In 

— V  Sn 

8.11 

6.24 

Sn  stable 


Decay  Schemes: 


Pd 

—  Ag 

Undetermined 

Ag 

_ ,Cd 

Undetermined 

Cd 

— »In 

Undetermined 

In 

— Sn 

Undetermined 
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125. 


Half-lives  (sec): 

Mass  Diff.(Mev):  Cameron 

Corvell 

Ag 

2 

Ag - ►Cd 

9.45 

7.20 

Cd 

5  ^ 

Cd - •In 

7.18 

5.70 

In 

15  ^ 

In - ►Sn 

6.60 

4.70 

Sn^ 

582 

Sn - .Sb 

2.69 

2.70 

Sng 

8.12x10® 

Sb - ►Te 

.837 

.800 

Sb 

8 . 51x10'^ 

Te, 

5.01x10® 

stable 

Decay  Schemes: 

Ag-^ 

Cd 

Undetermined 

Cd — ' 

■In 

Undetermined 

In — . 

Snj^ 

Undetermined  (50/6) 

In — . 

Sng 

Undetermined  (50%) 

Sn^ — ► 

Sb 

.65  (2.2%),  2.04  (97.8%) 

Sng — » 

Sb 

.37  (2.1%),  .47(1.5%), 

.93  ( 

.1%),  1.3  (1.3%), 

2.33  (95%) 

Sb — ► 

Te, 

,09  (.5%),  .118  (1.3%), 

.124  (6.4%),  .295  (8.6%i 

.437  (1.3%),  .612  (2.9%) 

Sb  — ^ 

Teg 

.09  (1.7%),  .118  (4.8%) 

1,  .124  (24.6%), 

.295  (32.5%),  .437  (4.7%), 

612  (10.7%) 
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126. 


Half- 

-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Ag 

1.5  * 

Ag 

— ^  Cd 

11.9 

10.0 

Cd 

3  ^ 

Cd- 

-^In 

5.65 

3.71 

In 

7  * 

In 

— >-  Sn 

9.01 

7.02 

Sn 

6.30x1013 

Sn— 

—  Sb 

1.19 

.729 

Sbj 

1.14x103 

Sb 

-♦  Te 

3.93 

4.04 

Sbg 

1.08x10® 

Te 

stable 

Decay  Schemes: 

Ag — ►  Cd 

Undetermined 

Cd - ►In 

Undetermined 

In  — >  Sn 

Undetermined 

Sn - ►Sbj 

Undetermined 

Sbg— >-Te 

1.9  (10056) 
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127. 


*7 


Ag 


48 


Cd - :^ln 


Half- 

-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Ag 

1  ^ 

Ag - -  Cd 

10.3 

7.99 

Cd 

1.5 

Cd - »  In 

8.10 

6.51 

In 

3  * 

In - ►  Sn 

7.48 

5.53 

Sn 

9.36x10^ 

Sn - *  Sb 

3.60 

3.55 

Sb 

3.2x10^ 

Sb - ►  Te 

2.42 

1.67 

Te, 

9.07x10® 

Te - >  I 

.738 

.592 

Tea 

3.35x10* 

I 

stable 

Decay  Schemes: 


Ag 

Cd 

In. 

Sn 

Sb- 

Sb 

Te, 

Teo 


►Cd 
'In 
Sn 
Sb 
Te, 
■Tea 
►  I 
»I 


Undetermined 
Undetermined 
Undetermined 
Undetermined 
1.50  (5056) 

.80  (2256),  .86  (656),  1.11  (2256) 
.73  (1.5%),  IT  (98.5^) 

.27  {1%),  .695  (9956) 
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128. 


48 


Cd- 


49 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Cd 

1  ^ 

Cd - ►In 

6.50 

4.56 

In 

2 

In - >  Sn 

9.94 

7.79 

Sn 

3.72x10'* 

Sn - >  Sb 

2.07 

1.62 

Sb^ 

3 . 46x10'* 

Sb - »  Te 

4.84 

4.85 

Sbg 

600 

Te 

stable 

Decay  Schemes: 

Cd - »In 

Undetermined 

In - »Sn 

Undetermined 

Sn  — *Sbj^ 

Undetermined  (3%) 

Sn - Sbg 

Undetermined  (97^ 

Sbj — ►Te 

1.0  (10056) 

Sbg-^  Te 

2.9  (100%) 
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129. 


Half-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Coryell 

Cd 

1  ^ 

Cd - - 

In 

8.95 

7.30 

In 

1.5  * 

In - » 

Sn 

8.34 

6.34 

Sn 

3.60x10® 

Sn - r 

Sb 

A.  53 

4.38 

Sb 

1.55x10* 

Sb - r 

Te 

3.32 

2.52 

Te, 

2.85x10® 

Te - ► 

I 

1.65 

1.46 

Teg  4.02x10® 
I  stable 


Decay  Schemes: 


Cd — »  In 
In  — >  Sn 
Sn — ^  Sb 
Sb— »  Te^ 
Sb — »  Te- 
Teg-^I 


Undetermined 

Undetermined 

Undetermined 

1.87  (4.856),  3.32  (19.2^) 

1.87  (15%),  3.32  (61%) 

.3  (1%),  .7  (3%),  .99  (16%),  1.453  (80%) 


in 


GNE/Phys/65-3 


130. 


48Cd 


49 


In 


54 


Xe 


Half- 

-lives  ( sec) : 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Cd 

1  ^ 

Cd 

-^In 

8.67 

12.9 

In 

1.5 

In 

Sn 

10.8 

8.71 

Sn 

156 

Sn  — 

-»■  Sb 

2.94 

2.67 

Sb^ 

420 

Sb  — 

-^Te 

5.77 

5.79 

Sbg 

1.98x10^ 

Te - »  I 

.130 

7.05 

Te 

stable 

I 

—  Xe 

2.71 

1.02 

I 

4. 5x10^ 

Xe 

stable 

Decay  Schemes: 

Cd - -  In 

Undetermined 

In  — ►  Sn 

Undetermined 

Sn  — ^  Sbj^ 

Undetermined  (50/f) 

Sn - ^  Sbg 

Undetermined  (50%) 

Sb^—  Te 

Undetermined 

Sbg — '■Te 

Undetermined 

Te  — ►  I 

No  Transition  -  Shielded  Nuclide 

I  — »  Xe 

.6  (47%),  1.02  (53%) 
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131. 


Half- 

-lives  (sec): 

Mass  Diff.(Mev) : 

Cameron 

Coryell 

In 

1.5  ^ 

In  — 

—  Sn 

10.5 

7.46 

Sn 

204 

Sn 

-*Sb 

5.39 

5.51 

Sb 

1380 

Sb  — 

-^Te 

4.18 

3.66 

Te^ 

1.04x10® 

Te 

2.59 

2.61 

Tea 

1.5x10® 

I  — 

—  Xe 

1.19 

1.16 

I  6.96x10® 
1.04x10® 
XSg  stable 

Decay  Schemes: 


In - *Sn 

Sn — >Sh 
Sb — ►Tej 
Sb — ^Tea 
T9,^I 

Tea — *1 

I  .Xej 

I  ►Xea 


Undetermined 
Undetermined 
Undetermined  (15%) 

Undetermined  (85$) 

.215  (3.6$),  .420  (43$),  .570  (30$), 

2.457  (4.4$),  IT  (19$) 

1.15  (10$),  1.36  (5$),  1.68  (25$),  2.14  (60$) 
.25  (1.4$),  .33(4.6$),  .606  (  43.7$),  .81  (.3$) 
.25  (1.4$),  .33  (4.6$),  .606  (43.7$),  .81  (.356) 
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49 


In- 


>«Sn- 


60 


51 


Sb- 


>,Te- 


52 


53 


54 


Half-lives  (sec): 

Mass  Diff.(Mev):  Cameron 

Coryell 

In  1.5  ^ 

In  — — '  Sn 

13.5 

9.44 

Sn  132 

Sn - ►Sb 

5.12 

3.52 

Sb  126 

Sb - .Te 

6.63 

6. 56 

Te  2.81x10® 

Te - -I 

1.00 

.640 

I  8.28x10^ 

Xe  stable 

I  - -Xe 

3.65 

3.68 

Decay  Schemes: 

In - ^Sn 

Undetermined 

Sn - ►Sb 

Undetermined 

Sb  — ►  Te 

Undetermined 

Te  —  I 

.22  (100%) 

I  — ►Xe 

.80  (21%),  1.04  (15%), 
1.61  (21%),  2.14  (18%) 

1.22  (12%) 

,1.49  (12%), 
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133. 


Half- 

-lives  ( sec)  : 

Mass  Diff.(Mev)s 

Cameron 

Corvell 

Sn 

1.8 

Sn - ►Sb 

8.07 

6.29 

Sb 

246 

Sb - ►Te 

6.37 

4.46 

Te, 

3.0x10^ 

Te - ►! 

3.46 

3.43 

Tea 

120 

I  ^  Xe 

2.06 

2.00 

I 

7.56x10* 

Xe - ►  Cs 

.310 

.572 

Xe^ 

1.99x10® 

Xea 

4.55x10® 

Cs 

stable 

Decay  Schemes: 

Sn - ►Sb 

Undetermined 

Sb  — ►  Te^ 

3.67  (92%) 

Sb  — ►Te- 

4.00  (8%) 

Te^^I 

1.3  (705t),  2.4  (17%),  IT  (13%) 

Te--^  I 

2.4  (100%) 

I  - ^Xe^ 

1.54  (4%) 

I  — ‘Xeg 

.5  (5.8%),  1.22  (90.2%) 

Xe^  — *  C  s 

.268  (1%),  .347  (99%) 
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134. 


6oSn— ^Sb— ^Te— ^3! - g^Xe 


Half-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Sn  1.3 

Sn - '■Sb 

8.04 

4.63 

Sb  48 

Sb  - -Te 

9.32 

7.59 

Te  2.52x10"^ 

Te - »I 

3.19 

1.79 

I  3.18x10^ 

I  - ^Xe 

4.52 

4.75 

Xe  stable 

Decay  Schemes s 

Sn - *  Sb 

Undetermined 

Sb  — ►Te 

Undetermined 

Te  — ►! 

Undetermined 

I  — »Xe 

.5  (6.556),  1.05  (1^),  1.25  (2356),  1.49  (1556), 

1.68  (7.556),  1.81  (9.556),  2.21  (1256), 

2.41  (2556) 
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135. 


Half-lives  (sec); 


Sn 

1  ^ 

Sb 

24 

Te 

60 

I 

2.4x10^ 

Xe 

3.31x10* 

Csj^ 

3.18x10-' 

CSg 

6.3x10^-' 

Ba 

stable 

Decay  Schemes: 

Sn  — Sb 
Sb — .Te 
Te  —  I 

I  - -Xe 

Xe  — 

Xe  ^  ^^2 
Csj^ — *•  Ba 
C  Sg — »  Ba 


Mass  Diff .  (Mev) ; 

Sn - 'Sb 

Sb - ^Te 

Te - I 

I  - Xe 

Xe - ►  Cs 

Cs - ►  Ba 


Undetermined 
Undetermined 
Undetermined 
.47  (35^),  1.0  (40%) 
Undetermined  ( 50%) 
Undetermined  (50%) 
Undetermined  (16%), 
.21  (100%) 


Cameron 

Corvell 

9.75 

7.33 

9.29 

5.92 

6.15 

4.51 

4.26 

3.10 

1.18 

.282 

-.169 

-1.13 

,  1.4  (25%) 
IT  (84%) 
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136. 


ei 


Sb- 


^oTe- 


52 


53- 


66*^® 


Half-lives  (sec): 

Sb  3  * 

Te-  6 

I  83 

Xe  stable 

Cs  1.2x10® 

Ba  stable 


Mass  Diff.(Mev): 

Sb - -Te 

Te - '  I 

I  - -  Xe 

Xe - -Cs 

Cs - ^Ba 


Cameron 

Corvell 

11.0 

8.58 

6.12 

2.90 

7.22 

5.78 

.924 

.100 

2.29 

2.98 

Decay  Schemes: 

Sb - -  Te 

Te  — .  I 
I  - ,Xe 

Xe — -Cs 
Cs  — -Ba 


Undetermined 

Undetermined 

2.73  (6%),  3.7  (445^),  4-16  (5^),  4-37  (24%), 
5.60  (15%),  7.00  (6%) 

No  Transitions  -  Shieldeii  Nuclide 
.341  (92.6%),  -  .657  (7.4%) 
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61 


Sb 


6^'^® 


136 

Xe+n 

f 

64 


”^4^®‘ 


96% 


^■=Cs 


/ 


66 


66^63 


Half- 

-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Sb 

1.5  * 

Sb - - 

Te 

10.4 

6.55 

Te 

3  ^ 

Te - * 

I 

7.83 

5.56 

I 

24 

I  - ^ 

Xe 

7.19 

4.17 

Xe 

234 

Xe - ^ 

Cs 

3.89 

2.78 

Cs 

9.45x106 

Cs - » 

•  Ba 

2.03 

1.39 

Baj  156 
Bag  stable 


Decay  Schemes: 


Sb — ‘Te 

Undetermined 

Te — ‘I 

Undetermined 

I  — ‘Xe 

Undetermined 

Xe  — »Cs 

3.5  (96%) 

Cs  — ►Ba^ 

.514  (95%) 

Cs  — 

1.18  (5%) 
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138. 


137 

g^Xe+n 

1^ 


52 


Te— 


54 


Xe- 


99 


% 


55 


Cs- 


56 


iBa 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Te 

2 

Te  —I 

7.19 

4.07 

I 

6.3 

I  - -Xe 

8.91 

6.88 

Xe 

1.02x10^ 

Xe - ♦  Cs 

3.86 

1.33 

Cs 

1.93x10^ 

Cs - -  Ba 

5.00 

4.14 

Ba 

stable 

Decay  Schemes; 


Te - •  I 


Undetermined 


I  — »Xe 


Undetermined 


Xe — >  Cs 
Cs - »  Ba 


2-4  995S) 

2.2  (lO%j,  2.39  (36$),  2.53  (5$),  2.62  (16$) 
2.94  (12$),  3.40  (21$) 
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139. 


54 


138 

Xe+n 


53^®  53^ 


64 


small% 
4  00% 


le¬ 


ss 


Cs 


66 


Ba- 


67 


La 


Half-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Te  1.3  * 

Te - *1 

8.97 

6.80 

I  2 

I  - »Xe 

8.26 

5.44 

Xe  41 

Xe - ►  Cs 

5.58 

4.08 

Cs  570 

Cs - » Ba 

4.97 

2.72 

Ba  4.98x10^ 

Ba - »La 

1.77 

1.36 

La  stable 

Decay  Schemes: 

Te - »I 

Undetermined 

I  - .Xe 

Undetermined 

Xe  — >Cs 

3.5  (50%),  4.6  (50%) 

Cs  — ►  Ba 

2.7  (16%),  3.4  (4%), 

4.0  (80%) 

Ba  — ^La 

2.17  (28%),  2.34  (72%) 
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140. 


53 


Te- 


64 


Xe- 


56 


Cs- 


56 


Ba- 


07 


La- 


6a 


Ce 


Half- 

-lives  ( sec) : 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Te 

1  * 

Te 

7.80 

5.37 

I 

1.8  ^ 

I  — 

—  Xe 

10.0 

8.12 

Xe 

16 

Xe  — 

—  Cs 

4*94 

2.67 

Cs 

66 

Cs  — 

Ba 

6.69 

5.42 

Ba 

1.11x10® 

Ba 

—  La 

1.75 

-.025 

La 

1.45x10® 

La  — 

—  Ce 

3.17 

2.72 

Ce 

stable 

Decay  Schemes: 

Te — •  I 

Undetermined 

I  - .Xe 

Undetermined 

Xe - »Cs 

Undetermined 

Cs  — » Ba 

6.0  ’(10056) 

Ba - ►La 

.48  (25^),  .6  (1056),  .9  (556),  1.01  (6056) 

La  ^  C  e 

.83  (1256),  1.10  (2656),  1.38  (4556),  1.71 
2.20  (756) 
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141. 


63  ^  6*^®  ^6^®  5%^®  3?^® 


Half- 

-lives  (sec); 

Mass  Diff.(Mev): 

Cameron 

Corvell 

I 

1.5  ^ 

I 

—  Xe 

8.88 

6.70 

Xe 

2 

Xe  — 

—  Cs 

6.72 

5.36 

Cs 

24 

Cs 

Ba 

6.06 

4.02 

Ba 

1.08x10^ 

Ba 

-  La 

3.47 

2.68 

La 

1.40x10* 

La — 

-  Ce 

3.15 

1.69 

Ce 

2.81x10® 

Ce 

-  Pr 

-.163 

.0 

Pr 

stable 

Decay  Schemes: 

I - •  Xe 

Undetermined 

Xe  —  Cs 

Undetermined 

Cs  — r  Ba 

Undetermined 

Ba  — *La 

2.8  (100%) 

La  — -f  Ce 

1.06  (2%),  2.43,(98%) 

0 

CD 

1 

.435  (70%),  .580  (30% 
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142. 


53I - 6-4^® - 


Half-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

I  1  ^ 

I  - ►  Xe 

10.6 

9.31 

Xe  1.5 

Xe - ►  Cs 

5.56 

3.98 

Cs  2.3 

Cs  - — Ba 

7.84 

6.65 

Ba  660 

Ba - ►  La 

2.83 

1.32 

La  5.04x10'=' 

La - ►  Ce 

4.88 

3.99 

Ce  stable 

Decay  Schemes: 

I  — ►Xe 

Undetermined 

Xe  —  Cs 

Undetermined 

Cs  — ►  Ba 

Undetermined 

Ba  — «■  La 

4.0  (100%) 

La  — ^  Ce 

2.5  (30%),  4.0  (30%) 

,  4.5  (40%) 
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143. 


64^® TbCs - ^Ba — -^La — -igPr - ^oNd 


Half- 

-lives  (sec): 

Mass 

Diff . (Mev) : 

Cameron 

Corvell 

Xe 

1 

Xe 

- -  Cs 

7.33 

6.60 

Cs 

2 

Cs 

- ^Ba 

6.68 

5.28 

Ba 

12 

Ba 

- >  La 

4.62 

3.96 

La 

840 

La 

- ^Ce 

4.24 

2.99 

Ce 

1.18x10® 

Ce 

- »  Pr 

1.54 

1.32 

Pr 

1.18x10® 

Pr 

- »  Nd 

1.57 

.350 

Nd 

stable 

Decay  Schemes: 

Xe  — '■Cs 

Undetermined 

Cs  — *-Ba 

Undetermined 

Ba  — vLa 

Undetermined 

La  — *•  Ce 

3.3  (100^) 

Cr  — rPr 

.22  (  6%),  .52  (I25t),  .72  (5%),  1.09  (40%) 

1.38  (37%) 

Pr — -Nd 

.933  (100%) 
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144. 


54^®  ^8*^® 


Half- 

-lives  (sec): 

Mass  Diff . (Mev) : 

Cameron 

Corvell 

Xe 

1 

Xe 

—  Cs 

6.08 

5.10 

Cs 

1.5  ^ 

Cs  — 

— ^Ba 

8.45 

7.73 

Ba 

3.5 » 

Ba 

— >  La 

3.46 

2.48 

La 

15  ^ 

La 

—  Ce 

6.03 

5.11 

Ce 

2.46x10’^ 

Ce  - 

—  Pr 

.912 

-.136 

Pr 

1.04x10^ 

Pr 

—  Nd 

3.29 

2.49 

Nd 

stable 

Decay  Schemes: 

Xe  — -Cs 

Undetermined 

Cs  — -Ba 

Undetermined 

Ba  — -La 

Undetermined 

La  -  C  e 

Undetermined 

Cr  — -Pr 

.186  (20$),  .320  (72$),  2.40  (8$) 

Pr — -Nd 

.80  (1$),  2.29  (1.3$),  2.98  (97.7 
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145. 


lo^'^ 


Half-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Cs  1  » 

Cs - * 

Ba 

7.21 

6.37 

Ba  2  * 

Ba - ► 

La 

5.24 

5.07 

La  9  * 

La - - 

Ce 

4.87 

4.12 

Ce  180 

Ce - ' 

Pr 

2.70 

2.47 

Pr  2.12x10* 

Nd  stable 

Pr - » 

•Nd 

2.66 

1.52 

Decay  Schemes: 

Cs - >-Ba 

Undetermined 

Ba - >La 

Undetermined 

La - *-Ce 

Undetermined 

Ce — -Pr 

2.0  (10056) 

Pr  — »Nd 

1.8  (10056) 
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146. 


Half- 

-lives  (sec): 

Mass  Diff.(Mev): 

Cameron 

Corvell 

Cs 

1 

Cs 

— ^Ba 

9.28 

8 . 84 

Ba 

1.5  * 

Ba 

-^La 

3.99 

3.70 

La 

4  * 

La 

-^Ce 

6.65 

6.28 

Ce 

840 

Ce  — 

-^Pr 

1.55 

1.14 

Pr 

1.44x10^ 

Pr 

-^Nd 

4.45 

3.72 

Nd 

stable 

Decay  Schemes? 


Cs- 

—  Ba 

Undetermined 

Ba 

-^La 

Undetermined 

La  — »-Ce 

Undetermined 

Ce — ^Pr 

.7  (10056) 

Pr — ^Nd 

2.3  (4456),  3.7  (5656) 
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147. 


66^ - — g-^Pm — ^Sm 


Half- 

-lives  (sec); 

Mass  Diff . (Mev) ; 

Cameron 

Corvell 

Ba 

1  * 

Ba  — 

—  La 

6.06 

6.22 

La 

2  * 

La  — ►  Ce 

5. a 

5.30 

Ce 

72 

Ce 

-^Pr 

3.33 

3.68 

Pr 

720 

Pr  — 

Nd 

3.30 

2.76 

Nd 

9.59x10^ 

Nd  — 

--  Pm 

.986 

1.14 

Pm 

7.88x10’ 

Pm 

Sm 

.9a 

.223 

Sm 

stable 

Decay  Schemes: 

Ba  — »•  La 

Undetermined 

La  — «-  Ce 

Undetermined 

Ce  ^Pr 

Undetermined 

Pr  — Nd 

Undetermined 

Nd  —  Pm 

.212  (3%),  .368  (20%),  .810  (77%) 

Pm  — »■  Sm 

.225  (100%) 
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14s. 


Half- 

-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Coryell 

Ba 

1 

Ba- 

- -La 

4.69 

4.91 

La 

2  ^ 

La 

- -Ce 

7.48 

7.43 

Ce 

42 

Ce 

- -Pr 

2.09 

2.39 

Pr 

120 

Pr 

- -Nd 

5.08 

4.91 

Nd 

stable 

Nd 

- -Pm 

-.167 

-.126 

Pnij 

3.72x10® 

Pm 

- -  Sm 

2.74 

2.39 

Pnig 

4.67x10® 

Sm 

stable 

Decay  Schemes: 

Ba  — -  La 

Undetermined 

La  — -Ce 

Undetermined 

Ce  — -Pr 

Undetermined 

Pr  — -Nd 

Undetermined 

Nd  — -Pm^ 

No  Transitions  -  Shielded 

Nuclide 

Nd  — -  Pnig 

No  Transitions  -  Sheilded 

Nuclide 

Pfflj^  — Sm 

.39  (42%),  .49  (22^),  .56 

(7%),  68  (22%) 

IT  {!%) 

Pnia — *Sm 

.99  (4156),  1.9  (14^),  2.45  (45%) 
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149. 


B7La - g-gC® - 


Half-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

La 

1.5  * 

La - -  Ce 

6.11 

6.35 

Ce 

3.5  ^ 

Ce - '  Pr 

4.16 

4.75 

Pr 

30  » 

Pr - -  Nd 

3.84 

3.85 

Nd 

6 . 48x10^ 

Nd - ►Pm 

1.61 

2.25 

Pm 

1.91x10^ 

Pm - ►  Sm 

1.59 

1.35 

Sm 

stable 

Decav  Schemes: 

La- 

— ►  Ce 

Undetermined 

Ce 

^  Pr 

Undetermined 

Pr 

—  Nd 

Undetermined 

Nd 

— 'Pm 

.95  (1656),  1.1  (43^) 

,  1.5  (3156) 

Pm 

— »  Sm 

.784  (356),  1.071  (9756) 
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150. 


57 


La 


"Ts 


Ce 


-TgPr — ToNd - ffiPm — e^sSm 


Half-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

La 

1.5  ^ 

La - ^  Ce 

8.52 

8.23 

Ce 

2.5  ^ 

Ce - -  Pr 

2.79 

3.26 

Pr 

15  ^ 

Pr - -Nd 

5.92 

5.73 

Nd 

stable 

Nd - -  Pm 

.378 

.765 

Pm 

9.72x10^ 

Pm  - -Sm 

3.37 

3.23 

Sm 

stable 

Decay  Schemes: 

La 

— -Ce 

Undetermined 

Ce 

^Pr 

Undetermined 

Pr 

—  Nd 

Undetermined 

Nd 

— *■  Pm 

No  Transition  - 

Shielded  Nuclide 

Pm 

— -Sm 

.35  (.4%),  .52  ( 

.4%),  .9  (1%),  1. 

02  (1.4; 

1.36  (22%),  1.46  (8.5%),  1.8  (22%),  2.28  (29.3%) 
2.77  (2.5%),  3.10  (12.5%) 
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151. 


gy  La  ^8^®  60^^  6*1^®  63^®  63^^ 


Half-lives  ( sec) : 

Mass  Diff  (Mev) :  Cameron 

Corvell 

La 

1  » 

La - ^Ce  6.99 

7.05 

Ce 

2 

Ce - *"Pr  5*20 

5.46 

Pr 

7 

Pr - r  Nd  4-  55 

4.57 

Nd 

720 

Nd - ►Pm  2.46 

2.98 

Pm 

1.02x10^ 

Pm - ►  Sm  2.13 

2.09 

Sm 

2.84x10® 

Sm - *  Eu  .  182 

.496 

Ehi 

stable 

Decay  Schemes s 

La 

— ^Ce 

Undetermined 

Ce 

— ^Pr 

Undetermined 

Pr 

— »Nd 

Undetermined 

Nd 

-^Pm 

.24  (1%),  .54  (3%),  1-19  (30%), 

1.65  (7%) 

1.82  (1556),  2.06  (45^) 

Pm 

— ►  Sm 

.374  (5%),  .454  (3%),  .75  (6%), 

.85  (39%) 

.871  (2%),  .968  (1%),  1.03  (10^) 

,  1.09  (6 

1.125  (856),  1.129  (5/6),  1.20  (11%) 

Sm 

— »  Eu 

.0543  (1.6%),  .0759  (98.4%) 

133 


GNE/Phys/65-3 


152. 


68^® - - 6lPf“ — s^Sm 


Half-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

Ce  1.5  ^ 

Ce - ^  Pr 

3.68 

3.99 

Pr  5 

Pr - -  Nd 

6.96 

6.43 

Nd  180  ^ 

Nd - *-  Pm 

1.09 

1.51 

Pm  360 

Sm  stable 

Pm - Sm 

4.21 

3.95 

Decay  Schemes: 

Ce — *.Pr 

Undetermined 

Pr — »  Nd 

Undetermined 

Nd — >Pm 

Undetermined 

Pm — »Sm 

2.2  (10056) 
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153. 


ss 


Ce- 


69 


Pr- 


60 


Nd- 


6  1 


Pm - ^gSm  — 


63 


Eu 


Half-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

Ce  1.5  ^ 

Ce - -Pr 

6.22 

6.15 

Pr  3 

Pr - *  Nd 

5.44 

5.27 

Nd  18  ^ 

Nd - -  Pm 

3.50 

3.69 

Pm  330 

Pm - >  Sm 

2.85 

2.81 

Sm  1.68x10® 

Sm - -  Eu 

1.03 

1.23 

Eu  stable 

Decay  Schemes: 

Ce - ►Pr 

Un(determined 

Pr — -Nd 

Undetermined 

Nd  — -Pm 

Undetermined 

Pm - *Sm 

1.65  (10056) 

Sm - -Eu 

.641  (33.456),  .696  (45.6%),  .707  [.5%) 

.721  (.356),  .803 

(2056) 
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154. 


ggCe - ^gPr  — g^oNd - g-^Pm  — gigSm 


Half- 

-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

Ce 

1  ^ 

Ce - r 

Pr 

4.66 

4.79 

Pr 

2  * 

Pr - - 

Nd 

7.98 

7.17 

Nd 

12 

Nd - ^ 

Pm 

1.98 

2.35 

Pm 

150 

Pm - ^ 

Sm 

5.26 

4.73 

Sm 

stable 

Decay  Schemes: 

Ce — -Pr 

Undetermined 

Pr  _  Nd 

Undetermined 

Nd  — »■  Pm 

Undetermined 

Pm  — »•  Sm 

2.5  (100%) 
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155. 


69 


Pr - ^oNd g^iPm  — ^gSm g^Eu 64 


Half- 

-lives  (sec): 

Mass 

Diff  (Mev) : 

Cameron 

Corvell 

Pr 

2  ^ 

Pr- 

- -Nd 

6.43 

6.16 

Nd 

5  ^ 

Nd 

- -Pm 

4.53 

4.60 

Pm 

60  * 

Pm 

- Sm 

3.74 

3.74 

Sm 

1.32x10^ 

Sm- 

- »  Eu 

2.08 

2.18 

Eu 

5.36x10’^ 

Eu 

- -  Gd 

.882 

1.17 

Gd 

stable 

Decay  Schemes; 


Pr- 
Nd 
Pm 
Sm 
Eu  ■ 


^Nd 

'Pm 

Sm 

Eu 

■Ge 


Undetermined 
Undetermined 
Undetermined 
1.39  (5%),  1.53  (95$) 

.140  (43$),  .160  (37$),  .185  (10$) 
.247  (15$) 
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156. 


59 


Pr 


60 


Nd 


61 


^Pm- 


^  Sin¬ 


es 


j^Eu- 


Ti 


Gd 


Half-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

Pr  1.5 

Pr - ,  Nd 

9.07 

7.94 

Nd  4.5  * 

Nd - »  Pm 

2.97 

3.19 

Pm  25  ^ 

Pm - ^  Sm 

6.29 

5.52 

Sm  3.24x10* 

Sm - ^  Eu 

.566 

.771 

Eu  1.3x10® 

Eu - ^  Gd 

3.30 

3.10 

Gd  stable 

Decay  Schemes*. 

Pr - »Nd 

Undetermined 

Nd - ►Pm 

Undetermined 

Pm — ^Sm 

Undetermined 

Sm  — >Eu 

.43  {55%),  .715  (45%) 

Eu  — ►  Gd 

.3  (23%),  .48  (34%), 

1.21  (11%), 

2.447  ( 
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157. 


59^^  ^  Eu  6’4^^ 


Half- 

-lives  ( sec) : 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

Pr 

1.5  * 

Pr  — 

-^Nd 

6.71 

Nd 

3  * 

Nd 

— >  Pm 

5.61 

5.16 

Pm 

13  * 

Pm 

— »  Sm 

4.74 

4.31 

Sm 

30 

Sm  — 

— f  Eu 

3.11 

2.76 

Eu 

5 . 40x10'* 

Eu 

—  Ge 

1.79 

1.76 

Gd 

stable 

Decay  Schemes: 


Pr  — '■Nd 

Undetermined 

Nd  — '-Pm 

Undetermined 

Pm  — '-Sm 

Undetermined 

Sm — '■Eu 

Undetermined 

Eu  — -^-Gd 

1.27  {100%) 
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158. 


60 


Nd- 


61 


Pm— 


62 


Sm- 


63 


Eu” 


64 


Gd 


Half- 

-lives  (sec): 

Mass 

Diff  (Mev) : 

Cameron 

Coryell 

Nd 

2 

Nd 

- =•  Pm 

4.03 

3.78 

Pm 

6  ^ 

Pm 

- >  Sm 

7.38 

6.06 

Sm 

900 

Sm- 

— Eu 

1.56 

1.38 

Eu 

3.6x10^ 

Eu 

- ^  Gd 

4.33 

3.66 

Gd 

stable 

Decay  Schemes: 

Nd - ►Pm 

Undetermined 

Pm — ►  Sm 

Undetermined 

Sm — -►  Eu 

Undetermined 

Eu - ►  Gd 

2.65  (10056) 

140 


GNE/Fhys/65-3 


159. 


60 


Nd 6^  Sm ga  ^4^ - 65 


Tb 


Half-lives  (sec): 

Mass  Diff  (Mev) : 

Cameron 

Corvell 

Nd 

1.8  ^ 

Nd - »-Pm 

6.68 

5.76 

Pm 

3.5 

Pm - >  Sm 

5.80 

4.91 

Sm 

40  * 

Sm - *•  Eu 

4.21 

3.36 

Eu 

1 . 14x10^ 

Eu - -Gd 

2.79 

2.36 

Gd 

6 . 48x10^ 

Gd - *Tb 

1.46 

.960 

Tb 

stable 

Decay  Schemes: 

Nd — . 

Pm 

Undetermined 

Pm - ^ 

'Sm 

Undetermined 

Sm — ' 

>Eu 

Undetermined 

Eu - • 

.Gd 

2.2  (100$) 

Gd - - 

.Tb 

.59  (13%),  .89  {21S) 

,  .947  (63%) 
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Appendix  C 
Computer  Programs 

This  appendix  contains  the  information  necessary  to  understand 
the  operations  of  the  four  computer  programs  referred  to  in  the 
body  of  the  thesis.  This  information  consists  of  a  listing  of  the 
program' s  Fortran  statements,  a  list  of  the  definitions  of  all  terms 
used,  and  a  flow  chart.  The  information  for  each  of  the  four  pro¬ 
grams  is  presented  in  this  order.  In  programs  two  and  four,  the 
omission  of  the  definition  of  a  variable  in  a  function  or  subroutine 
implies  that  the  definition  found  earlier  still  applies  and  has  been 


omitted  for  brevity. 


GNE/Phys/65-3 


C  PROGRAM  I 

C  PKOGKAM  ro  COMPUTE  AND  CHECK  PZ  VALUES 

dimension  CHAIN! 85) , AI 11 ) , NOK! U ) , ANO! 85 ) , Y ! 85 ) , P ( 8 5 ) 

READ  16,ITF,YM 

N=l 

1  GU  TO! 2, 3, 8,20) , I TF 

2  READ  1 7, ANC! M) ,NI, ZMINtZP, WtCHAININ) , WW, ( 1 1 ) 

READ  18,WlvW 

ZP  =  ZP<-.  316 
GO  TO  5 

20  READ  17,AN0!N) ,NI ,ZMIN,ZP, W,WW,WWW,( in 

READ  18,CHAIN!N) 

ZP=ZP+.023 
GO  TO  5 

3  READ  17,ANG!N) ,NI,ZMIN,ZP,CHAIN!N),W,WW, ( 11) 

READ  18, WWW 

GO  TO  5 

4  READ  I  7, ANO! N) ,N I , ZMIN,ZP, W ,WW,CHA IN! N ) , (  1 1  ) 

READ  18, WWW 

ZP=ZP-.475 

5  READ  52, (NOK! I ) , I=l,NI ) 

Z=ZMIN 

X=CHA IN! N) 

DU  9  I=l,NI 
Z=Z+1. 

IF!1-N0K! I ) )6,7,8 

6  Z=Z-l. 

7  C=ABS!Z-ZP) 

IF!2.2-C)22,23,23 

23  Y  I  £00=6.84  4462  I* X*eXP!-.,  7*C<‘C)*YM 
GO  TO  24 

22  YIELD  =  2.6  7  5E27*X<'EXP!-7. 38*C)*YM 

24  A!  I  ) =YIEL0/2. 

GO  TO  9 

8  C=A8S!Z-ZP) 

IF!2.2-C)27,28,28 

28  YIEL0  =  6.84  44E21*X<'6XP  !-.7*C*C)*YM 
GO  TO  29 

27  Y  ieLD  =  2.6  7  5E2  7«^X*EXP!-7.  TB^O’f'YM 

29  A!I)=YIFLD 

9  CONTINUE 

PUNCH  100,ANO!N) ,ZMIN,NI 
PUNCH  lOl , ! A! I  )  , I  =  l  ,NI ) 

Y!N) =0.0 
DO  10  1=1, NI 

10  Y!N)=Y!N)+A! I )/1.45E22 
N  =  N+l 

GO  TO  1 

11  DO  12  N=1,H3 

12  P!N)=!CHAIN!N)-Y!N))<=100./CHAIN!N) 

PRINT  103 
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PKINl  10'-. 
no  13  ,\i=l,4  1 
M  =  + 

13  PK  1  M  I  0  ‘j  ,  AND  (N)  ,CHAtNJ(N),Y(N),P(\'),  ANP  (  M )  ,  C  H  A  I N  (  M  )  ,  Y  (  M  )  ,  P  (  M  ) 

PK  INT  106,  AN(1(  )  tCIlAI  N(  4?  )  ,  Y('.?  )  ,  P(  4?  ) 

STOP 

16  FOkPAK  I^,P8.4) 

17  FURMAKFIO.O,  I  10,5F10.4) 

18  FORMAT (F 10.4) 

52  FnPMAT(llI2) 

100  FfJRMAT(  2F10.0,  I  10) 

ICl  FURMAT(  1P6F12.4) 

IC'3  format  (  IHl  ,'39X  ,  52HCnMPAR  1  SON  OF  CHAIN  YIFLO  AND  SUM  OF  ISOTOPIC  YI 
IFLOS) 

104  FORMAT!  IH0,2( 8X, 5HA  NO , , 4X , 1 1 HCH A  1 N  Y I  FLO , 3X , 1 4H I SOTOP I C  YIELD, 2X, 
laUPCT  DIFF ,5X) , IH  ) 

105  FnRMAr(2(5X,F10.0,2F15.4,F10.2,5X)  ) 

1G6  F0KMAT(5X,F10.0,2F15-.4,F10.2) 

ENO 
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A,  ANO 
C 

CHAIN 

ITF 


NI 

NOK 


P 


Y 

YIELD 

YM 

Z 

ZMIN 


ZP 

X,  W,  WW, 
WWW 

M,  N,  I 


PROGRAM  I 

Atomic  weight  of  mass  chain 
Absolute  value  of  Z  -  Zp 
Chain  yield 

Indicator  of  the  type  of  fission;  1  indicates  Pu-239 
thermal  neutron  fission;  2  indicates  U-235  thermal 
neutron  fission;  3  indicates  thermonuclear  neutron 
fission;  4  indicates  U-235  fission  spectrum  neutron 
fission 

Number  of  isotopes  in  mass  chain 

Indicator  of  isomeric  split:  0  indicates  no  isomeric 

pair;  1  indicates  first  member  of  isomeric  pair;  2 

indicates  second  member  of  isomeric  pair 

Percent  difference  between  chain  yield  and  sum  of 

isotopic  yield 

Sum  of  isotopic  yields 

Isotopic  yield 

Number  of  megatons  of  fission 
Atomic  number  of  isotope 

One  unit  less  than  the  atomic  number  of  the  first  isotope 
in  the  chain 

Most  probable  charge  formed  in  fission 

Dummy  variables 
Indicies  of  "DO"  loops 
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PROGRAM  I 


Fi-.  17 
Flov:  Cliart, 


for  Pro"ram  I 
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PROGRAM  I  I 

THIS  PROGRAM  COMPUTES  AND  WRITES  ON  TAPE  THE  NORMALIZED  SPECTRAL 
SHAPES 

COMMON  A,Z,StSO,LtTHALF,GIRLY,FRAC,EOP,FOG 
READ  91, MF, ALPHA 
READ  17,AN0,NI,ZMIN,(16) 

WRITE  TAPE  3,ANO,NI 
KA=0 
A  =  ANO 
Z=ZMIN 

DO  15  L=1,NI 
READ  18, IO,THALF,JO 
IF  (KA)  5,4,5 
Z=Z+1. 

IF  (JO)  6,7,6 
KA  =  1 
GO  TO  8 
KA=0 

IF(THALF)9,30,9 

30  MI=1 

FUNN=0.0 

WRITE  TAPE  3,MI,THALF 
WRITE  TAPE  3,FUNN 
GO  TO  15 

9  SO=(Z/137.0371)**2 

S=SQRTFC 1.0-S0)-1.0 
I I=IABS( 10) 

00  15  L0LA=1,II 
IF( 10)20,20,21 

20  IF(MF)22,22,23 

22  READ  90,DIFF,W 
GO  TO  24 

23  READ  90,W,DIFF 

24  EOP=OIFF*ALPHA 
FRAC  =  1.000 
GO  TO  14 

21  READ  19,E0P,FRAC 

14  GIRLY=SHAPE( 2) 

CALL  BINERILOLA, II ) 

15  CONTINUE 
GO  TO  1 

16  END  FILE  3 
REWIND  3 
STOP 

17  FORMAT(F10.0, I5,F10.0) 

18  FORMAT! 15, 1PE10.2, 15) 

19  FORMAT(2F10.3) 

90  FORMAT(2E10.3) 

91  FORMAT( I5,E10.4) 

END 
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FUNCTION  SHAPE  (R) 

♦  THIS  SUBROUTINE  COMPUTES  LOG  10  FT  VALUES  FOR  THE  FISSION  PRO-  ♦ 
♦DUCTS.  THESE  LOG  10  FT  VALUES  ARE  IN  TURN  USED  IN  DETERMINING  ♦ 
♦THE  FORBIDENESS  SHAPE  FACTORS  TO  BE  USED  WITH  THE  FERMI  BETA  OE-^ 
♦CAY  THEORY  IN  FORMULATING  SPECTRAL  SHAPES.  ♦ 

COMMON  A,Z,S,SC,L,THALF,GIRLY,FRAC,EOP,FOG 
FP=GORTI 1) 

FP I =ALOG I FP^THALF) 72.30258 
PSQ=5.46AE03/IFP^THALF ) 

IF  I6.-FP1)  2,1,1 

1  SHAPE=20. 

RETURN 

2  IF  I8.-FP1)  4,3,3 

3  SHAPE=-1. 

RETURN 

4  IF  (lO.-FPl)  6,5,5 

5  SHAPE=0.0 
RETURN 

6  IF  I15.5-FP1)  8,7,7 

7  SHAPE=1.0 
RETURN 

8  SHAPE=2.0 
RETURN 
END 
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FUNCTICN  GGRT  (BO) 

*  THIS  SUBROUTINE  COMRUTFS  THE  FERMI  INTEGRAL  FUNCTION  USING  THF  ♦ 
«T«ARLZOIl)AL  RULE.  THt  INTEGRATION  IS  CALCULATED  IN  TWO  PARTS.  ♦ 
■“FOk  T  LESS  THAN  OK  EUUAL  TO  l.O  MEV  EACH  0.1  mev  INTERVAL  IS  ♦ 

♦SEPERATEU  INTO  TEN  INCREMENTS.  FOR  T  GREATER  THAN  l.O  MEV  EACH  ♦ 
*1.0  MEV  INTERVAL  IS  SEPFRATEO  INTO  TEN  INCREMENTS.  ♦ 

CU»-’MON  A,Z  ,  S,SO,L  ,  T  HAL  F  ,  G  I  RL  Y  ,  FR  AC  ,  E  CP  ,  FOG 
GnftT=0. 

E=0. 

Cl AMY= .01957 
E0  =  ( (EOP/.Sl  !)♦!.) 

NnT=lOl 
L  nG=PHi ( 1 ) 

IF  (l.-FOP)  2,2,1 

1  NOT=(E0P/.01+l. ) 

2  no  6  K=1,NUT 
E=E+CL AMY 

IF  (K-l)  T,3,A 

3  E=1.0 

A  X=SORTF( E*E-1.0) 

Y  =  .0^»5H5*Z*E/X 
Z AM=En-E 

f  =Y<=X*F  JG*ZAW*ZAM*E*CLAMY*(  (  S0*F*E  +  X*X/4.  )**S)/ABS(  l.-EXPF(-Y)  ) 

IF  (NOT-K)  5,5,6 

5  F=F/2. 

6  GORT=GaKTfr 

IE  (  l.-EOP)  8,  T,7 

7  RETURN 

8  NOT  =  { FOP-l . ) /.  I 
CLAMY=, 195 ( 

DO  13  K=1,NUT 
IF  ( K-  1  )  9,9,10 

9  E=2. 95 7 
GO  TO  11 

10  E=E+CLAMY 

11  X=SOKTF (E*E-1.0) 

Y=.045d5*Z*F/X 
ZAM  =  EIJ-E 

F  =  Y*X*FnG*ZAMtZAM*E*CLAMY* ( ( Sn*E*E+X*X/4 . ) **S ) /ABS ( l.-EXPF (-Y )  ) 

IE  (K-l)  12,12,13 

12  F=E/2. 

13  GURT=GnRT+F 
RETURN 

END 
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FUNCTION  PHI  (BB) 

♦  THIS  SUBROUTINE  COMPUTES  PHKZl,  A  FUNCTION  OF  ATOMIC  CHARGE  AND* 

♦  MASS  NUMBER  ONLY.  NOTE  THAT  PHKZ)  IS  UNITLESS.  ♦ 

COMMON  A, Zf S»SOtL  f THALF,GIRLY,FRAC  f EOPfFOG 

C=5. 1794E12 

0=A^^0. 33333333 

R=I. 11 2E-15^(0-. 62025/D) 

R  S  ^  S 

PHI=( (C*R)^^RS)*( (2./GAMMA(3.+RS) )**2)*( 1.+S/2.) 

RETURN 

END 
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FUNCTION  GAMMA  (X) 

♦THIS  FUNCTION  COMPUTES  THE  GAMMA  FUNCTION  OF  X  FOR  POSITIVE  ANO^ 
♦NEGATIVE  VALUES  OF  X.  ♦ 

COMMON  A,Z,S,SO,L,THALF,GIRLY,FRACf EOP.FOG 
IF  (X)  8,2,3 
PRINT  10, X 
GAMMA=l.70141l79E38 
RETURN 

Y  =  X 

GLN  =  0.918938568+(Y-0.5)^ALOG(Y)-Y+(  1.0/1  12.0^Yn-f  1.0/ ( 360.0^Y^^3 » 
l)+( 1.0/( 1260.0^Y^^5) )-(1.0/( 1680.0^Y^^7) )♦( 1 . 0/( 1 1 88.0^Y^^9 I ) 

IF  (GLN-88.0)  5,2,2 

5  G=EXPF(GLN) 

IF  (X)  7,1,6 

6  GAMMA=G 
RETURN 

7  GAMMA=3.1415927/(SINFI3.1415927^Y)^GI 
RETURN 

8  K=-X 

Y  =  K 

IF  (Y+X)  9,2,9 

9  Y=1.0-X 
GO  TO  4 

10  FORMAT  (18H  GAMMA  FUNCTION  OF , IPE 1 5. 7, 27H  EXCEEDS  ♦/- I .70141 179E+3 
18.  I 

END 
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SUBKOUTINL  H I N ER ( L C L A ,  I  I  ) 

C 

C  *  THIS  SUCRUUriNE  CALCULATES  AND  SUMS  ALL  OF  THE  INDIVIDUAL  IN-  * 

C  <=STANT  ANEOUS  AND  INTEGRATED  SPECTRA  FROM  THF  FISSION  PRODUCTS.  « 

C  *ir  combines  all  of  the  activity,  half-life,  EPF,  AND  FORDIDF-  ♦ 

C  ^NFSS  INFCRMATIfiN  IN  CALCULATING  THE  EFFECTIVE  BETA  SPFCTRUM  OF  ♦ 

C  ■('THE  FISSION  PRODUCTS.  * 

C 

COVMON  A,Z,S,Sn,L, THALF ,GIRLY,FRAC ,EGP,FnG 
DIMENSION  FUN(-?10)  ,FUNM210,  ID 

4  SUM^O.O 
E=0,0 
N0T=10l 
CLAMY=. 01957 

E0  =  (  (E0P/.5ll)«-l.  ) 

IE  (l.-EOP)  6,6,5 

5  NnT  =  EOP/.OUl. 

6  DO  18  K=1,N0T 
E=E+CLAMY 

IF  (K-1)  7,7,8 

7  E=1.0 

8  X=S0RTF(E*E-1.0) 

Y  =  .04585«'Z*E/X 
7AM=EC-E 

F  =  Y*X*F0G*Z  AM'C'ZAM«‘E«CL  AMY#  I  I  SO#E  #E  ♦X  ♦X /4  .  )##S)/ABS(1.-FXPF(-Y)) 

IF  (GIRLY)  10,11,9 

9  IF  (GIRLY-2.0)  12,13,14 

10  FUN(K)=F#( X##2*ZAM##2) 

GO  TO  15 

11  FUNIK  I  =F#{  X##2  +  ZAM##2) /12. 

GO  TO  15 

12  FUN( K) =F#( X##4+10. ♦( X*#2 )# ( Z AM#*2 ) /3.+Z AM##4 ) 

GC  TO  15 

13  FUN(K)  =F#(  X##6<-7.#X#X#ZAM#ZAM#{X#X«-ZAM#ZAMD-ZAM#*6) 

GO  TO  15 

14  FUN(K)=F 

15  IF  (NOT-K)  16,16,17 

16  SUM=SUM+FUN(K)/2. 

GO  TO  18 

17  SUM=SUM+EUN(K) 

18  FUN(K) =FUN(K)/CLAMY 
IF  (l.-EOP)  19,33,33 

19  M0T=(F0P-1.)/.U1. 

CLAMY=. 1957 

DO  32  K=l,MnT 

M0T0=K+100 

IF  (K-l)  20,20,21 

20  F=2.957 
GO  TO  22 

21  £=E+CLAMY 

22  X=SQKTF(E#E-1.0) 
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Y=.04585*Z*E/X 

ZAM=EO-E 

F  =  Y*X*F0G*ZAM*ZAM*E*CLAMY*( ( S0*E*E  +  X*X/4. )**SI/ABS( l.-EXPF(-Y)  I 
IF  (GIRLY)  24,25,23 

23  IF  (GIRLY-2.0)  26,27,28 

24  FUN(MOTO)=F*(X**2+ZAM**2) 

GO  TO  29 

25  FUN(MOTO)=F*(  X**2<-ZAM*«2)/12. 

GO  TO  29 

26  FUN(M0T0)=F«(X**4+10.*(X**2)*( ZAM**2 )/3.+ZAM**4» 

GO  TO  29 

27  FUN( MOTO)=F*( X**6+7.*( X**2 )*(ZAM**2>*( X**2+Z AM**2 ) +2 AM**6 ) 

GO  TO  29 

28  FUN(MOTO)=F 

29  IF  (K-1)  30,30,31 

30  SUM=SUM+FUN(M0T0I/2. 

GO  TO  32 

31  SUM=SUM+FUN( MOTO) 

32  FUN(MOTO)=FUN(MOTO)/CLAMY 
MI=MOTO 

GO  TO  34 

33  MI=NOT 

34  IF(L0LA-1)90,77,90 

77  DO  79  K0=l,210 

79  FUN1(KO,L)=0.0 

90  00  83  K0=1,MI 

83  FUNI(K0,L)=FUN1(K0,L)+FUN(K0)*FR AC/SUM 

84  IF(LOLA-II >75,78,75 

78  WRITE  TAPE  3,MI,THALF 

WRITE  TAPE  3,(FUN1(K0,L),K0=1,MI) 

75  RETURN 

END 
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A,ANO 

ALPHA 

DIFF 

EOP  = 

FHAC  = 

FUNN  = 

10,  II 
JO 

MI 

MF 

NI 

SHAPE 

THALF  = 

W 

Z 

ZMIN 

FOG,  GIRLY, 
S,  SO, 

KA,  L 


PROGRAM  II 

Main  Program 

Atomic  weight  of  mass  chain 

Fraction  of  transition  energy  going  by  beta  decay 
Mass  difference  between  two  isotopes 
End  point  energy 

Fraction  of  decays  with  specific  end  point  energy 
Normalized  spectral  value 
Number  of  different  betas 

Branching  indicator:  0  indicates  no  branching; 

1  indicates  branching 
Number  of  energy  increments 

Mass  Formula  indicator:  0  indicates  Coryell’s  formula; 

1  indicates  Cameron’ s  formula 

Number  of  isotopes  in  mass  chain 

See  Function  SHAPE 

Half-life  of  isotope 

Dummy  variable 

Atomic  number  of  isotope 

One  unit  less  than  the  atomic  number  of  the  first  isotope 
in  the  chain 

Interim  Calculation  Variables 
Indicies  of  ”D0”  loops 
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SHAPE 

FP 

GORT 

FPI,  PSQ 
R 

GORT 

BO 

CLAMY 
E,  EO,  F, 
Y,  ZAM 
K,  NOT 

PHI 

BB 

C,D,R,RS 

GAMMA 

X 

G,GLN,K,Y 


Function  SHAPE 

A  fvinction  which  calculates  the  "log  ft"  values 

Fermi  Integral  Function 

See  Function  GORT 

Interim  Calculation  Variables 

Dummy  variable 


Function  GORT 

A  function  which  calculates  the  Fermi  Integral  Function 
Dummy  variable 
Energy  Interval 

Interim  Calculation  Variables 
Indicies  of  "DO"  loops 

Function  PHI 

A  fvinction  which  computes  a  term  used  in  GORT 
Dummy  variable 

Interim  Calculation  Variables 

Function  GAMMA 

A  function  which  computes  the  value  of  the  gamma  function 
Dummy  variable 

Interim  Calculation  Variables 
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BINER 

FUN 
FUN  1 
SUM 

NOT,  MOT, 
MOTO 
LOLA,  K, 
KO 


Subroutine  BINER 

A  subroutine  which  computes  and  stores  the  values 
of  the  normalized  spectral  shapes 
Value  of  the  spectral  shape 
Normalized  value  of  the  spectral,  shape 
Integral  value  of  the  spectral  shape 

Interim  Calculation  Variables 

Indicies  of  "DO"  loops 
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PROGRAM  II 


_  lov/  Cliart  Tor  Proc:rar;i  II 
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Fic:.  16  (Continued) 


15S 
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PROGRAM  III 

PROGRAM  TO  GENERATE  NEW  INPUT  DATA  TAPE 
DIMENSION  A( 11) ,THALF( 11) tBI 11) ,FUNl(2I0,ll) ,MI(  11 ) 
NX=0 

READ  101«ANO,ZM1N,NI ,(8) 

READ  102»(A( I),I=1»NI) 

READ  TAPE  1,AN0,NI 
00  7  I=1,NI 

READ  TAPE  1 , MI ( I ) ,THALF ( I ) 

MM=MI ( I ) 

READ  TAPE  1 , ( FUNl I K, I ) , K=1 , MM) 

DO  4  1=1, NI 
IF(THALF( I ))2,3,2 
B(  I  )=0.693/THALF( I ) 

GO  TO  4 
B(  I  )=0.0 
CONTINUE 

WRITE  TAPE  3yAN0,NI,(A(I),I=l,NI),(B(I),I=l,NI) 

DO  5  L=1,NI 
WRITE  TAPE  3,MnL) 

MM=MI (L) 

5  WRITE  TAPE  3  ,  ( FUNl (K ,L ) , K= 1 , MM ) 

GO  TO  1 

8  END  FILE  3 

NX=NX+1 
PRINT  103, NX 
REWIND  1 
IF(NX-4) 1,6,6 

6  REWIND  3 
STOP 

101  F0RMAT(2F10. 0,110) 

102  FORMAT! 1P6E12. 4) 

103  F0RMAT(I5) 

END 
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PROGRAM  III 

A 

=  Concentration  of  isotope  immediately  after  fission 

ANO 

=  Atomic  weight  of  mass  chain 

B 

=  Decay  constant  of  isotope 

FUN  1 

=  Normalized  spectral  shape  value 

MI 

=  Number  of  energy  increments 

NI 

=  Niomber  of  isotopes  in  mass  chain 

THALF 

=  Half-life  of  isotope 

ZMIN 

=  One  unit  less  than  the  atomic  number  of  the  first 

isotope  in  the  chain 

MM,  NX 

=  Interim  Calculation  Variables 

I,  K,  L 

=  Indicies  of  "DO"  loops 
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PROGRAM  III 


Do  4 
UI.NI 


■  .  1 

Calculate 

B 


Flc.  19 

Flov.'  Chart  for  ProGrtam  III 
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C  PROGRAM  IV 

C  THIS  PROGRAM  SOLVES  THE  DECAY  EQUATIONS  EXACTLY 

C 

DIMENSION  A( 11 ) ,AA( 11) ,B( 11) ,BB( 1 1 ) , FUNl ( 2 10, 1 1 ) , MM( 1 1 ) , BC ( 1 1 ) ,L ( 1 
11)  ,AB( 11 ,20) ,DC( 10,10) ,E( 10) 

COMMON  EVE (20) , SUN (20, 2 10) , ITF,SUM( 2 10) 

READ  990,ITF,YM 

ITT=ITF-1 

IF( ITT) 12,  13,12 

12  DO  lA  NX=1,ITT 

lA  CALL  WINDER 

13  READ  990, NT 

READ  991, (EVE(K) ,K=1,NT) 

DO  97  KK=1,83 

READ  TAPE  1 , ANO, NI , ( AA ( I ) , I  =  1 , NI ) , ( BB( I ) , 1  =  1 ,N I ) ,  (  2000 ) 

DO  98  LX=1,NI 
AA(LX)=AA(LX)^YM 
READ  TAPE  1,MI 
MM(LX)=MI 

98  READ  TAPE  1 , ( FUN  1 ( KO,LX ) ,K0= 1 , MI ) 

DO  99  IX=1,NI 

READ  992, ANO, NP, (2000) 

DO  99  JJ*1,NP 

READ  993,NC,(8C(K) ,L(K) ,K=1,NC) 

DO  99  11=1, NT 
T  =  EVE( 1 1 ) 

AB(  IX, I  I )=0.0 

GO  TO  ( 10,20,30, AO, 50, 60, 70, 80, 90, 110) ,NC 

10  DO  11  K=1,NC 
LL=L(K) 

B(K)=BB(LL) 

11  A(K)=AA(LL) 

P1=BC( 1) 

1  XNl  =  P1*A( 1 )*EXP(-B(1 )*T) 

AB( IX, I  I  )=AB(  IX, II )^XN1 

GO  TO  99 

20  DO  21  K=1,NC 
LL=L(K) 

B(K)=BB(LL) 

21  A(K)=AA(LL) 

P1=BC( 1) 

P2=BC(2) 

2  DO  201  I  =  1,2 
DO  201  J  =  1,2 

201  0C(I,J)  =  B(I)-B(J) 

DO  202  I  =  1,2 

202  E(  I )  =  EXP(-B(  I  X'T) 

T1  =  E(1)/0C(2,1) 

T2  =  E(2)/0C( 1,2) 

XN2  =  P1*B( 1 )*A' 1)*(T1+T2)+P2*A(2)*E(2) 

AB(IX,II)=AB(I.(,II)+XN2 
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GO  TO  99 

30  DO  31  K=l,NC 

LL=L(K) 

B(K)=8B(LL) 

31  A(K)=AA{LL) 

Pl=BCi 1) 

P2=BC(2) 

P3=BC(3) 

3  DO  301  I  =  1,3 

DO  301  J  =  1,3 

301  OC(I,J)  =  8(I)-B(J) 

DO  302  I  =  1,3 

302  E(  I )  =  EXP(-B( I )*T) 

T1  =  E(1)/(DC(2,1)*DC(3,1» ) 

T2  =  E(2)/(0C(1,2)*0C(3,2) ) 

T3  =  E(3)/(0C(1,3)*0C(2,3) ) 

XN3  =  P1*B{1)*B(2)*A(1)*(T1+T2+T3) 

T1  =  E<2)/0C(3,2) 

T2  =  E(3)/DC(2,3) 

XN3  =  XN3+P2*B(2)*A(2)*(T1+T2)+P3*A( 3)*E(3) 
AB( [X,I I )=AB( IX, II )^XN3 
GO  TO  99 

40  DO  41  K=1,NC 

LL=L(K) 

B(K)=BB(LL) 

41  A(K)=AA(LL) 

P1=BC( 1) 

P2=BC(2) 

P3=BC(3) 

P4=BC(4) 

4  DO  401  I  =  1,4 

DO  401  J  =  1,4 

401  OC(I,J)  =  B(n-B(J) 

DO  402  I  =  1,4 

402  Ell)  =  EXP(-B(I)*T) 

T1  =  E( l)/(DC(2,l)*0C(3,l)*0C(4,l) ) 

T2  =  E{2)/(0C(1,2)*0C(3,2)*0C(4,2)) 

T3  =  E(3)/I0C(1,3)*DC(2,3)*DC(4,3)) 

T4  =  E(4)/(0C( 1,4)*DC(2,4)*0C13,4) ) 

XN4  =  P1*B(1)*B(2)*BI3)*A(1)»(T1+T2+T3+T4) 
T1  =  E(2)/(0C( 3,2)*DC<4,2) ) 

T2  =  E(3)/(0C(2,3)*DCI4,3) ) 

T3  =  E(4)/(DC(2,4)*DC(3,4) ) 

XN4  =  XN4+P2*B12)*B(3)*A(2)*(T1+T2+T3) 

T1  =  E(3)/DC(4,3) 

T2  =  EI4)/0C(3,4) 

XN4  =  XN4+P3*B(3)*A(3)*(T1+T2)+P4*A(4)*EI4) 

AB(IX,II)=AB(IX,II)+XN4 

GO  TO  99 

50  DO  51  K=1,NC 
LL=L(K) 
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B(K)=BB(LL  ) 

51  A(K)=AA(LL) 

P1  =  BC(  I ) 

P2=BC(2) 

P3=BC(3) 

P4=BC(4) 

P5=BC(5) 

5  00  501  I  =  lf5 
DO  501  J  =  lt5 

501  0C(  I  f J)  =  8(1 )-8( J) 

DO  502  I  =  1.5 

502  E(  I  )  =  EXP(-B( I )*T) 

T1  =  F(1)/(0C(2,1)*DC(3,1»*DC(4,1)*0C(5,1)) 

T2  =  E(2)/(0C(l,2)*0C(3,2)*0C(4,2)*DC(5,2)) 

T3  =  E(3)/(0C( 1,3)*0C(2.3)*0C(4,3)*DC(5,3) » 

T4  =  E(4)/(0C(  1,4)*DC(2,4)*0C(3,4)*0C(5,4) ) 

T5  =  E(5)/(DC(1,5)*0C(2,5)*0C(3,5)*0C(4,5) » 

XN5  =  Pl»B(l)«B(2)*8(3)»8(4J*A(l)*(Tl+T2+T3+T4+T5) 
T1  =  E(2)/(DC(3.2)*0C(4,2)*0C(5,2) ) 

T2  =  E(3)/(0C(2,3)*0C(4,3)*0C(5,3) ) 

T3  =  E(4)/(0C(2,4)*0C(3,4)*0C(5.4n 
T4  =  E(5)/(0C(2,5)*0C(3.5)*0C(4,5) 1 
XN5  =  XN5  +  P2*B(2)*B(3)*B(4)*A(2)*CTH-T2+T3+T4) 

T1  =  E(3)/(0C(4,3)*0C(5,3) ) 

T2  =  E(4)/(0C(3,4)*0C(5,4) ) 

T3  =  E(5)/(DC(3,5)*0C(4,5» ) 

XN5  =  XN5+P3*B(3)*B(4)*A(4) 

T1  =  E(4)/0C(5,4) 

T2  =  E(5)/0C(4,5) 

XN5  =  XN5+P4*B(4)*A(4)*(T1+T2)+P5*A(5)*E(5» 

AB(  IX, II  )=AB(  IX, II  )-»-XN5 

GO  TO  99 

60  DO  61  K=1,NC 

LL=L(K) 

B(K) =8B(LL  » 

61  A(K)=AA(LL) 

P1=BC{ 1) 

P2=BC(2) 

P3=BC(3) 

P4=BC(4) 

P5=BC( 5) 

P6=BC(6) 

6  DO  601  I  =  1,6 

DO  601  J  =  1,6 

601  DC(I,J)  =  B(I)-B(J) 

DO  602  I  *  1,6 

602  E( I )  =  EXP(-B( I )*T) 

T1  =  E{l)/(OC(2,l)*DC(3,ll*OC(4,l)*OC(5,l»*OC<6,l) ) 
T2  -  E(2)/(0C( 1,2)*0C(3,2)*0C(4,2)*0C(5,2)*DC{6,2) ) 
T3  =  E(3)/(0C(’ ,3»*0C(2,3)*DC(4,3)*0C(5,3)*0C(6,3n 
T4  =  E(4)/(0C( l,4)*DC(2f4)*DC(3,4)*DC(5,4)*0C(6,4n 
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=  t  (  s  )  /  (Of  (  1 ,  •>)  *nc  (  51  *IJC(  ‘i)  «r)C(  4, 5)*DC.  (  “j )  ) 

T6  =  !  (6  )  /(L1C(  I  ,6)  «OC(  2,4)  «OC{  6)  <=nr.  (  4,6)  *OC(  5, 4)  ) 

=  Pl<‘lHl)<'^(2)<'fM?)<‘i')(4)*f'(5)*A(S)*(ri  +  T?fTi<-T4«-I5fT6) 

T1  ^  1.  (2)  /  (DC  (  3, 2  )  *OC(  4, 2)  *OC(  5,2  )  *00  (  6,  ?)  ) 

T2  =  E  (  3  )  /  (DC  (  2  ,  i)  *DC  (  4, 3)  *r>C(  5, -^  )  «nC(  6, 1)  ) 

T3  =  t (4 ) / (DC ( 2,4 )*DC ( 3,4) *DC ( 5,4) *nC( 6,4) ) 

T/,  =  fc  (  S  )  /  (  DC  (  2 , 5)  <‘DC(  3,5  )  «'0C(  4, 5  )  <‘()C(  6, 5)  ) 

T5  =  t  (  6  )  /  (DC  (  2, 6  )  *DC  (  3 ,6)  <'DC(  4,6  )  *0C  (  5 ,6  )  ) 

XN6  =  XN6  +  P2<'b  (  2  )  <=B(  3)  *B(4  )  ♦(3(  5)  «A(  ?  )  <=  (  T  1  ♦T2«-T  3i-T4*^T5) 

Tl  =  (  (  3)/(DC(4,3)*DC(5,3)*nC(6,3)  ) 

T2  =  h  (4)  /  (DC(  3,4)  «()€(  5,4)  *0C(6,4  )  ) 

T3  =  L(5)/(DC(3,5)<‘DC(4,5)*0C(6,5)  ) 
r4  =  h  (6  )  /  (  0C(  3,6  )  «DC(  4,6)  <=DC(  5, 6  )  ) 

XN6  =  Xl\i64-P3*(3(  3)  *8(4) -frBC  5  )*A(  3)  *(  Tl  ♦■T2«-T34-T4) 

Tl  =  h(4)/(DC(5,4)*DC(6,4)  ) 

T2  =  r ( 5 ) / (DC( 4,5 ) *DC( 6,5)  ) 

T3  =  r (6 ) / (DC ( 4,6) *DC ( 5,6)  ) 

XiN6  =  XN6  +  P4*H  (  4  )  *!3  (  5  )  *A  (4  )  *(  T  1<-T2 +T  3  ) 

Tl  =  (■(5)/DC(6,5) 

T2  =  C( 6) /DC (5,6) 

XN6  =  XN6+P5*3(  5  )  *A  (  5  )  *  (  T  I +T2  ) +P6*A(  6  )  *L- (  6) 

AB (  IX,  [ 1  )  =  AB(  I  X,  I  1  )  +XN6 

GO  TO  -19 

70  DO  71  K=l,NC 
I  L=l  (  K  > 

fa(K)  =f'y(LL) 

71  A(K)=AA(LI.  ) 

Pl=BC( 1 ) 

P2=8C( 2 ) 

P3=BC( 3) 

P4=BCI 4) 

P5=RC( 5) 

P6=BC ( 6 ) 

P  7  =  BC(  7) 

7  no  701  I  =  1  ,  7 

DC  701  J  =  1,7 
7Gl  0C( 1 ,J)  =  3( I )-B( J) 

DO  702  I  =  1,7 
7G2  E(  I )  =  CXPl-Bl  I  )*T) 

T 1  =  M I ) / ( DC ( 2, I ) *nc( 3,1) *DC (4,1) *DC (5,1) *0C (6,1) *DC (7,1)) 

T2  =  e  (2)/(DC(  l,2)*DC(  3 , 2 ) *0C ( 4 , 2 ) *DC ( 5 , 2 ) *0C ( 6, 2 ) *DC ( 7,2)  ) 

T  3  =  (:(  3)  /  (DC(  1 , 3)  *nC(  2 ,3)  *0C(4,  3  )  *nc  (  5 ,  3  )  *DC  (  6 , 3  )  *nc  (  7,3)) 

T4  =  t-  (4  )  /  (DC(  I  ,4)  *DC(  2,4)  *0C(  3,4)  *DC(  5,4  )*  DC  (  6,4  )*DC  (  7,4  )  ) 

T  5  =  '  ( 5 ) / ( DC (  1 , 5 ) *DC (2,5) *DC (3,5) *DC ( 4 , 5 ) *DC ( 6 , 5 ) *DC (  7,5)) 

T6  =  t(6) /(0C(  1 ,6)*DC( 2,6) *0C( 3 , 6 ) *nc ( 4 , 5 ) *0C ( 5 , 6 ) *DC ( 7,6)  ) 

T7  =  1  (  7) / (DC(  1 , 7)*0C( 2,7) *DC( 3,  7 )*DC( 4,  7 ) *DC ( 5 , 7 ) *DC ( 6 ,  7)  ) 

XN7  =  2  1*'3(  1  )*B(2)*B(3)*B(4)*B(5)*B(5)*A(  l)*(Tl<-T2+T34-T4  +  T5+T6«-T7) 
Tl  =  (  (2)  /  (DC  (  3, 2)  *DC(4,2)  *DC(  5,2)  *0C(6, 2)<'DC(  7, 2)  ) 

T2  =  L  (  3  )  /  (  DC  (  2,  3  )  *0C  (  4 , 3)  *DC  (  5,  3  )  *DC  (  6, 3)  *DC(  7,  3  )  ) 

T3  =  ^  (4 ) / (DC ( 2,4) *DC( 3 ,4) *DC( 5, 4) *0C( 6,4) *0C( 7,4) ) 

T4  =  I (5)/(DC( 2,5)*DC(3,5)*0C(4,5)*0C(6,5)*DC( 7,5) ) 
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=  e { 6  ) / (DC ( 2, 6) *DC ( 3»6) *00(4, 6) *0C( 5,6) *0C( 7,6)  ) 

T6  =  iH7)/(DC(2,7)<‘0C(3,7)*0C(4,7)«DC(5,7)4‘nC(6,7)  ) 

XN7  =  XN7  +  P2<‘n(2)*B(3)*B(4)*B(5)*R(6)«A(2)*(  T1  +  T2  +  T3  +  T4  +  T5  +  T6) 

T1  =  t  (3)/(nC(4,3)*0C(5,3)=e=0C(6,3)*nC(7,3)  ) 

T2  =  E (4 ) / (0C( 3,4)#0C( 5,4) 000(6,4) ♦0C( 7,4) ) 

T3  =  F(5)/(0C{3,5)O0C(4,5)*0C(6,5)*0C(7,5) ) 

T4  =  F (6) / {0C( 3,6) ODC( 4,6) ODC(5,6) *0C(7,6) ) 

T5  =  c(7)/(L)C(3,7)0DC(4,7)O0C(5,7)ODC(6,7)) 

XN7  =  Xi\7+P3OR(3)Ofi(4)0B(5)O8(6)0A(3)*ITl  +  T2*^T3<-T4  +  T5) 

Tl  =  F(4)/(0C(5,4)ODC(6,4)ODC(7,4)  ) 

T2  =  e(5)/(0C(4,5)*0C(6,5)ODC(7,5)  ) 

T3  =  E(6)/(0C(4, 6)000(5, 6)*0C(7, 6) > 

T4  =  F(7)/(D0(4, 7)000(5, 7)000(6, 7) ) 

XN7  =  XN7+P40B  (  4)  OB(  5)  OB(6  )OA(  4)  0(  T1  +  T2  +  T3-H4) 

Tl  =  F(5)/{0C(6, 5)000(7, 5)  ) 

T2  =  E(6)/(00(5, 6)000(7, 6) ) 

T3  =  e(7)/(D0(5, 7)000(6,7)  ) 

XN7  =  XN7  +  P50B(5)OB(6)oA(5)o(TH-T2  +  T3) 

Tl  =  L-(6)/00(7,6) 

T2  =  &( 7)/00(6,7) 

XM7  =  XN7+P608(6)0A(6)0(T1+T2) +P70A( 7)oE( 7) 

AB  (  I  X,  I  I  )  =A8(  IX,  I  I  )+XNi7 

GO  TO  99  I 


80  DO  81  K=1,N0 

LL=L(K) 

6 ( K) =8H ( LL  ) 

81  A(K)=AA(LL) 

Pl=ftO( 1 ) 

P2=B0( 2) 

P3^BL ( 3) 

P4=B0 ( 4 ) 

P5=BC ( 5 ) 

P6^ftC (6  ) 

6  7  =  2.0  (  1 ) 

PH=BO  (  .3  ) 

8  DO  801  I  =  1,8 

on  HOI  J  =  1,8 

801  nc( 1 , J )  =  8( I )-B( J) 


802 


00  802  I  =  1,8 

E(  1 )  =  EXP(-B(  1  )0T) 

T  1  =  t-  (  1 )  /  (00(  2,  1  )*nC  (3,  1)  000(4,  1  )  of)C  (5,  1)000(6, 1 )  ODC  (  7,1)  00C(  8,1) 
1  ) 

T2  =  [.  (  2)  /  (U0(  1 , 2)  O00(  3, 2)  OD0(  4, 2  )O0C(  5, 2)  O00(  6, 2)O00(  7,2  )  000(8,2  ) 
1  ) 

T  3  =  F  (  3  )  /  {i)0(  1,  3)000(2 , 3  )  OQO  (  4 , 3  )  OOC  (  5 , 3  )  *00  (  6 , 3  )  *00  (  7,3)O0C(  8,  3) 
1  ) 

T4  =  E ( 4) / (DC ( 1 ,4) O00( 2,4) OOC( 3,4) OOC( 5,4)O0C( 6,4)O0C( 7,4)OOC( 8,4) 
1  ) 

TS  =  r ( S) / (DC (1,5) 000(2,5) O0C( 3,5 ) OOC{ 4, 5 )O0C( 6, 5 ) *00 ( 7 , 5 ) ODC ( 8, 5 ) 
L  ) 

To  =  ■'  (5  )  /  (;)C(  I,  6)  0UC(  2,6)  *00  (  3,6  )  one  (4, 6)  one  (  5,6)  *00  (  7,6)  *00  (  8,6) 
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I ) 

T  7  =  f  (  7)  /  {  (  1 , 7)*r>C(  ?,  7)  *DC(  1,  7  )  *nc  (  4,  7  )  ♦OC  (  5, 7  )  *00  (  6,  7  )  *OC  (  8,7) 

I  ) 

To  ^  (-  (  8  )  /  (i.)C  (  I  ,  8)  *DC  (  2,3)  *nC(  3, 8)  *0C(  A,8  )*DC(  S,  8  )  *00  (  6, 8  )  *0C  (7,8) 
1  ) 

X.\8  ^  tU*LT{  I  )  *8(  2  )  *rt(  3  )*3(  4)*R(  8)*8(  6)*8{  7  )*A(  I  )*(  T  l  +  T?<-T3<-T4-»- T8<-T 
1 6  +  T  7  +  r  ) 

T 1  =  t  ( 2) / { OC ( 3, 2) *00 ( 4,2 ) *00 ( 5,2 ) *00 ( 6 , 2 ) *00 ( 7 , 2 ) *00 ( 8 , 2 ) ) 

T2  =  ►  (  n  /  (00  (  2, 3)  *1)0  (  4, 3)  *00(5,3)  *00  (  6 , 3  )  *00  (  7 , 3  )  *00  (  8 , 3  )  ) 

T3  =  E(4 ) / (n0( 2,4 )*D0 ( 3,4) *00( 5, 4) *DC( 6,4) *0C( 7,4)*n0( 8,4  )  ) 

T4  =  I.  (  5)  /  (0C(  2, 5)  *00(  3,5)  *00(4, 5)  *00(  6, 5)*0C(  7, 5)#00(  8,5)  ) 

T5  =  t  (6)  /  (f)0(  2,6)  *0C(  3,6)  *00(4,6  )  *00 (  5,6)  *00 (  7, 6)  *00  (8,6)) 

Tt.  =  i:  (  7  )  /  (00  (  2, 7)  *00(  3, 7)  *00(4, 7)  *DC(  5, 7)  *00  (6, 7)  *00  (  8, 7)  ) 

T7  =  ( rj )  /  (00  (  2,8)  *nO(  3,8)  *00  (4, 8)  *00  (  5,  B)  *00  (  6, 8)  *00  (  7,8)  ) 

X\;5  =  '<Ns  +  P2*(5(2)*B(3)*8(4)*R(5)*6(6)*8(7)*A(2)*(T  1  +  T  2+T  3  fTA  T  5+T  6 
H-T7) 

T1  =  i-(  3)/(OC(4,3)*OC(5,3)*DO(6,3)*nO(7,3)*OC(a,  3)  ) 

T2  =  t (4) / (00 ( 3,4 ) *00 ( 5 ,4) *0C( 6,4) *0C( 7,4)*0C( 8, 4) ) 

T  3  =  r-  (5  )  /  (0C(  3,5  )*00(4,5)  *00(6,  5)  *00(7,5)  *0C(  8,5)  ) 

T4  =  L (6) / (00 ( 3,6 )*0C( 4 ,6) *00( 5, 6) *00( 7,6)*DC( 8,6) ) 

T5  =  c ( 7) / (00 (  3, 7) *00( 4, 7) *00( 5, 7 )*n0 (6, 7)*D0( 8, 7)  ) 

T5  =  ;- (  8  )  /  (0G(  3,  8)  *00  (  4,8)  *D0(  5, 8)  *00(6,8  )*00{  7,8)  ) 

Xiv3  =  XiN.8+P3*8(  3)*B(4)*3(5)*R(6)*P(7)*A(3)*{T1+T2  +  T3+T4  +  T5  +  T6) 
ri  =  ■=  (4  )  /  {0C(  5,4  )  *D0(  6,4)  *00(  7,4)*0C(  8,4)  ) 

T2  =  t ( 5) / (00 (4, 5 ) *DC( 6,5) *OC( 7,5 ) *0C( 8, 5 ) ) 

T3  =  6 (6)/ (00(4,6) *00(5,6 )*0C( 7,6) *0C( 8,6) ) 

14  =  t(  7)/ ('10(4, 7)  *00  (5, 7)  *00  (6,  7)  *00  (8,7)  ) 

T5  =  t ( 8 ) / (00 ( 4,8) *00( 5, 8) *0C( 6, 8 ) *0C( 7,8 )  ) 

X.\B  =  X'>ic,  +  P4*8(  4)  *B(  5)  *rt  (  6  )*R(  7)  *A  (4  )  *(  T  l-t-T2  +  T3*T4  +  T5) 

11  =  h  (  5  )  /  (00(  6, 5)  *L!0.(  7,5)  *0C(  8 ,5  )  ) 

12  =  l(6)/(00.(5,6)*00(  7,6)*0C(8,6)  ) 

13  =  1- (  7  )  /  ((30  (  5, 7)  *00  (6 ,7)  *D0(  8,  7)  ) 

14  =  e  ( 8 ) / (0C( 5 ,8 )*OC ( 6 ,H) *00( 7, 8 )  ) 

X'JS  =  XNrt  +  P5*8(  5)*8(6)  *B(  7  )*A(  5)  *(  ll+12«-13  +  14) 

T  1  =  L(6)/(0r.  (  7,6)  *00(8,6)  ) 

1?  =  P ( 7 ) / (00(6, 7) *0C{ 8, 7) ) 

13  =  t (8 ) / (00 (  6, 8) *00(  7 , 8)  ) 

X'!8  =  Xrvi8  +  P6*8  (6)  *8  (  7)  *A  (6  )  *(  1 1 +  12  +  T3) 

11  =  c(  7)/0C(8,7) 

12  =  P(8)/00(7,8) 

=  <N.3  +  P /*B  (  7  )  *A  (  7  )  *  {  11 +12) +P8*A(  8  )  *F  (  6  ) 

A  1  (  1  A,  I  1  )=A(-,  (  I  X,  I  I  )  +XN8 

G.;  10 

90  i;i:  >-1  -<=1,60 

L  t.  =L  (  K  ) 

:.'(n)=2H(LL  ) 

Sr  1  A  (  .N  >  =  4  A  (  L  L  ) 

PL=PO(  1) 

P2=r'0  (  2  ) 

P  3  =  "•  C  (  3  ) 

P4=0C ( s  ) 
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^3) 

PA  =  HC  (  ') ) 

P  7  =  BC (  7  ) 

PH=BC ( H ) 

P  =  H  C  (  <^3  ) 

^3  DO  901  I  =  1,9 
DO  9Cl  J  =  1,9 

901  CC(  I  ,  J)  =  15  (  I  )-H(  J) 
no  90?  I  =  1,9 
90?  fctn  =  exp(-H(i)*r) 

T  l  =  i  (  1  )  /  (  DC(  ?,  I  )  *0C(  3,  1  )  *0C  (  4,  I  )  <<00  (  9,  1  )*0C(  6, 1  )  *00  (7,1)  *0C  (  B,  I  ) 
i*nc( 9, 1 ) ) 

T?  =  p  (  ?)  /  (DC(  1,?)*DC(  3,?)  (  4,?)  <-|)C(S,?)*nC(6,  ?)  V^DCI  7,2)*0C(  «,?) 

1<'0C(9,?)  ) 

T  3  =  i:  (  3)  /  (0C(  1,  3)*0C(  2,  3)  (4,  3)  *0C(  9,3)«DC  (6, 3)4DC{  7,  OX-nC  (  fi,  3) 

l*DC(9, 3)  ) 

T4  =  F(4)/(DC(l,4)  4UC  (  ?,  4)  *0C(  3 , 4  )  <=00  (  5,4)  '['DC  (  6,4)  *0C  (  7,4)*nC(8,4) 
l«DC (9,4) ) 

T5  =  t  (  5  )  /  (DC  (  1 , 5  )  ’*=!)C(  2,5)  *00  (3,5)  *00(4,5  )  *00  (  6,5)  *00  (7,5)  *00  (8,5) 
l*OC( 9,5)) 

T6  =  t ( 6 ) / (00 (  1 , 6 ) *00  (  ? ,6) *00 ( 3, 6) *00 ( 4,6) *00 (5,6)  *00  (  7,6 ) *00 ( 8,6) 
1*00(9,6)) 

T  7  =  E ( 7 ) / (00 (  1 , 7) *00 (2,7) *00 ( 3,  7) *00 ( 4 , 7 ) *00 ( 5 , 7 ) *00(6,  7 ) *00 (  8,  7) 
1*00(9,7) ) 

Ta  =  h(8)/(nC( I, 8 ) *00 ( 2, a ) *00 ( 3, 8) *00 (4, 3) *00 ( 5, 8 ) *00 ( 6 , 8 ) *00 ( 7,8) 
1*00(9,8) ) 

T9  =  E (9) / (00 ( 1,9) *00 ( ?,9) *DC( 3 , 9 ) *00 ( 4 , 9 ) *00 ( 5 , 9 ) *00 ( 6 , 9 ) *00 (7,8) 
1*00(9,9)) 

X\9  =  Pl*l3(l)*B(2)*R(3)*ri(4)*H(5)*a(6)*a(  7)*R(8)*A(l)*(Tl  +  T?*T3*T4 
l+T5+T6*T7f TH+T9) 

n  =  F(?)/(0C( 3,2)*OC(4,2)*DC(5,2)«OC(6,?)*OC( 7 , 2 ) *00 ( 8 , 2 ) *00 ( 9, ? ) 
I ) 

T?  ^  E(3)/(00(2,3)*00(4,3)*OC(5,3)*00(6,3)*OC( 7 , 3 ) *0C ( 8 , 3 ) *00 ( 9, 3 ) 
L  ) 

T3  =  F (4) / (00( 2,4) *00(  3,4) *0C( 5, 4) *00( 6,4)*0C( 7 , 4 ) *00 ( 8 , 4 ) *00 (9,4) 
1  ) 

T4  =  fc(5)/{OC(2,5)*OC(3,5) *00 ( 4 , 5 ) *0 0 ( 6 , 5 ) *00 ( 7, 5 ) *0C( 8, 5 ) *00 (9,5) 
I  ) 

T5  =  0(6 ) /(0C( 2,6) *00 ( 3,6) *00 ( 4,6) *00 ( 5,6) *00 ( 7,6) *00 ( 8,6)*0C( 9,6 ) 
I  ) 

T6  =  F ( 7) /(0C( 2,  7) *00(  3 , 7  )  *00  (  4 , 7 ) *00 ( 5 , 7 ) *00 ( 6 , 7 ) *00 ( 8, 7)  *00 (  9,  7) 
1  ) 

T7  =  f:  (  8  )  /  (00  (  2 ,8)  *0C(  3,8)  *00  (  4,8)  *00  (  5 , 8  )  *O0  (  6 , 8  )  *00  (  7, 8)  *00  (9,8) 
I  ) 

T8  =  F(9)/(OC(2,9)*DO(3,9) *00 (4,9) *00 ( 5 , 9 ) *00 ( 6 , 9 ) *00 ( 7,9) *00 ( 8,9 ) 
1  ) 

XND  =  XN9+P2*P( ?)*R( 3) *8(4 )*e(5)*P(6)*P(7)*3(8)*A(2)*(Tl*T2*T3vT4* 
IT55-T6+T7+T8) 

T1  =  F(3)/(O0(4,3)*UC(5,3)*D0(6,3)*DO(7,3)*OC(8,3)*DC(9,3) ) 

T?  =  F (4 ) / (DC( 3, h) *00( 5,4) *00(6,4) *00( 7,4)*0C( 8,4) *00( 9,4) ) 

T3  =  F(5)/(00(3,5)*00(4,5)*DC(6,5)*OC(7,5)*OC(8,5)*OC(9,5) ) 
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X4  =  F  (  6  )  /  (DC  (  3  t  6  )  *DC  (  4 ,6)  *0C  ( *>  t  6  )  ♦l)C(  7,  <S  )  *0C(  8, 6  )  *DC(  9, 6  )  ) 

T5  =  t ( 7 ) / ( OC ( 3» 7 ) *0C (  4 ♦ 7 ) *00 ( 5 , 7 ) ^OC I  6 1 7) *0C( 8 f 7 ) *0C ( 9, 7  )  ) 

T6  =  F  (  8  )  /  (DC  (  3t  8)  ♦DC  (  4  f  8)  *DC(  5,  8  )  *nC(6 , 8  )  *0C(  7, 8  )  ♦OCi  9, 8  )  ) 

T7  =  F  (  9)  /  (DC  (  3,9  )  «0C  (  4, 9)  <=0C  (  8, 9  )  ♦DC  (  6, 9  )  ♦DC(  7 ,9  )  ^OC  (  8 , 9  )  ) 

Xr\(9  =  XNl9  +  83*8(  })^B(4)*B(5)*B(6)*B(7)*B(8)*A(3)^(Tl+-T2+T3  +  T4-»-T  8+  T  6 

1  +  77) 

Tl  =  F (4) / (DC( 5,4) ♦0C(  6, 4) ♦0C( 7,4 ) ♦0C( 8,4)*0C( 9, 4) ) 

T2  =  t(5) /(DC(4,5)*0C(6,5)*'0C(7,5)<'0C(8,5)<=0C(9,5)  ) 

73  =  t (6 ) / (OC ( 4,6) ♦0C( 5 ,6) ♦0C( 7,6 ) ♦DC ( 8,6 ) ♦DCC 9,6) ) 

74  =  fc(7)/(CC(4,7)^OC(5,7)^OC(6,7)^DC(8,7)*OC(9,7) ) 

75  =  F  (  8  )  /  (  DC(  4 , 8)  ♦DC  (  5 , 8  )  ♦OCC  6,  8  )  ♦DC  (  7, 8  )  ♦0C(  9,  8)  ) 

76  =  E(9)/(UC(4,9)*0C(5,9)«0C(6,9)^0C(7,9)^0C(8,9) ) 

XN9  =  XN9+P4^B (4)^8(5)*B(6)^8(7)^B(8)*A(4)^(7H-T2+73+74+75+76) 

71  =  t(5)/(DC(6,5)*0C(7,5)«DC(8,5)*DC(9,5) ) 

72  =  £ (6 ) / (DC ( 5 ,6 ) ♦DC( 7,6) *DC( 8,6) ♦DC) 9, 6) ) 

73  =  e  (  7) /(DC ( 5, 7) ♦DC (6 , 7) ♦DC ( 8, 7) ♦0C(9, 7) ) 

74  =  E(8)/(DC(5,8)^0C(6,8)*0C(7,8)^0C(9,8) ) 

75  =  £  (  9 ) / (DC  (  5, 9) ♦DC ( 6, 9) ♦DC{ 7, 9 ) ♦DC) 8 ,9 ) ) 

XN9  =  XN9  +  P5«B(5)^B(6)*B(7)^B(a)^A(5)*(  7 1 +  7  2<-73+T4  +  T  5  ) 

71  =  F ( 6 ) / (DC ( 7,6 ) ♦DC) 8 ,6) ♦DCM  9,6)  ) 

72  =  fc(7)/(DC(6,7)*DC(8,7)^DC(9,7)  ) 

73  =  F  (  8  )  /  (DC  (  6,  8)  ♦DC)  7, 8)  ♦DC(  9, 8  )  ) 

74  =  F(9)/{0C(6,9)^DC( 7,9)^DC(8,9) ) 

XN9  =  XN9  +  P6*B ( 6 ) *8 ( 7 ) ♦B ( 8 ) *A ( 6 ) ♦ ( 7  I +72+ 73+74) 

71  =  E( 7)/(0C(8,7)^DC(9,7) ) 

72  =  F(8)/(DC( 7,8)*CC(9,8)  ) 

73  =  E (9) /(DC( 7,9)^0C(8,9) ) 

XN9  =  XN9+P7+B(7)*B(a)*A(7)+(7l+72+73) 


71  =  E(8)/DC(9,8) 

72  =  E(9)/DC(8,9) 

XN9  =  XN9+P8«B(8)*A(8)*( 71+72) +P9^A(9)+E(9) 

AB(IX, I  I)=A6( IX, II  )+XN9 

GC  70  99 

110  DO  111  K=1,NC 
LL=L(K) 

B(K)=6R(LL) 

111  A(K)=AA(LL) 

P1=BC(  1) 

P2=BC( 2) 

P3=BC( 3) 

P4=BC(4) 

P5=BC(5) 

P6=BC(6) 


P7=BC( 7) 
P8=BC(8) 
P9=BC(9) 
P10=BC( 10) 
1000  DO  1001  I 
DO  lOCl  J 
1001  0C(1,J)  = 
DO  1002  I 


=  1,10 
=  1,10 
B( I )-e( J) 
=  1,10 
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100?  r- ( I )  =  EXP(-R(  n-s'T ) 

Tl  =  F(l)/(nC(2,l)<'0C(3,l)  *00(4,1  )  *0C(  5,  I  )  *nc(  6,  1  )  *0C(  7,1)  *0C(  «,  1  ) 
1*0CI9,1 )*0C(  10,1)  ) 

T?  =  h{?)/(0C(  1, ?) ♦DC( 3,2) *00(4, 2 ) *DC ( 5 , ? ) *0C ( 6 , 2 ) *00 ( 7,2 )*00( 8,2) 
1*00( 9,2)*00(  10,2)  ) 

7  3  =  t ( 3  )  / (DO (  1 , 3) *nO( 2, 3) *00(4, 3) *00( 5,  3)*00( 6, 3 ) *00 ( 7,3) *00(8,  3 ) 
1*00(  9,  3)*00( 10,3) ) 

T4  =  L (4 ) / (00 (  1 ,4) *00 ( 2 ,4) *00 ( 3 , 4 ) *00 ( 5 , 4 ) *00 ( 4 , 4 ) *00 ( 7,4)*00( 8,4) 
1*00(9, 4)*00(  10,4)  ) 

T5  =  F (4) / (00(  1,51*00(2,5) *00 ( 3, 5) *00 (4,5) *00 ( 6,5) *00 ( 7, 5) *00 ( 8,5) 
1*00(9, 5)*DC(  10,5)  ) 

06  =  F (6) / (00( I ,6) *00 ( 2,6) *00 ( 3 , 6 ) *00 ( 4 , 6 ) *00 ( 5, 6 ) *00 ( 7,6) *00 ( 8,6) 
1*00(9,61*00(10,6) ) 

T7  =  C(7)/(00(  1,  7) *00 (2,7) *00 { 3,71*00(4,  7) *00 ( 5,  7) *00 (6,7) *00 ( 8, 7) 
1*00(9, 7)*00(  10,7)  ) 

T8  =  E ( 8)  / (00(  1, 8) *00 ( 2 ,8) *00 ( 3,8 ) *00 ( 4 , 8 ) *00 ( 5 , 8 ) *00( 6 , 8 ) *00 ( 7,8) 
1*00(9, 8)*00( 10,8) ) 

T9  =  £(9)/{00(  1,9)*DO(2,9)*00(3,9)*00(4,9)»00( 5 , 9 ) *00 ( 6, 9 ) *00 ( 7,9) 
1*00( 9,9)*00( 10,9) ) 

TIO  =  b(10)/(00(l,  10)*00( 2, 10) *00 (  3,  1 0 ) *00 ( 4 ,  1 0 ) *00 (  5 , 1 0) *00 ( 6 , 10 ) 
1*00(7,  10) *00 (8, 10) *00 (9, 10) ) 

XNIO  =  P1*H(  1)*B(  2)*B(  3)*R(4)*B(  5)*B(6)*(3(  7)*B(8)*B(9)*A(  1  )*(  TUT? 
1  +  T3+T4  +  T5+T6+T  7+T8  +  T9+T10) 

Tl  =  E(2)/(00( 3,2)*00(4,2)*OC(5,2)*OC(6,2)*00( 7,2)*00( 8,2)*00( 9,2) 
1*00( 10,2) ) 

T2  =  E(3)/(00(2,3)*00(4,3)*OC(5,3)*00(6,3)*00( 7 , 3 ) *00 ( 8 , 3 ) *00 ( 9, 3 ) 
1*D0( 10,3) ) 

T3  =  E (4) / (00( 2,4)*00( 3 , 4 ) *00 ( 5 , 4 ) *00 ( 6 , 4 ) *00 ( 7 , 4 ) *00 ( 8 , 4 ) *00 ( 9, 4 ) 
1*00( 10,4) ) 

T4  =  F(5)/(0C(2,5)*O0(3,5)*O0(4,5)*OC(6,5)*OC( 7 , 5 ) *DC ( 8 , 5 ) *00 ( 9 , 5 ) 
1*001 10,5)) 

T5  =  E  (6)/ (00 (2, 6) *00(3,6) *00(4, 6) *00(5, 6) *00 (7, 6) *00(  8, 6) *00 (9, 6) 
1*00 ( 10,6 ) ) 

T6  =  E 1 7) /(00(  2, 7) *001  3,7) *00(4,  7) *00 ( 5,  7) *00 ( 6,  7) *00(  8,71*00(9,  7) 
1*D0(  10,7)) 

T7  =  E (8)/ (00(2,8 ) *00 ( 3, 8) *00(4, 8) *001 5,8) *00(6,8) *00 ( 7,8 )*00( 9,8) 
1*00( 10,8) ) 

T8  =  E(9)/(00(?,9)*0C(3,9) *00 (4,9)*00(5,9)*00(6,9)*00(7,9)*00(8,9) 
1*00(10,9) ) 

T9  =  E(  10) /(DO (2, 10) *00(3,  10 ) *00 ( 4 , 1 0 ) *00 ( 5 , 10) *00 (6,  10)*00( 7,  10)* 
100(8, 10)*00(9, 10) ) 

XNIO  =  XN10*P?*B(2)*B(3)*R(4)*B(5)*B(6)*0(7)*R(8)*B(9)*A(2)*(T1*T2 
1  +  T3*  r4*T5*T6  +  T7+T8*T9) 

T  1  =  E (  3 ) / ( 00 ( 4 , 3 ) *00 (5,3) *00 (6,3) *00 ( 7 , 3 ) *00 ( 8 , 3 ) *00 (9,3) *00 (10,3 
1)  ) 

T2  =  E(4)/ (D0(  3,4) *00(  5, 4) *00  (  6,4 ) *00 (  7,4)*O0( 8 , 4 ) *00 ( 9, 4 ) *00 ( 1 0, 4 
1  )  ) 

T3  =  F (5) /(00( 3,5)*00( 4,5)*nO(6, 5) *00( 7,5) *00( 8 , 5 ) *00 ( 9, 5 ) *00 ( 10,5 
1 )  ) 

T4  =  E(6)/(00(3,6’*0C(4,6)*D0(5,6)*D0{7,6)*00(8,6)*00(9,6)*00(10,6 
1  )  ) 
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T5  =  E(7)/(0C(3,7)*DC(A,7)*0C(5,7)*0CI6,7)*DCI 8, 7) *00 ( 9, 7 1 ♦DC ( 10,7 
1)  ) 

T6  =  E(8)/(DC( 3,8)*DC{A,8)*0CC5,8)*0C(6,8)*DC(7,8)*DC(9,8»*DC( 10,8 
1)  ) 

T7  =  E(9)/(0C( 3,9)*0C(4,9)*0C( 5,9)*DC(6,9»*0C( 7,9>*0C<  8,9)*0CI 10,9 
1)  ) 

T8  =  E ( 10) /(0C(3, 10) *00 (4, 10)*0C( 5, 10) *00(6, 10 )*DC( 7, 10)*0C( 8, 10)* 
1DC(9, 10) ) 

XNIO  =  XN10+P3*B(3)*B(4)*8(5)*B(6)*B(7)*8(8)*B(9)*A(3)*(Tl+T2+T3+T 
14+T5+T6+T7+T8) 

T1  =  E(4)/(DC(5,4)*0C(6,4)*0C(7,4)*0C(8,4)*0C(9,4)*0CC 10,4) ) 

T2  =  E(5)/(DC(4,5)*0C(6,5)*0C(7,5)*0CI8,5)*0C(9,5)*0C( 10,5) ) 

T3  =  E(6)/(0C(4,6)*DC(5,6)*0CI7,6)*0C(8,6)*0C(9,6)*0C( 10,6) ) 

T4  =  E(7)/(0C(4,7)*0C(5,7)*0C(6,7)*0C(8,7)*0C(9,7)*0C(10,7) ) 

T5  =  E(8)/(OC(4,8)*OC(5,8)*OC(6,8)*OC(7,8)*OC(9,8)*OC( 10,8) ) 

T6  =  E(9)/(0C(4,9)«0C(5,9)*0C(6,9)*0C(7,9)*0C(8,9)*0C( 10,9)) 

T7  =  E( 10) /(OC (4, 10) *00 (5, 10 ) *00 ( 6 , 10 ) *00 ( 7, 10 ) *DC ( 8 , 10 ) *00 ( 9, 10) ) 
XNIO  =  XN10+P4*B(4)*B(5)*B(6)*B(7)*B(8)*8(9)«A(4)*(T1+T2+T3+T4+T5+ 
1T6+T7) 

Tl  =  E(5)/(OC(6,5)*OC(7,5)*DC(8,5)*DCC9,5)*OC( 10,5) ) 

T2  E(6)/(DC(5,6)*0C(7,6)*DC<8,6)*DC(9,6)*0C(  10,6)  ) 

T3  =  E(7)/(0C<5,7)*0C(6,7)*0C(8,7)*0C(9,7)*0C( 10,7) ) 

T4  =  E(8)/I0C(5,8)*0C(6,8)*0C(7,8)*0C(9,8)*DC( 10,8) ) 

T5  =  E(9)/(DC(5,9)*0C(6,9)*0C(7,9)*0C(8,9)*0Cn0,9)  ) 

T6  =  E( 10) /(0C(5, 10) *00(6, 10)*0C( 7, 10)*0CI 8, 10)*0C(9, 10) ) 

XNIO  =  XN10+P5*8 (5)*B(6)*B(7)*8(8)*B{9)*A(5)*(T1+T2+T3+T4+T5+T6) 

Tl  =  E(6)/(DC(7,6)*0C(8,6)*DC(9,6)*0C(  10,6) ) 

T2  =  E(7)/(0C(6,7)*0C(8,7)*0C(9,7)»DC( 10,7) ) 

T3  =  E(8) /(0C(6,8)*0C(7,8)*0C(9,8)*DC( 10,8) ) 

T4  =  E(9)/(0C(6,9)*0C(7,9)*DC(8,9)*0C( 10,9) ) 

T5  =  E(10)/(OC(6,lO)*OC(7,lO)*OC(8,lO)*OC<9,lO)) 

XNIO  =  XN10+P6*B(6)*B(7)*B(8)*B(9)*A(6)«(T1+T2+T3+T4+T5) 

Tl  =  E(7)/(DCI8,7)*0C(9,7)*0C( 10,7) ) 

T2  =  E(8)/(0C(7,8)*0C(9,8)*0C(10,8) ) 

T3  =  E(9)/(0C{ 7,9)*0C(8,9)*DC( 10,9) ) 

T4  =  E(10)/CDC(7,10)«DC(8,10)*0C(9,10) ) 

XNIO  =  XN10  +  P7»B(7)*B(8)*8(9)*A(7)*(T1  +  T2+T3*^T4) 

Tl  =  E{8)/(0C(9,8)*0C(10,8) ) 

T2  =  E(9)/(0C(8,9)*0C(10,9)) 

T3  =  E( 10)/(0C(8,10)*0C(9,10)) 

XNIO  =  XN10+P8*B(8)*B(9)*A(8)*(T1+T2+T3) 

Tl  =  E(9)/0C( 10,9) 

T2  =  E( 10) /DC(9, 10) 

XNIO  =  XN10+P9*8(9)*A(9)*(T1+T2)+P10*A( 10)*E(10) 

A8( IX, I  I )  =  A8( IX,1 I )+XN10 
99  CONTINUE 

DO  96  M=l,NI 
MMM=MM(M) 

00  96  N=1,NT 
DO  96  1=1, MPM 

96  SUNIN, I )=SUN(N, I )+FUNl ( I , M) ♦ AB ( M , N ) ♦BI M ) 


171 


GNE/Phys/65-3 


97  CONTINUE 
2000  CALL  PRINT 
STOP 

990  F0RMAT(I2,E8.3) 

991  FORMAT (8E10. 4) 

992  F0RMAT(F5.0,15) 

993  F0RMAT(I3,10(F4.3,I3)) 
END 
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39 

4C 

41 

42 

43 

44 

144 

244 

344 

444 

145 


36 


35 

37 

47 

147 

247 

347 

447 

14B 


45 


49 

46 


SUBROUTINE  PRIM 

CCViMON  EVE  (20)  ,  SUN  (20,210)  ,  I  IF  ,  SUM  (  2  10  ) 
DIMENSION  ENERGY(210) 

DO  43  WO^I ,200 
IF(MC-1)40,39,40 
ENERGY(Mn)=().0 
GO  TO  43 

IF(M0-101 ) 41 ,41 ,42 

ENERGY  (MO)  =ENERG  Y  (  MO- 1  )  ^-.  0  I 

GO  TO  43 

ENERGY  (MO)  =ENERGY(  MG-I  )  1 

CONTINUE 

KAR  =  0 

GO  TO  ( 144,244,344,444) , ITF 

PRINT  155 

GO  TO  145 

PRINT  255 

GO  TO  145 

PRINT  355 

GO  TO  145 

PRINT  455 

PRINT  57, ( EVE( KP) ,KC=1 ,7) 

PRINT  53 
PRINT  65 
Mr3ST=0 

IF(KArN)  37,  36,37 
DO  37  KAR=1,?00 

PRINT  bl , ( SUN( MG.KAR ) , M0= 1 ,7) , ENERGY ( KAR ) 

MCST=Mr:sr^i 

1F(KAk-200) 35,47,35 

IF(M0ST-50) 37,44,44 

Cf^MT  INUE 

GO  TC  ( 14 7,247,347,447) , ITF 

PR  INI  155 

GO  TO  148 

PRINT  255 

GO  TO  148 

Print  355 

GO  T(]  14  8 

PRINT  455 

PRINT  59, ( E Vt ( RC )  ,K0  =  8 , I  3) 

PRINT  60 
PRINT  65 
MOST  =0 

IF(KAR-2J;))46,45,46 
DO  40  rAR=  1,200 

PRINT  62, ( SoN(MO,KAR ),M0=8,13) ,FNERGY(KAR) 

MCST=Mr]STf  1 

IF ( K AR-2U0 ) 49,54,49 

IF  (M[:sr-50)46,4  7,4  7 

CONT INUE 
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54  I^KITE  TAPE  3  ,  (  ENE  RGY  <  MO )  ,  M0=  1 ,  195  ) 

DO  70  M0=l,13 

70  WRITE  TAPE  3 , ( SUN ( MG , NI T ) , N I T=l t 1 95 ) 

END  FILE  3 
REWIND  3 
REWIND  1 
RETURN 

57  FORMAT  ( /, IX,5HTIME=,F5. 1, 6F15. 1 ) 

58  FORMAT  (/t6X,7HBIN  (1),8X,7HBIN  (2)f8X,7HBIN  (3»t8X,7HBIN  I4),8X,7 

IHBIN  (5>,8X,7HBIN  I6),8Xt7H8IN  17),8X,1HE) 

59  FORMAT  ( / , I X , 5HT I ME  = , F6 . 1 , 5F 15 . 1 ) 

60  FORMAT  (/,6X,7HBIN  (8),8X,7HBIN  (9),8X,8HBIN  ( 10 > , 7X, 8HBIN  C11»,7X 
l,8H8IN  ( 12 ) t7X,8HBIN  (13),7X,1HEI 

61  FORMAT  ( 1P7E15.5,0PF9.2) 

62  FORMAT  I 1P6E 15 . 5 » 0PF9. 2 ) 

65  FORMAT  (/) 

155  FORMAT! IHl ,24X,67HINSTANTANE0US  BETA  RAY  SPECTRA  FROM  THERMAL  PU-2 
139  FISSION  PRODUCTSt/ t29Xf 58HAT  VARIOUS  TIMES  AFTER  DETONATION  (BE 
2TAS/ENERGY  RANGE/SEO) 

255  FORMAT  (IH I , 24X, 66H I NST ANT ANEOUS  BETA  RAY  SPECTRA  FROM  THERMAL  U-2 
135  FISSION  PRODUCTS, /f29X,58HAT  VARIOUS  TIMES  AFTER  DETONATION  (BE 
2TAS/ENERGY  RANGE/SEO) 

355  FORMAT  ( IHl ,24X,66HINSTANTANE0US  BETA  RAY  SPECTRA  FROM  THERMONUCLE 
IR  FISSION  PRODUCTS,/, 14X,58HAT  VARIOUS  TIMES  AFTER  DETONATION  (BET 
2AS/ENERGY  RANGE/SEO) 

455  FORMAT! IHl ,20X,75HINSTANTANE0US  BETA  RAY  SPECTRA  FROM  U-235  FISSIO 
IN  SPECTRUM  FISSION  PRODUCTS ,/, 29X , 58HAT  VARIOUS  TIMES  AFTER  DETONA 
2TI0N  !  BETAS/ENERGY  RANGE/SEO) 

END 
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SliHROUTINE  k^INDER 

1  READ  TAPE  I, AND, (2) 
GC  TO  1 

2  RETURN 
END 
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PROGRAM  IV 

Main  Program 

A,AA 

Concentration  of  isotope  immediately  after  fission 

AB 

Concentration  of  isotope  at  the  time  of  interest 

ANO 

Atomic  weight  of  mass  chain 

B,  BB 

Decay  constant  of  isotope 

BC 

Branching  indicator  which  shows  the  fraction  of  the  isotope 

which  contributes  to  the  decay 

EVE,  T 

The  time  of  interest 

FUN  1 

Value  of  the  normalized  spectral  shape 

ITF 

Indicator  of  the  type  of  fission:  1  indicates  Pu-239 

thermal  neutron  fission;  2  indicates  U-235  thermal  neutron 

fission;  3  indicates  thermonuclear  neutron  fission; 

4  indicates  U-235  fission  spectrum  neutron*  fission 

L,  LL 

Indicators  which  show  which  isotopes  in  the  chain  con¬ 
tribute  to  the  production  of  another  isotope 

MI,  MM 

Number  of  energy  increments 

NC 

Number  of  isotopes  in  the  decay  path 

NI 

Number  of  isotopes  in  the  mass  chain 

NT 

Number  of  times  of  interest 

SUN 

Value  of  the  spectrum  at  the  time  of  interest 

YM 

Number  of  megatons  of  fission 

SUM 

Dummy  variable 

DC,E,Pi,F2, 


F3,P4,Pf>,P6, 

P7,P8,P9, 
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PiO,Tl,T2, 

T3,T4,T5,T6, 

T7,T8,T9, 

T10,XN1,XN2, 

XN3,XN4,XN5, 

Xn6,XN7,XN8, 


XN9,XN10 

=  Interim  Calculation  Variables 

I, II, IX, 

ITT,J,JJ, 

K,KK,K0,LX, 

M,MMM,N,NX 

=  Indicies  of  "DO"  loops 

Subroutine  PRINT 

PRINT 

=  Subroutine  to  print  results 

ENERGY 

=  Value  of  energy 

MOST 

=  Interim  Calculation  Variable 

K0,KAR, 

MO, NIT 

=  Indicies  of  "DO"  loops 

Subroutine  WINDER 

WINDER  =  Subroutine  to  advance  tape  to  read  the  next  record 

ANO  =  Dummy  variable 
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PROGRAM  IV 


Fir.  20 

?lov<-  Claart  for  Proc;ra:n  IV 
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